Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,040,126 G→T L277L (CTG→CTT appB → cytochrome bd‑II subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,040,1260GT87.9% 191.6 / 3.4 66L277L (CTG→CTTappBcytochrome bd‑II subunit 2
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (26/32);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGGAGCTA  >  NC_000913/1039879‑1040368
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ggATTCAGCTGCGGTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACg                                                                                                                                                                                                                                                 <  2:65741/250‑1 (MQ=255)
      atttgggTTCACAGCTGGCGACCAACCGCGCTGCACTGGTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATg                                                                                                                                                                                                                                           <  2:265650/245‑1 (MQ=255)
      aTCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGAGATGCTGTGCTTTTTGCTGGCGGGTTACTGGGTGTGGGGCGGTATTGATGGCTTTGTAATGCCCGCCCAGGATGCTAACGGTCCTTTCAATCCGTTAATGAAAATTGTGGCAGGGCTCCCTGGTGCCTGGATGAAAAATTTTGTCGCGTCGCCCGTTTTGTGGGTCTTCCCGCTGCCGGGGTTTTTCCGCCCATTTCTTTCGGTGATg                                                                                                                                                                                                                                           >  2:345857/1‑251 (MQ=255)
          gcggtCACAGCTGGCGACCAAACGCGCTGCTCTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTTCCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATgg                                                                                                                                                                                                                                          <  1:437789/244‑1 (MQ=255)
               cacaGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGCTGTGGGGCGGTATTGATGGCTTTGTAATGGTCGCCCAGGATGCTAACGGGCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTAACTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGGTTCTTCTGCCCATTGCTGACGGTGATGGCGGTTTAt                                                                                                                                                                                                                                  >  1:357358/1‑251 (MQ=255)
                      ggCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTc                                                                                                                                                                                                                           <  2:485426/251‑1 (MQ=255)
                                gctgctggatagTTGGTGTGGCTGTGCTTGTTGCTGGCGGGTTGCTGGCTGTGGGTGGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGGTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAAGAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCGGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTg                                                                                                                                                                                                                 <  2:167813/240‑1 (MQ=255)
                                cgcgCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTg                                                                                                                                                                                                                 <  2:226099/251‑1 (MQ=255)
                                   gcTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGGTGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTTGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGACGCCCGGGTTggg                                                                                                                                                                                                               >  1:287822/1‑250 (MQ=255)
                                   gcTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTgggg                                                                                                                                                                                                              <  1:69786/251‑1 (MQ=255)
                                   gcTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGCTGTGGGTCGGTATTGCGGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGGGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGGGTCGCCCGTTTTGTGGATCTTCCGGCTGCTGGGGTTCTTCTGCCCATTTCTGACGGGGGTGGCGATTTATCTGGGTCGCCCggggtgagg                                                                                                                                                                                                              >  1:167813/1‑245 (MQ=255)
                                   gcTGCACTGTTGGTGAGGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTTGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGACGCCCGGGTTg                                                                                                                                                                                                                 <  2:287822/248‑1 (MQ=255)
                                    cTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCCGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTggggg                                                                                                                                                                                                             >  1:295871/1‑250 (MQ=255)
                                     tGCACTGTTGGTGATGATGTGCTTTTTGCTGGCGGGTTAATGGCTGTGGGTCGGTATTGATGGCTTTGTTCTGCTCGCCCCGGGTGCTAACGGTCCTTCCAATCCGTTTATGTAACTGGGGGCAGTGCTACCCGGTGCCTGGGTGGATAATTTTGTCGGGTTGCCCGGTTTGTGGGTCTTTCCGCTGCTGGGGTTTTTCTGCCCATTGCTTACGGTTGTGGCGATTTATCGTGGTCGCCCGGGTTGGGGAt                                                                                                                                                                                                            >  2:219160/1‑251 (MQ=255)
                                      gcccggtgtgTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCGGTGGGTCGGTATTGATGGCTTTGTCCTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGAtt                                                                                                                                                                                                           <  2:509269/242‑1 (MQ=255)
                                      gCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGAtt                                                                                                                                                                                                           <  1:17844/251‑1 (MQ=255)
                                                       gtgCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTg                                                                                                                                                                                           <  1:116488/250‑1 (MQ=255)
                                                                    gCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGt                                                                                                                                                                             <  1:115951/251‑1 (MQ=255)
                                                                      gggTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGCTTTTTGATGGCATCATTGATGCAATTTGGCGTg                                                                                                                                                                            >  1:503423/1‑250 (MQ=255)
                                                                           aCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTc                                                                                                                                                                                                                           <  2:322724/198‑1 (MQ=255)
                                                                           aCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTc                                                                                                                                                                                                                           >  1:322724/1‑198 (MQ=255)
                                                                           aCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGGGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTTTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTc                                                                                                                                                                      <  2:154662/251‑1 (MQ=255)
                                                                           aCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTGCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGAtttt                                                                                                                                                                       <  1:345857/250‑1 (MQ=255)
                                                                             tGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGGGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCAc                                                                                                                                                                    >  1:53056/1‑251 (MQ=255)
                                                                              ggcggggggTCGGTATTGTTGGCTTTGTACTGCTAGCCCAGGATGCTCTCGGTCCTTCCAATCCGTTAATGTAACTGGTGGCAGTGCTACCTGGTGCCTGGATGCATATTTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCAc                                                                                                                                                                    <  1:474632/244‑1 (MQ=255)
                                                                              ggCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCAc                                                                                                                                                                    <  1:214072/250‑1 (MQ=255)
                                                                                               gATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGAATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGtt                                                                                                                                                  >  2:467638/1‑251 (MQ=255)
                                                                                                             gCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATg                                                                                                                                     <  1:283287/250‑1 (MQ=255)
                                                                                                                     aGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACCGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGc                                                                                                                            >  1:6247/1‑251 (MQ=255)
                                                                                                                      ggATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCg                                                                                                                           <  1:570251/251‑1 (MQ=255)
                                                                                                                           cTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCCCGCGGTTCCCCTTTGTCCTTCCGTCAAGCGTgggt                                                                                                                      >  2:43559/1‑248 (MQ=255)
                                                                                                                               cGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCgg                                                                                                                  >  1:247645/1‑250 (MQ=255)
                                                                                                                                      tCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCa                                                                                                            <  1:563908/250‑1 (MQ=255)
                                                                                                                                      tCCAATCCGTTAATGAAACTGGTGGCAGTGATACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTTGCGTGATTTTAACGGCAGGCATAACGCTGTTCCCCTTTGTCTTGCCGTTAGGCGTGAGTCCGATCTcc                                                                                                             >  2:472282/1‑249 (MQ=255)
                                                                                                                                         aaTCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCt                                                                                                        >  1:374789/1‑251 (MQ=255)
                                                                                                                                          aTCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTg                                                                                                       <  1:512309/251‑1 (MQ=255)
                                                                                                                                                   tGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTgg                                                                                              >  1:5868/1‑251 (MQ=255)
                                                                                                                                                     aaaCTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTg                                                                                                                                                                            >  1:307319/1‑171 (MQ=255)
                                                                                                                                                     aaaCTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTg                                                                                                                                                                            <  2:307319/171‑1 (MQ=255)
                                                                                                                                                             ggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGGGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACtt                                                                                    >  1:130397/1‑251 (MQ=255)
                                                                                                                                                                aGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCCTCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGGGTCCGATCTCCAGCCTTACGTTGTGGGACCGTCCTTcc                                                                                  >  2:236028/1‑250 (MQ=255)
                                                                                                                                                                aGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGAAGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGCCAAGCGTGAGTCCGATCTCCCACCTGTCGGTGTGGGGAAGTACTTccc                                                                                 >  2:537759/1‑250 (MQ=255)
                                                                                                                                                                  tGCTACCTGGTGCCTGGATGAATACTTTTGTCGCGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCCAGCGTGCGACCGATCTCCCGCCTGAAGTTTTGGGACCGTACTTcccgt                                                                               >  2:147564/1‑248 (MQ=255)
                                                                                                                                                                   gttcccTGGTGCCTGGATGAATAATTTTGCCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGt                                                                               <  2:294601/246‑1 (MQ=255)
                                                                                                                                                                      aCCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCCTGCCGTCAAGCGTGAGTCCGGTCTCCCGCCTGCCGTTGTGGGGCAGGACTTCCAGTCAGc                                                                           >  2:461204/1‑251 (MQ=255)
                                                                                                                                                                                 ggATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGc                                                                <  1:5130/251‑1 (MQ=255)
                                                                                                                                                                                   aTGAATAAATTTGTCGAGTCGCCCGTTTGGTGGAGCTTCCCGCTGCTTGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTGTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCAt                                                              <  2:123542/251‑1 (MQ=255)
                                                                                                                                                                                      aataatTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTAt                                                           <  1:236028/251‑1 (MQ=255)
                                                                                                                                                                                                  gAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGCCGTTGTGGGGCAGTACTTCCAGTCAGCTGACGCGGAGCATTATGTTGGTAGTCGt                                               >  2:434213/1‑251 (MQ=255)
                                                                                                                                                                                                   aGTCGCCCGTTGTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGACTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGtt                                              <  2:325259/251‑2 (MQ=255)
                                                                                                                                                                                                      cGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTTTTTTGCCCATTGCTGACGGTTATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTg                                           <  2:357358/251‑1 (MQ=255)
                                                                                                                                                                                                          cGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGatat                                       <  1:166940/251‑1 (MQ=255)
                                                                                                                                                                                                           gTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATAtt                                      <  1:40743/251‑1 (MQ=255)
                                                                                                                                                                                                             tttGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCTATTTAGCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCTTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATAttt                                     <  1:537759/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGGTATGGCGAGTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGTAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGccc                               <  1:531645/251‑1 (MQ=255)
                                                                                                                                                                                                                   gATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTACCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGccc                               >  1:540334/1‑250 (MQ=255)
                                                                                                                                                                                                                        tCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCAtt                            <  1:434213/248‑1 (MQ=255)
                                                                                                                                                                                                                        tCCCGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATtgtg                         >  1:471975/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGCTGCTGGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGGACTTCCAGTCAGCTGACGCTGAGCATTTTGTTGGTAATCGGGCTGATTTTTTTGCCCATTGTGttt                      >  2:156132/1‑250 (MQ=255)
                                                                                                                                                                                                                             ctgctgGGATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTAGCGTGGTCGCCCGGGTTGGGGCTTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACCATGTTCCCATTTGTCATGCCGTCAAGCGTGAGTCCGATATCCCACCTGCCGTTGTGGGACAGTACTTCCAGTCAGCTGAAGCTGAGCATTTTGTTTGTCATCGTGCTGTTCTTTTTGCCCCTTGTGTTGct                    >  2:106629/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gggATTCTTCTGCCCATTGCTTTCGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTa                 <  2:399343/249‑1 (MQ=255)
                                                                                                                                                                                                                                  gggATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGGCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGCAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTacc               >  1:300973/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ggATTCTTCTGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCAAGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTacac              >  1:499600/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tcttgccATTGCTGACGGTGAGGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAAGCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACActctct         <  2:295871/244‑1 (MQ=255)
                                                                                                                                                                                                                                          tctGCCCATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTg        <  1:258488/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ccATTGCTTACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGGAGcta  <  1:565105/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
TGATTCATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATCATTGATGCAATTTGGCGTGATTTTCACGGCAGGCATCACGCTGTTCCCCTTTGTCATGCCGTCAAGCGTGAGTCCGATCTCCAGCCTGACGTTGTGGGACAGTACTTCCAGTCAGCTGACGCTGAGCATTATGTTGGTAATCGTGCTGATATTTTTGCCCATTGTGTTGCTCTACACTCTCTGGAGCTA  >  NC_000913/1039879‑1040368

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: