Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,149,176 C→A 7.2% intergenic (+84/+33) 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,1760CA 7.2% 652.1 / 27.1 193intergenic (+84/+33)murB/rneUDP‑N‑acetylenolpyruvoylglucosamine reductase/ribonuclease E
Reads supporting (aligned to +/- strand):  ref base C (83/96);  new base A (5/9);  total (88/105)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 196 total aligned reads displayed.

AAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCA  >  NC_002947/2149027‑2149324
                                                                                                                                                     |                                                                                                                                                    
aaGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGc                                                                                                                                                      >  1:3514340/1‑150 (MQ=255)
 aGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCa                                                                                                                                                     <  1:3887684/150‑1 (MQ=255)
 aGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATGTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCa                                                                                                                                                     <  2:3555294/150‑1 (MQ=255)
       aGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCg                                                                                                                                               >  1:2675715/1‑150 (MQ=255)
        ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCgg                                                                                                                                              <  1:1669076/150‑1 (MQ=255)
        ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                              <  2:1988846/150‑1 (MQ=255)
        ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                              <  1:3461930/150‑1 (MQ=255)
        ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                              <  1:2134258/150‑1 (MQ=255)
        ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGCGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                              <  2:68131/150‑1 (MQ=39)
         gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTTGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggt                                                                                                                                             <  2:908958/150‑2 (MQ=255)
         gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                             >  1:2845207/1‑150 (MQ=255)
         gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                             >  1:2845104/1‑150 (MQ=255)
         gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                             <  1:636843/150‑1 (MQ=255)
         gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                             <  1:585527/150‑1 (MQ=255)
            gATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCt                                                                                                                                          >  1:3102553/1‑150 (MQ=255)
               atCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtttt                                                                                                                                       <  1:1317641/150‑1 (MQ=255)
               atCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtttt                                                                                                                                       <  1:4078760/150‑1 (MQ=255)
                tCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCttttt                                                                                                                                      <  2:428120/150‑1 (MQ=255)
                 cTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttc                                                                                                                                     >  2:2940987/1‑150 (MQ=255)
                  ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                    >  1:215251/1‑150 (MQ=255)
                  ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                    <  1:317556/150‑1 (MQ=255)
                  ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                    >  1:4432951/1‑150 (MQ=255)
                   tCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttcat                                                                                                                                   >  1:607260/1‑150 (MQ=255)
                     gAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttcattc                                                                                                                                 >  2:4455642/1‑150 (MQ=255)
                             tGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTgg                                                                                                                         <  1:2962153/150‑1 (MQ=255)
                             tGGGGTGGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGTAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTgg                                                                                                                         <  1:1507853/150‑1 (MQ=255)
                              ggggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGc                                                                                                                        <  2:2446067/150‑1 (MQ=255)
                              ggggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGc                                                                                                                        <  2:2871021/150‑1 (MQ=255)
                               gggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCt                                                                                                                       >  2:1696893/1‑150 (MQ=255)
                                ggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTa                                                                                                                      >  2:4256609/1‑150 (MQ=255)
                                 gTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTAt                                                                                                                     <  1:508450/150‑1 (MQ=255)
                                  ttGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCtgat                                                                                                                    <  1:1480879/150‑4 (MQ=255)
                                  ttGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATc                                                                                                                    <  2:3983390/150‑1 (MQ=255)
                                   tGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCa                                                                                                                   <  2:2897766/150‑1 (MQ=255)
                                       cTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCg                                                                                                               >  2:2622773/1‑150 (MQ=255)
                                        tGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGt                                                                                                              >  2:3976509/1‑150 (MQ=255)
                                         ggAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAGAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtg                                                                                                             >  1:224560/1‑150 (MQ=255)
                                         ggAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtg                                                                                                             >  1:4351775/1‑150 (MQ=255)
                                          gagtTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAATATGCCCCGCCAGTTTGCTCTGGCGGGGGTTTTTCATTCAGCTTTGGCTATCAGCCGtgt                                                                                                            >  2:4293445/1‑149 (MQ=39)
                                          gagaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtga                                                                                                            <  1:99111/150‑1 (MQ=255)
                                           agaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgat                                                                                                           <  2:3482262/150‑1 (MQ=255)
                                            gaTGGAGCCTAATCTTTACTGATTTGCCGGCCGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgatg                                                                                                          >  2:3043281/1‑150 (MQ=255)
                                             aTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgatga                                                                                                         <  1:3438594/150‑1 (MQ=255)
                                             aTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgatga                                                                                                         <  2:738051/150‑1 (MQ=255)
                                                gAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGt                                                                                                      <  2:2778012/150‑1 (MQ=255)
                                                gAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGt                                                                                                      >  2:2075930/1‑150 (MQ=255)
                                                  gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                    <  2:2181600/150‑1 (MQ=255)
                                                  gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                    <  1:1388879/150‑1 (MQ=255)
                                                   ccTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtg                                                                                                   >  1:2996255/1‑150 (MQ=255)
                                                    cTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgt                                                                                                  <  2:2257317/150‑1 (MQ=255)
                                                    cTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgt                                                                                                  <  2:121919/150‑1 (MQ=255)
                                                     tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                                 <  2:146347/150‑1 (MQ=255)
                                                     tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                                 <  1:1180234/150‑1 (MQ=255)
                                                             aCTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcgg                                                                                         <  1:1248043/150‑1 (MQ=255)
                                                             aCTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcgg                                                                                         >  1:4469784/1‑150 (MQ=255)
                                                                      ccGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCa                                                                                >  2:2113639/1‑150 (MQ=255)
                                                                       cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                               <  2:1072674/150‑1 (MQ=255)
                                                                       cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                               <  1:1932718/150‑1 (MQ=255)
                                                                       cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                               >  1:765140/1‑150 (MQ=255)
                                                                        ggCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAgg                                                                              >  2:4286129/1‑150 (MQ=255)
                                                                        ggCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAgg                                                                              >  2:536508/1‑150 (MQ=255)
                                                                          cTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGc                                                                            <  2:1574235/150‑1 (MQ=255)
                                                                              gAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGt                                                                        >  1:2898506/1‑150 (MQ=255)
                                                                               aaCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGGTTGTGCTCATCGGCGGCTTCCAGGGCTGGtt                                                                       >  1:4343721/1‑150 (MQ=255)
                                                                                aCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTc                                                                      >  2:2840195/1‑150 (MQ=255)
                                                                                  ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                    >  1:2160246/1‑150 (MQ=255)
                                                                                  ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                    <  1:3254422/150‑1 (MQ=255)
                                                                                   cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                   <  1:3843278/150‑1 (MQ=255)
                                                                                    ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcgg                                                                                         <  2:1635953/127‑1 (MQ=255)
                                                                                    ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTgg                                                                  <  2:701607/150‑1 (MQ=255)
                                                                                    ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTgg                                                                  <  1:1103290/150‑1 (MQ=255)
                                                                                           aGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTg                                                           <  1:3201660/150‑1 (MQ=255)
                                                                                                gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                      >  2:2996996/1‑150 (MQ=255)
                                                                                                     cccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                 <  2:1926483/150‑1 (MQ=255)
                                                                                                     cccGCGAAGAGGCCAGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                 <  1:2078828/150‑1 (MQ=255)
                                                                                                        gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                              <  1:4087331/150‑1 (MQ=255)
                                                                                                        gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                              >  1:2794288/1‑150 (MQ=255)
                                                                                                         cgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGAGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                             <  2:3730370/150‑1 (MQ=255)
                                                                                                                gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                      >  1:4005362/1‑150 (MQ=255)
                                                                                                                 cccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCttt                                     <  2:2015636/150‑1 (MQ=255)
                                                                                                                  ccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTg                                    >  1:700476/1‑150 (MQ=255)
                                                                                                                    gTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgc                                  >  2:175849/1‑150 (MQ=255)
                                                                                                                        ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                              >  2:3236641/1‑150 (MQ=255)
                                                                                                                        ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                              >  2:1380010/1‑150 (MQ=255)
                                                                                                                         gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGcagcag                             >  1:807353/1‑150 (MQ=255)
                                                                                                                           aaaCCAAAAAATGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGcc                           >  2:3488912/1‑150 (MQ=255)
                                                                                                                           aaaCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGcc                           >  2:3743461/1‑150 (MQ=255)
                                                                                                                               cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                       <  2:3514340/150‑1 (MQ=255)
                                                                                                                               cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                       >  1:2524244/1‑150 (MQ=255)
                                                                                                                               cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                       >  1:1775680/1‑150 (MQ=255)
                                                                                                                                aaaaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGc                      <  1:925302/150‑1 (MQ=255)
                                                                                                                                aaaaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGc                      >  1:244659/1‑150 (MQ=255)
                                                                                                                                   aaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTGTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTc                   <  2:4213686/150‑1 (MQ=255)
                                                                                                                                    aaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTcc                  <  1:1861256/150‑1 (MQ=255)
                                                                                                                                     aaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                 <  1:2996080/150‑1 (MQ=255)
                                                                                                                                         cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt             >  1:167382/1‑150 (MQ=255)
                                                                                                                                          ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg            >  2:2741747/1‑150 (MQ=255)
                                                                                                                                          ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg            <  1:4256609/150‑1 (MQ=255)
                                                                                                                                           cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc           <  1:84318/150‑1 (MQ=255)
                                                                                                                                                    gCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCa  >  2:3678274/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCA  >  NC_002947/2149027‑2149324

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