Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,149,174 T→C 13.2% intergenic (+82/+35) murB → / ← rne UDP‑N‑acetylenolpyruvoylglucosamine reductase/ribonuclease E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,149,1740TC13.2% 339.1 / 34.3 130intergenic (+82/+35)murB/rneUDP‑N‑acetylenolpyruvoylglucosamine reductase/ribonuclease E
Reads supporting (aligned to +/- strand):  ref base T (47/65);  new base C (6/11);  total (53/77)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 133 total aligned reads displayed.

CAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACG  >  NC_002947/2149026‑2149322
                                                                                                                                                    |                                                                                                                                                    
cAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTg                                                                                                                                                     <  2:4623574/150‑1 (MQ=255)
cAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTg                                                                                                                                                     <  2:4623446/150‑1 (MQ=255)
   gATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCAc                                                                                                                                                  >  2:4781879/1‑150 (MQ=255)
     ttCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTg                                                                                                                                                >  1:1799631/1‑150 (MQ=255)
       cAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGc                                                                                                                                              <  2:2339708/150‑1 (MQ=255)
        aGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCg                                                                                                                                             >  2:2504505/1‑150 (MQ=255)
        aGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCg                                                                                                                                             >  2:1267267/1‑150 (MQ=255)
         ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                            <  1:3848794/150‑1 (MQ=255)
         ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                            <  2:4388712/150‑1 (MQ=255)
          gCGGATATCTTCGCGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATGTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGACGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                           <  1:31359/150‑1 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                           <  2:4642233/150‑1 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                           <  1:1027179/150‑1 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATGTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                           <  2:1177051/150‑1 (MQ=255)
            ggATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGc                                                                                                                                         >  1:459910/1‑150 (MQ=255)
              atatCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCGGGGCtt                                                                                                                                       <  1:1560204/150‑1 (MQ=255)
               tatCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCttt                                                                                                                                      >  1:4098150/1‑150 (MQ=255)
                atCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtttt                                                                                                                                     <  1:325569/150‑1 (MQ=255)
                  cTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttc                                                                                                                                   >  2:1845249/1‑150 (MQ=255)
                   ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                  >  1:2496324/1‑150 (MQ=255)
                   ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                  >  1:2496487/1‑150 (MQ=255)
                   ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                  >  1:507932/1‑150 (MQ=255)
                    acGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCGGGGCtttttttg                                                                                                                                 <  1:3003683/149‑4 (MQ=39)
                     cGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttcatt                                                                                                                                <  2:1229276/150‑1 (MQ=255)
                      gAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttcattc                                                                                                                               <  2:2292496/150‑1 (MQ=255)
                              tGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTgg                                                                                                                       <  2:3746691/150‑1 (MQ=255)
                                  gTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTAt                                                                                                                   >  2:3783734/1‑150 (MQ=255)
                                   ttGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATc                                                                                                                  <  2:2570118/150‑1 (MQ=255)
                                       gCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGcc                                                                                                              <  1:1368578/150‑1 (MQ=255)
                                           gagaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtga                                                                                                          >  1:127692/1‑150 (MQ=255)
                                              aTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgatga                                                                                                       >  2:213300/1‑150 (MQ=255)
                                                 gAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGt                                                                                                    <  1:745067/150‑1 (MQ=255)
                                                   gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                  <  2:127692/150‑1 (MQ=255)
                                                   gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTCTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                  <  1:1662057/150‑1 (MQ=255)
                                                    ccTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtg                                                                                                 <  1:854882/150‑1 (MQ=255)
                                                      tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTCCCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                               <  2:585376/150‑1 (MQ=255)
                                                      tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                           <  2:138860/106‑1 (MQ=38)
                                                      tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                           >  1:138860/1‑106 (MQ=38)
                                                           tttACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCAt                                                                                          <  1:1244981/150‑1 (MQ=255)
                                                            ttACTGATTTGCCGGCTGTGACCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATc                                                                                         <  1:3917556/150‑1 (MQ=255)
                                                             tACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcg                                                                                        >  2:696045/1‑150 (MQ=255)
                                                             tACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcg                                                                                        >  1:3018241/1‑150 (MQ=255)
                                                                tGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcggcg                                                                                     <  2:459910/150‑1 (MQ=255)
                                                                 gATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcggcgg                                                                                    >  2:301197/1‑150 (MQ=255)
                                                                  aTTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCggcggc                                                                                   >  1:2317375/1‑150 (MQ=255)
                                                                       ccGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCa                                                                              <  1:3421230/150‑1 (MQ=255)
                                                                       ccGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCa                                                                              <  1:21988/150‑1 (MQ=255)
                                                                        cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                             >  1:2958164/1‑150 (MQ=255)
                                                                        cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                             >  1:2958078/1‑150 (MQ=255)
                                                                            tgtgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCCTTCAGCGTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCt                                                                         >  2:1650010/1‑150 (MQ=255)
                                                                            tgtgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCt                                                                         >  2:3435285/1‑150 (MQ=255)
                                                                               gAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGt                                                                      >  1:4966859/1‑150 (MQ=255)
                                                                               gAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGt                                                                      >  2:364811/1‑150 (MQ=255)
                                                                               gAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGt                                                                      >  2:2536064/1‑150 (MQ=255)
                                                                                aaCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGtt                                                                     >  1:2896629/1‑150 (MQ=255)
                                                                                   ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                  <  2:4219131/150‑1 (MQ=255)
                                                                                   ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                  >  1:2022442/1‑150 (MQ=255)
                                                                                    cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                 <  2:680999/150‑1 (MQ=255)
                                                                                    cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                 >  2:5073266/1‑150 (MQ=255)
                                                                                    cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                 <  2:456336/150‑1 (MQ=255)
                                                                                    cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                 <  2:5004211/150‑1 (MQ=255)
                                                                                     ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                             >  2:1681036/1‑137 (MQ=255)
                                                                                     ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                             <  1:1681036/137‑1 (MQ=255)
                                                                                      cGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                           <  1:1043210/74‑1 (MQ=25)
                                                                                      cGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                           >  2:1043210/1‑74 (MQ=25)
                                                                                        tAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGaa                                                             <  2:3230017/150‑1 (MQ=255)
                                                                                           gAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGCTTCCCAGGaacca                                                          <  1:781360/150‑4 (MQ=255)
                                                                                             gCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgc                                                        <  2:507932/150‑1 (MQ=255)
                                                                                             gCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgc                                                        >  2:2547828/1‑150 (MQ=255)
                                                                                              cGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgca                                                       <  2:4571831/150‑1 (MQ=255)
                                                                                               gggTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcag                                                      <  1:2537336/150‑1 (MQ=255)
                                                                                                ggTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagc                                                     <  1:422641/150‑1 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                    <  2:1289943/150‑1 (MQ=255)
                                                                                                      cccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                               >  2:363241/1‑150 (MQ=255)
                                                                                                       ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTcc                                              >  2:2016132/1‑150 (MQ=255)
                                                                                                          cgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                           <  2:569994/150‑1 (MQ=255)
                                                                                                          cgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                           <  1:1151516/150‑1 (MQ=255)
                                                                                                           gAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgtg                                          >  2:3917636/1‑150 (MQ=255)
                                                                                                                 gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                    >  2:4878078/1‑150 (MQ=255)
                                                                                                                    cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                 <  1:4904569/150‑1 (MQ=255)
                                                                                                                    cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                 <  1:4878078/150‑1 (MQ=255)
                                                                                                                    cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                 <  1:804298/150‑1 (MQ=255)
                                                                                                                      tCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgca                               <  2:2617359/150‑1 (MQ=255)
                                                                                                                      tCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgca                               <  1:1924911/150‑1 (MQ=255)
                                                                                                                         ggaaaaaCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                            <  1:1004654/147‑1 (MQ=255)
                                                                                                                             aaCCAAAAAAAGCCCACCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTc                       <  2:1812694/150‑1 (MQ=255)
                                                                                                                                cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTGTCATTCAGCTTTGGCTATCCGCCGTGATGAGGGTTGTGCTCAGCGGCGGCTTCCAGGGCTGGTTCCTGGGACGCTGCAGCAGCGTCCTGTGCCGCTTTGGCAGCAGCCTCgg                     >  2:147544/1‑150 (MQ=255)
                                                                                                                                      aaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc               <  1:783036/150‑1 (MQ=255)
                                                                                                                                      aaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc               <  1:4521677/150‑1 (MQ=255)
                                                                                                                                        gCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTCCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGc             <  1:2169810/150‑1 (MQ=255)
                                                                                                                                         ccccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTgg                                                                <  1:1317877/98‑1 (MQ=255)
                                                                                                                                         ccccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCt            >  1:1510046/1‑150 (MQ=255)
                                                                                                                                         ccccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCt            >  1:1111488/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt           >  1:3144721/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt           >  2:4238438/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt           >  1:294010/1‑150 (MQ=255)
                                                                                                                                              ccAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcgg       <  2:2940589/150‑1 (MQ=255)
                                                                                                                                               cAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcggc      >  1:261938/1‑150 (MQ=255)
                                                                                                                                                  tttGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGACTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGAc   <  1:4068893/150‑1 (MQ=255)
                                                                                                                                                  tttGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATAGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGAc   >  2:1242251/1‑150 (MQ=255)
                                                                                                                                                   ttGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACg  <  1:2657669/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACG  >  NC_002947/2149026‑2149322

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