Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,125 C→T 5.4% T162T (ACC→ACT)  ycdY → redox enzyme maturation protein (REMP) chaperone for YcdX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,1250CT 5.4% 737.1 / 21.4 220T162T (ACC→ACT) ycdYredox enzyme maturation protein (REMP) chaperone for YcdX
Reads supporting (aligned to +/- strand):  ref base C (95/113);  new base T (6/6);  total (101/119)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 220 total aligned reads displayed.

GTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCC  >  NC_000913/1099862‑1100368
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
 tgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGACGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGACGCTTAGAGCGAAGCACTGGAAACCCTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAACGTGGGGGGCCATTCAACCACGCCTTTCTGGCCCACCC‑TTGCAccccca                                                                                                                                                                                                                                         >  4:241992/1‑271 (MQ=255)
                   cATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCCGTGCGTATCCGTTCCCCTGGTGTGCTGCGTTCCTTTGCAACGTGGAGGCCCCTGCAACCCCGCCTTTCTGGCGCACCA‑TGCCACCGCTAACCCGCGATCCCCtt                                                                                                                                                                                                                          >  4:129352/1‑272 (MQ=255)
                   cATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGt                                                                                                                                                                                                                       >  2:192639/1‑275 (MQ=255)
                     tGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTTGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCCCCA‑TGGCCCCGCAAACCCGCGTAGCCATTAGTg                                                                                                                                                                                                                      >  4:320132/1‑274 (MQ=255)
                         ttggggggggcccggaagcggagaCGCGCCCCTGTCTTTCTGAGGGGTGGATGCCATTAGCGGTTAGGCCAGCCGATCACATCGGCACATTGCTTCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACAGTTCAGGGAGTATCTGTTACCCTAGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAAtg                                                                                                                                                                                                                 <  4:347522/263‑1 (MQ=255)
                                                             tttctGGGGTAGGGATGCAATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCCCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgt                                                                                                                                                                              <  4:532945/274‑1 (MQ=255)
                                                                              ccATTAGCGGATAAGGAGGCCGATACAATGGGGAAATTGGTGCTCGCAGCTTCCTGGCGGGAAGGTCAGGAAACGGAAGATGAGGGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCAt                                                                                                                                                            <  4:7488/275‑1 (MQ=255)
                                                                                      gggtaagCCAGCCGATGACACCGGCCCATTGCTGCTCGCAGCTTCCTAGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCCCTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATACAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCa                                                                                                                                                     <  4:411479/268‑1 (MQ=255)
                                                                                                      cacaTCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgt                                                                                                                                     <  1:2072/274‑1 (MQ=255)
                                                                                                      cacaTAGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgt                                                                                                                                     <  2:484759/274‑1 (MQ=255)
                                                                                                                      tgctCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCt                                                                                                                    <  2:320133/275‑1 (MQ=255)
                                                                                                                       gctCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAACTGTTTGCt                                                                                                                    >  2:398752/1‑274 (MQ=255)
                                                                                                                            cAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCAcc                                                                                                                                                                                                                                             >  5:957064/1‑148 (MQ=255)
                                                                                                                             aGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCg                                                                                                                                                                                                                                            <  6:887567/148‑1 (MQ=255)
                                                                                                                                    tggctggAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAAccc                                                                                                                                                                                                                                     >  5:677267/1‑148 (MQ=255)
                                                                                                                                    tagcagaaaacgcggccatgaatgagatcgagcacacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACCA‑TGGCACCGCTAAccc                                                                                                                                                                                                                                     <  5:145924‑M1/21‑1 (MQ=255)
                                                                                                                                     ggctggAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCcg                                                                                                                                                                                                                                    >  8:562754/1‑148 (MQ=255)
                                                                                                                                           aaGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGc                                                                                                                                                                                                                              >  6:602403/1‑148 (MQ=255)
                                                                                                                                               tCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          <  7:960116/148‑1 (MQ=255)
                                                                                                                                                  gTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCa                                                                                                                                                                                                                                                     >  8:305883/1‑119 (MQ=255)
                                                                                                                                                    cAACGGAAGATGAGAGCGAAGCACTGGAAACAGTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                                     >  7:721171/1‑148 (MQ=255)
                                                                                                                                                    cAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                                     >  5:232968/1‑148 (MQ=255)
                                                                                                                                                       cGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCAGCGATGCCATTAgatcgga                                                                                                                                                                                                                  >  5:376569/1‑142 (MQ=255)
                                                                                                                                                            gATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAATGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgcgc                                                                              <  2:115577/275‑1 (MQ=255)
                                                                                                                                                              atcgagcacacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATg                                                                                                                                                                                                                               <  5:668743‑M1/28‑1 (MQ=255)
                                                                                                                                                               gagaGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTggg                                                                                                                                                                                                             <  5:398412/145‑1 (MQ=255)
                                                                                                                                                                     gAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAAt                                                                                                                                                                                                                  <  2:565454/134‑1 (MQ=255)
                                                                                                                                                                      cacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCATTAg                                                                                                                                                                                                                        >  8:1129154‑M1/94‑128 (MQ=255)
                                                                                                                                                                      cacattttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCATTAg                                                                                                                                                                                                                        <  7:1129154‑M1/35‑1 (MQ=255)
                                                                                                                                                                      aaGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGga                                                                                                                                                                                                   <  7:562754/148‑1 (MQ=255)
                                                                                                                                                                         cACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAAcc                                                                                                                                                                                                                                      <  6:1051675/110‑1 (MQ=255)
                                                                                                                                                                          ttttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACCA‑TGGCACCGCTAACCCGCGATGc                                                                                                                                                                                                                              <  8:351795‑M1/28‑1 (MQ=255)
                                                                                                                                                                          aCTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGgaagaa                                                                                                                                                                                               <  7:782168/148‑1 (MQ=255)
                                                                                                                                                                           tttatgcaacaacaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATg                                                                                                                                                                                                                               >  7:676895‑M1/89‑116 (MQ=255)
                                                                                                                                                                            tGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAagaaga                                                                                                                                                                                             <  5:525895/148‑1 (MQ=255)
                                                                                                                                                                              gAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCa                                                                                                                                                                                                                                               <  8:244235/96‑1 (MQ=255)
                                                                                                                                                                               aaaCACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGga                                                                                                                                                                                                   <  7:299041/139‑1 (MQ=255)
                                                                                                                                                                               aaaCACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTc                                                                                                                                                                                          >  5:115877/1‑148 (MQ=255)
                                                                                                                                                                                 acacTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAAtg                                                                                                                                                                                                                 <  7:1024675/123‑1 (MQ=255)
                                                                                                                                                                                 acacTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGAt                                                                                                                                                                                            >  5:792569/1‑144 (MQ=255)
                                                                                                                                                                                   acTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAAt                                                                                                                                                                                  >  4:562413/1‑152 (MQ=255)
                                                                                                                                                                                      acaaattatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCAt                                                                                                                                                                                                                           >  6:785351‑M1/78‑109 (MQ=255)
                                                                                                                                                                                       ttCAGTGAGTATCTGTTAACCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATa                                                                                                                               <  2:108718/199‑1 (MQ=255)
                                                                                                                                                                                          aGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCt                                                                                                                    <  4:112528/207‑1 (MQ=255)
                                                                                                                                                                                           gTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgt                                                                                                                                                                              >  6:1065663/1‑148 (MQ=255)
                                                                                                                                                                                           gTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgt                                                                                                                                                                              >  8:847740/1‑148 (MQ=255)
                                                                                                                                                                                            tatccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACCA‑TGGCACCGCTAACCCGCGATg                                                                                                                                                                                                                               <  5:191821‑M1/27‑1 (MQ=255)
                                                                                                                                                                                            tGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          >  6:979456/1‑103 (MQ=255)
                                                                                                                                                                                            tGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgcgc                                                                              <  1:116570/243‑1 (MQ=255)
                                                                                                                                                                                               ccgttgcagttatgagtgacgtagatcacacttattcttcagaatcttcttccagctcgtcccacattGCACTA‑TGGCACCGCTAACCCGCGAt                                                                                                                                                                                                                                <  5:1031902‑M1/26‑1 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTgg                                                                                                                                                                                                    >  6:323318/1‑120 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGGTTCTga                                                                                                                                                                                       >  8:884002/1‑133 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  5:1061176/133‑1 (MQ=255)
                                                                                                                                                                                                 aTCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  7:884002/133‑1 (MQ=255)
                                                                                                                                                                                                   cTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      <  6:313345/148‑1 (MQ=255)
                                                                                                                                                                                                   cTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      >  7:552320/1‑148 (MQ=255)
                                                                                                                                                                                                     gTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCcgcg                                                                                                                                                                                                                                  <  6:188973/86‑1 (MQ=255)
                                                                                                                                                                                                     gTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGgaagaa                                                                                                                                                                                               <  6:463730/121‑1 (MQ=255)
                                                                                                                                                                                                      ttACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          >  5:561126/1‑93 (MQ=255)
                                                                                                                                                                                                      ttACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGga                                                                                                                                                                                                   <  5:315557/116‑1 (MQ=255)
                                                                                                                                                                                                      ttACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  5:87729/128‑1 (MQ=255)
                                                                                                                                                                                                       tACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTa                                                                                                                                                                                                                         <  7:1051564/93‑1 (MQ=255)
                                                                                                                                                                                                       tACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      >  6:719965/1‑144 (MQ=255)
                                                                                                                                                                                                       tACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                                                                      <  5:719965/144‑1 (MQ=255)
                                                                                                                                                                                                           cTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAAtg                                                                                                                                                                                                                 >  7:233833/1‑97 (MQ=255)
                                                                                                                                                                                                            tggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTc                                                                                                                                                                                          <  6:673163/119‑1 (MQ=255)
                                                                                                                                                                                                            tggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       >  6:1010941/1‑122 (MQ=255)
                                                                                                                                                                                                             ggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGATGAtt                                                                                                                                                                                           >  8:741284/1‑117 (MQ=255)
                                                                                                                                                                                                             ggtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATa                                                                                                                                                                                 <  8:1049619/127‑1 (MQ=255)
                                                                                                                                                                                                              gtgtggtgCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTg                                                                                                                                                       <  1:124331/152‑1 (MQ=255)
                                                                                                                                                                                                                        ttCCTTGGCAAAGTGGGGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGAAGAGCTGGAAGACGANTCTGAAGAATAAGTGTGATCTACNNCACCAATANCTGCAACGGATAATttg                                                                                                                                         >  3:30548/1‑156 (MQ=255)
                                                                                                                                                                                                                        ttCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACgg                                                                                                                                                 >  7:900899/1‑148 (MQ=255)
                                                                                                                                                                                                                         cacacttattcttcagaatcttcttccagctcgtcccacattGCACTAATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                                                                          <  6:513743‑M1/33‑1 (MQ=255)
                                                                                                                                                                                                                            ttGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  6:118844/106‑1 (MQ=255)
                                                                                                                                                                                                                              ggCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGAGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACCTACTTTCCTCTATTGCAATCAACACAGGCATTCTCTCCCGTTACTGGGGTTGGGTg            >  4:99721/1‑275 (MQ=255)
                                                                                                                                                                                                                               gCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAtt                                                                                                                                           <  8:291203/147‑1 (MQ=255)
                                                                                                                                                                                                                               gCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAtt                                                                                                                                           >  7:291203/1‑147 (MQ=255)
                                                                                                                                                                                                                                 aaaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgt                                                                                                                                        <  8:552320/148‑1 (MQ=255)
                                                                                                                                                                                                                                  aaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGa                                                                                                                                                >  5:33999/1‑139 (MQ=255)
                                                                                                                                                                                                                                  aaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgtt                                                                                                                                       <  6:850421/148‑1 (MQ=255)
                                                                                                                                                                                                                                  aaGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGTTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTgggg                   <  2:99722/264‑1 (MQ=255)
                                                                                                                                                                                                                                        aGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TTGCACCGCTAACCCCCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTccc  <  2:535913/275‑1 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGc                                                                                                                                                                                                       <  6:55677/73‑1 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGAtt                                                                                                                                                                                           <  6:99022/85‑1 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttg                                                                                                                                      >  5:261260/1‑138 (MQ=255)
                                                                                                                                                                                                                                             cATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttg                                                                                                                                      <  6:261260/138‑1 (MQ=255)
                                                                                                                                                                                                                                                gCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgtg                                                                                                                                                                             >  7:336475/1‑96 (MQ=255)
                                                                                                                                                                                                                                                gCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTc                                                                                                                                                                  >  5:479630/1‑107 (MQ=255)
                                                                                                                                                                                                                                                 cAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCTGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtgtg                                                                                                                        <  8:900899/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGa                                                                                                                                                                          >  7:858687/1‑95 (MQ=255)
                                                                                                                                                                                                                                                       cGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTg                                                                                                                                                                                        <  7:616595/78‑1 (MQ=255)
                                                                                                                                                                                                                                                       cGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTg                                                                                                                                                                                        >  8:616595/1‑78 (MQ=255)
                                                                                                                                                                                                                                                       cGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCg                                                                                                                  <  6:232968/148‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATa                                                                                                                                              <  5:883667/118‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgttg                                                                                                                                   <  6:359927/129‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTga                                                                                                                                                                                       <  6:877432/76‑1 (MQ=255)
                                                                                                                                                                                                                                                           tttCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCa                                                                                                                                 <  8:174873/129‑1 (MQ=255)
                                                                                                                                                                                                                                                            ttCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATaaa                                                                                                                             >  8:677509/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                    CACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgt                                                                                                                                        >  7:274918‑M1/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                      CCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttgttgt                                                                                                                                     <  5:986275‑M1/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCA‑TGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCC  >  NC_000913/1099862‑1100368

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑