Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,541 G→C 5.1% P28R (CCG→CGG)  yggN ← putative EcfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,5410GC 5.1% 260.8 / 4.2 79P28R (CCG→CGG) yggNputative EcfF
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base C (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTT  >  NC_000913/3101394‑3101689
                                                                                                                                                   |                                                                                                                                                    
tGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGGACGGTTTGCGGGCTGACAATCACATCGTCACgcggc                                                                                                                                                    >  1:871199/1‑150 (MQ=255)
 gATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACgcggcg                                                                                                                                                   >  1:1439713/1‑150 (MQ=255)
   tAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTg                                                                                                                                                 >  1:1336326/1‑150 (MQ=255)
    aaTCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGa                                                                                                                                                >  1:478088/1‑150 (MQ=255)
       ccTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGc                                                                                                                                             >  1:373655/1‑150 (MQ=255)
        cTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCt                                                                                                                                            <  2:1054407/150‑1 (MQ=255)
           cGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCTGCGTGACGCTGCa                                                                                                                                         <  2:333975/150‑1 (MQ=255)
           cGCCTGCTCGCGCCGCGCGGCACCCAGCCAACATTGCCTCCCGTCCTCCATCACGTTGCCGTCTGGCGCGATCACCAGATTGCCGCTTTCGCCCTTCACCTGCACGGCTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACGATgtg                                                                                                                                         <  2:204202/150‑3 (MQ=255)
              cTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCCCCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCAGTg                                                                                                                                      >  1:1166936/1‑150 (MQ=255)
                gCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGt                                                                                                                                    <  2:479169/150‑1 (MQ=255)
                gCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGt                                                                                                                                    <  2:1277276/150‑1 (MQ=255)
                gCTCGCGCGGGGCCGCATTCAGCGAATATTGCTCACCGCCATACACCACGTTGCCGGCTCGCGCGACCACCAGATTCCCGCTTTCGCCCTTCACCTGCACGGTTTGCGGGCGGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGt                                                                                                                                    <  2:420220/150‑1 (MQ=255)
                     cgcTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACg                                                                                                                                              <  2:1296359/135‑1 (MQ=255)
                     cgcTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCCGCGTGACg                                                                                                                                              >  1:1296359/1‑135 (MQ=255)
                               gaTTCAGGGATTATTGGTTACCGTTATACATCCCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTCGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                                     <  1:215995/149‑1 (MQ=255)
                               cATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCTGCGTGACGCTGCACTGGTAGTCGGCGTGAGCgg                                                                                                                     <  1:576593/150‑1 (MQ=255)
                                        aaTATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCg                                                                                                            >  2:761014/1‑150 (MQ=255)
                                             tGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACt                                                                                                       <  2:1336326/150‑1 (MQ=255)
                                              gCTTACCGTTATACATCACGTTGCCGTATTGCGTGATCACAAGATTGCCGTTTTCGCCCTTAACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTg                                                                                                      >  2:311042/1‑150 (MQ=255)
                                                 tACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCGCGGTTTGGGGGCTGACAGTCACATCGTCACGCGGCGTGACGCTGCGGTGGTTGTCGGCGTGAGCGGTCATTGCCGTCACTGaat                                                                                                   >  1:1272101/1‑149 (MQ=255)
                                                 tACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                   >  1:670629/1‑150 (MQ=255)
                                                 tACCGTTATACATCACGTTGCCGTCTGGCGTGATCAACAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGaaa                                                                                                   >  2:1135801/1‑150 (MQ=255)
                                                  aCCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAg                                                                                                  >  1:783238/1‑150 (MQ=255)
                                                     gTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgc                                                                                               >  1:447331/1‑150 (MQ=255)
                                                      ttATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGtgct                                                                                              <  2:1333872/150‑1 (MQ=255)
                                                           cATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTgccgc                                                                                         >  2:529674/1‑150 (MQ=255)
                                                              cACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGGcag                                                                                      >  1:1475686/1‑150 (MQ=255)
                                                               aCGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagc                                                                                     <  2:373655/150‑1 (MQ=255)
                                                                   tGCCGTATTCCGTGATCACCAGAATGCAGTTTTCGCCCTTAACCTCCACGGTTTGAGTGCTGACAATCACATCGTCACGCGGCGTGACGCTCCACTGGAACTCGGCGTCAGAGGACATTGCCGTCACTGAAAGTGCTGCCGCCAgcagca                                                                                 >  2:203021/1‑150 (MQ=255)
                                                                    gCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAt                                                                                >  1:385183/1‑150 (MQ=255)
                                                                     ccGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAtt                                                                               >  2:350101/1‑150 (MQ=255)
                                                                      ggTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAttt                                                                              <  2:592152/149‑1 (MQ=255)
                                                                        tCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCTGGCGGCAGCActttc                                                                            <  1:1240256/150‑2 (MQ=255)
                                                                        tCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAATGCTGCTGGCGGCAGCActttc                                                                            <  2:615869/150‑2 (MQ=255)
                                                                           ggCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTgcgc                                                                         <  1:371300/150‑1 (MQ=255)
                                                                           ggCGTCACCACCCCATTGCCGTTTCCGCCCTTCACCTCCACGTTTTGCGGGCTGACAATCACATCGCCACGCGGCGTCACGCTGCACTGGTAGTCGGCGTGCGCGGTCATTGCCGTCACTGAAACTGCTGCCGGCGGCAGCActttcagt                                                                         <  2:142715/150‑5 (MQ=255)
                                                                               tGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagc                                                                                     <  1:896420/134‑1 (MQ=255)
                                                                               tGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCcagc                                                                                     >  2:896420/1‑134 (MQ=255)
                                                                               tGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatca                                                                     >  2:431902/1‑150 (MQ=255)
                                                                                gATCACCTGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCAAGCGGCGTGACGCTGCACTGGTTTTCGGCGGTAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCTcatcct                                                                    <  2:1047337/150‑3 (MQ=255)
                                                                                gATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGcatcat                                                                    <  2:1439713/150‑1 (MQ=255)
                                                                                 aTCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATa                                                                   >  2:985948/1‑150 (MQ=255)
                                                                                  tCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAg                                                                  >  1:811441/1‑150 (MQ=255)
                                                                                  tCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAg                                                                  <  2:478088/150‑1 (MQ=255)
                                                                                  tCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAg                                                                  <  1:1400145/150‑1 (MQ=255)
                                                                                  tCACCAGATTGCCGTGTTCGCCCTGCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAg                                                                  <  1:555677/150‑1 (MQ=255)
                                                                                   cACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                 >  2:1052270/1‑150 (MQ=255)
                                                                                   cACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGt                                                                 >  1:1494760/1‑150 (MQ=255)
                                                                                    aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCATGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTc                                                                >  2:1021414/1‑150 (MQ=255)
                                                                                    aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTc                                                                <  1:1135801/150‑1 (MQ=255)
                                                                                    aCCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTGTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTc                                                                >  2:1219112/1‑150 (MQ=255)
                                                                                       aGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTc                                                             <  1:148207/150‑1 (MQ=255)
                                                                                         aTTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTccc                                                           <  1:1021113/150‑1 (MQ=255)
                                                                                         aTTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTccc                                                           <  2:670629/150‑1 (MQ=255)
                                                                                         aTTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTccc                                                           <  1:431902/150‑1 (MQ=255)
                                                                                         aTTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACCCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCCCATCATAGTCTTccc                                                           <  2:1272101/150‑1 (MQ=255)
                                                                                            gCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCa                                                        >  1:1054816/1‑150 (MQ=255)
                                                                                               gTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                     >  2:958269/1‑150 (MQ=255)
                                                                                                  ttCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAg                                                      <  2:1509351/146‑1 (MQ=255)
                                                                                                  ttCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAg                                                      >  1:1509351/1‑146 (MQ=255)
                                                                                                  ttCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaa                                                  >  1:375593/1‑150 (MQ=255)
                                                                                                   tCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                 >  1:270459/1‑150 (MQ=255)
                                                                                                    cGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                <  1:1272874/150‑1 (MQ=255)
                                                                                                        cTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAg                                            <  2:1535625/150‑1 (MQ=255)
                                                                                                           cACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                         >  2:1252654/1‑150 (MQ=255)
                                                                                                            aCCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                        >  2:54635/1‑150 (MQ=255)
                                                                                                             ccTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                       <  2:447331/150‑1 (MQ=255)
                                                                                                               tGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                     >  1:1257590/1‑150 (MQ=255)
                                                                                                                 cACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                   <  1:449217/150‑1 (MQ=255)
                                                                                                                        tGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                            >  1:617738/1‑150 (MQ=255)
                                                                                                                        tGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                            >  2:674749/1‑150 (MQ=255)
                                                                                                                              cTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                      >  1:1201787/1‑150 (MQ=255)
                                                                                                                                gACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                    >  2:7827/1‑150 (MQ=255)
                                                                                                                                 aCAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                   <  1:1126908/150‑1 (MQ=255)
                                                                                                                                 aCAATCACATCGTCACGCGGCGGGACGCTGCACTGGTAGTCGGCGTGAGCGGGCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                   <  2:218269/150‑1 (MQ=255)
                                                                                                                                    aTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                >  2:1415191/1‑150 (MQ=255)
                                                                                                                                       gcagcGTCACGCCGCGTGACGATGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACa             >  2:1529661/5‑150 (MQ=255)
                                                                                                                                        caTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa            >  1:787324/1‑150 (MQ=255)
                                                                                                                                        caTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa            <  1:1052270/150‑1 (MQ=255)
                                                                                                                                         aTCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa           >  2:265290/1‑150 (MQ=255)
                                                                                                                                          tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg          <  1:1434439/150‑1 (MQ=255)
                                                                                                                                                 gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt   <  2:385183/150‑1 (MQ=255)
                                                                                                                                                  cggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt  >  1:534807/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGATAATCCTTCGCCTGCTCGCGCTGGGCGGCATTCAGGGAATATTGCTTACCGTTATACATCACGTTGCCGTCTGGCGTGATCACCAGATTGCCGTTTTCGCCCTTCACCTGCACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTT  >  NC_000913/3101394‑3101689

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