Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,647,673 T→A 8.6% intergenic (+22/+32) gor → / ← dinQ glutathione reductase (NADPH)/UV inducible membrane toxin DinQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,6730TA 8.6% 299.6 / 11.9 102intergenic (+22/+32)gor/dinQglutathione reductase (NADPH)/UV inducible membrane toxin DinQ
Reads supporting (aligned to +/- strand):  ref base T (46/47);  new base A (6/3);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCA  >  NC_000913/3647524‑3647820
                                                                                                                                                     |                                                                                                                                                   
gCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGagtagt                                                                                                                                                     <  2:1265351/150‑1 (MQ=255)
 cTTTGGTATGGACGAAATGTGGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGCGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAgtg                                                                                                                                                    <  2:1482446/150‑1 (MQ=255)
      gTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGctc                                                                                                                                               >  2:33190/1‑150 (MQ=255)
        aTGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCtt                                                                                                                                             >  1:52274/1‑150 (MQ=255)
         tGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTa                                                                                                                                            >  1:463635/1‑150 (MQ=255)
          ggACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAg                                                                                                                                           <  1:1027316/150‑1 (MQ=255)
          ggACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAAATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAg                                                                                                                                           >  2:63847/1‑150 (MQ=255)
           gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                          >  1:1066714/1‑150 (MQ=255)
           gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                          >  2:1159186/1‑150 (MQ=255)
           gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                          <  1:168674/150‑1 (MQ=255)
           gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGc                                                                                                                                          <  1:866573/150‑1 (MQ=255)
            aCGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAATACTTCGACAATACCGTCGCCATTCACCTAACGGCTTCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGcc                                                                                                                                         <  2:384898/150‑1 (MQ=255)
              gAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCt                                                                                                                                       >  2:1626411/1‑150 (MQ=255)
               aaaTGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCtt                                                                                                                                      <  2:399385/150‑1 (MQ=255)
                aaTGTTGCAGGGCTTCGCGGGGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTa                                                                                                                                     <  1:828092/150‑1 (MQ=255)
                  tGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGCCCtttaa                                                                                                                                   <  2:238249/150‑4 (MQ=255)
                    ttGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTa                                                                                                                                 <  1:79153/150‑1 (MQ=255)
                      gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACg                                                                                                                               >  2:502008/1‑150 (MQ=255)
                      gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACg                                                                                                                               >  1:1270655/1‑150 (MQ=255)
                       cAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGt                                                                                                                              <  2:1509419/150‑1 (MQ=255)
                         gggCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGttt                                                                                                                            >  2:854230/1‑150 (MQ=255)
                                gcgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTAt                                                                                                                     >  1:1268706/1‑150 (MQ=255)
                                gcgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTAt                                                                                                                     >  2:1310952/1‑150 (MQ=255)
                                 cgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATc                                                                                                                    <  1:1056318/150‑1 (MQ=255)
                                 cgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATc                                                                                                                    <  2:1268706/150‑1 (MQ=255)
                                   ggggCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAg                                                                                                                  <  1:1311706/148‑1 (MQ=255)
                                      gcgcTGAAGATGGGGGCAACCAAAAATGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCTTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTc                                                                                                               <  1:635763/150‑1 (MQ=255)
                                       ggcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                              <  1:1029274/149‑1 (MQ=255)
                                       cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                              <  2:450818/150‑1 (MQ=255)
                                           gAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAg                                                                                                          >  2:62974/1‑150 (MQ=255)
                                            aaGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGc                                                                                                         <  2:877779/150‑1 (MQ=255)
                                              gATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGACGCTg                                                                                                       >  2:792472/1‑150 (MQ=255)
                                                tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGaa                                                                                                     >  2:846774/1‑150 (MQ=255)
                                                tGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGaa                                                                                                     <  2:982929/150‑1 (MQ=255)
                                                  ggggCAACCAAAAAAGACTTGGAGAATACCGTCGCCATTAACCCAAGGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCGTTAATTAGGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  2:709816/150‑1 (MQ=255)
                                                  ggggCAACCAAAAAAGACTTCGACAATACCGTCGGCATTCACCCAACGGCGGCAGTAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  1:1508486/150‑1 (MQ=255)
                                                  ggggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  1:1585159/150‑1 (MQ=255)
                                                  ggggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   >  2:1216179/1‑150 (MQ=255)
                                                  ggggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  2:950117/150‑1 (MQ=255)
                                                  ggggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                   <  1:846774/150‑1 (MQ=255)
                                                   gggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCa                                                                                                  <  1:945391/150‑1 (MQ=255)
                                                   gggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCa                                                                                                  <  1:566108/150‑1 (MQ=255)
                                                    ggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAg                                                                                                 <  1:1626411/150‑1 (MQ=255)
                                                          cAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCg                                                                                           <  1:1310952/150‑1 (MQ=255)
                                                           aaaaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCgg                                                                                          <  2:413428/150‑1 (MQ=255)
                                                           aaaaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCgg                                                                                          <  2:52274/150‑1 (MQ=255)
                                                             aaaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGAt                                                                                        <  1:854230/150‑1 (MQ=255)
                                                                aGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAg                                                                                     <  1:604818/150‑1 (MQ=255)
                                                                     tagACAATAACGTCTCCATTCACCCAACTGAAGCAGAAGAGTTCGTGACAATTAGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCca                                                                                >  2:1430750/3‑150 (MQ=255)
                                                                     tCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCca                                                                                <  2:917212/150‑1 (MQ=255)
                                                                      cGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcag                                                                               <  1:1159186/150‑1 (MQ=255)
                                                                      cGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcag                                                                               >  1:172870/1‑150 (MQ=255)
                                                                       gACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagc                                                                              <  1:470342/150‑1 (MQ=255)
                                                                        aCAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagca                                                                             >  2:604693/1‑150 (MQ=255)
                                                                        aCAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagca                                                                             >  1:991045/1‑150 (MQ=255)
                                                                         cAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagcag                                                                            >  1:1365264/1‑150 (MQ=255)
                                                                             aCCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCa                                                                        <  2:384859/150‑1 (MQ=255)
                                                                              ccGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCaa                                                                       >  1:1210805/1‑150 (MQ=255)
                                                                              ccGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGCGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCaa                                                                       >  2:36498/1‑150 (MQ=255)
                                                                               cGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCCGCGCAAt                                                                      >  1:154284/1‑150 (MQ=255)
                                                                                gTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAg                                                                                                          >  1:1557358/1‑113 (MQ=255)
                                                                                gTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAg                                                                                                          <  2:1557358/113‑1 (MQ=255)
                                                                                gTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAAtt                                                                     <  1:526800/150‑1 (MQ=255)
                                                                                 tCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTa                                                                    <  1:183780/150‑1 (MQ=255)
                                                                                     cATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGc                                                                >  2:1238981/1‑150 (MQ=255)
                                                                                      aTTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGcc                                                               >  2:383886/1‑150 (MQ=255)
                                                                                          aCCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTa                                                           <  1:502008/150‑1 (MQ=255)
                                                                                           cccAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAg                                                          >  1:223045/1‑150 (MQ=255)
                                                                                           cccAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTACGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAg                                                          >  1:535061/1‑150 (MQ=255)
                                                                                                cggcggCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACga                                                     >  1:1558996/1‑150 (MQ=255)
                                                                                                cggcggCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACga                                                     >  1:578333/1‑150 (MQ=255)
                                                                                                    ggcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatgat                                                 <  2:172870/150‑1 (MQ=255)
                                                                                                    ggcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatgat                                                 <  1:344688/150‑1 (MQ=255)
                                                                                                       agaagaGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGc                                              >  1:104924/1‑150 (MQ=255)
                                                                                                        gaagaGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCt                                             >  2:899345/1‑150 (MQ=255)
                                                                                                                cGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaat                                     >  2:65772/1‑150 (MQ=255)
                                                                                                                 gTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaatc                                    <  2:147096/150‑1 (MQ=255)
                                                                                                                         gCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGttt                            <  2:435798/150‑1 (MQ=255)
                                                                                                                           gTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTtc                          >  2:1414123/1‑150 (MQ=255)
                                                                                                                           gTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTtc                          >  1:1066392/1‑150 (MQ=255)
                                                                                                                             tAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTtctc                        >  2:388571/1‑150 (MQ=255)
                                                                                                                                 taattaaGGGCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATg                    >  1:441156/6‑150 (MQ=39)
                                                                                                                                 tGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATg                    >  2:59460/1‑150 (MQ=255)
                                                                                                                                   ttAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCg                  <  2:282038/150‑1 (MQ=255)
                                                                                                                                    tAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGa                 >  2:1087558/1‑150 (MQ=255)
                                                                                                                                       cgggCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg              >  1:599318/2‑150 (MQ=255)
                                                                                                                                       agcgCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTACGCTATCAGTTCAGAACATGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg              >  2:502763/4‑150 (MQ=255)
                                                                                                                                       aGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg              <  2:690675/150‑1 (MQ=255)
                                                                                                                                       aGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg              <  2:1365264/150‑1 (MQ=255)
                                                                                                                                        gggCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGa             >  2:1059629/1‑150 (MQ=255)
                                                                                                                                        gggCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGa             >  1:159373/1‑150 (MQ=255)
                                                                                                                                         ggCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAg            >  1:1169813/1‑150 (MQ=255)
                                                                                                                                         ggCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAg            >  1:22265/1‑150 (MQ=255)
                                                                                                                                         ggCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAg            >  1:1393925/1‑150 (MQ=255)
                                                                                                                                          gCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGt           <  1:1216179/150‑1 (MQ=255)
                                                                                                                                           cTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATCGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTg          <  1:63847/150‑1 (MQ=255)
                                                                                                                                               gagTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCATCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGaa      <  2:331536/150‑1 (MQ=255)
                                                                                                                                                 gtagtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATg    <  2:535061/150‑1 (MQ=255)
                                                                                                                                                  tagtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGc   <  2:1066714/150‑1 (MQ=255)
                                                                                                                                                   agtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCa  <  1:65772/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCA  >  NC_000913/3647524‑3647820

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: