Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,149,175 G→A 7.2% intergenic (+83/+34) 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,1750GA 7.2% 604.6 / 25.5 194intergenic (+83/+34)murB/rneUDP‑N‑acetylenolpyruvoylglucosamine reductase/ribonuclease E
Reads supporting (aligned to +/- strand):  ref base G (85/95);  new base A (5/9);  total (90/104)
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.99e-01

Only 100 of 198 total aligned reads displayed.

CAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCA  >  NC_002947/2149026‑2149324
                                                                                                                                                     |                                                                                                                                                     
 aaGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGc                                                                                                                                                      >  1:743691/1‑150 (MQ=255)
 aaGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGc                                                                                                                                                      >  1:3514203/1‑150 (MQ=255)
  aGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATGTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCa                                                                                                                                                     <  2:3555294/150‑1 (MQ=255)
     ttCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTg                                                                                                                                                  >  2:4394774/1‑150 (MQ=255)
     ttCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTg                                                                                                                                                  >  2:1691397/1‑150 (MQ=255)
        aGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCg                                                                                                                                               >  1:2675715/1‑150 (MQ=255)
         ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCgg                                                                                                                                              <  2:1754230/150‑1 (MQ=255)
         ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCgg                                                                                                                                              <  1:1669076/150‑1 (MQ=255)
         ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                              <  1:2134258/150‑1 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTTGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggt                                                                                                                                             <  2:908958/150‑2 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                             >  1:2845104/1‑150 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                             >  1:2845207/1‑150 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                             <  1:573109/150‑1 (MQ=255)
          gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCggg                                                                                                                                             <  1:955257/150‑1 (MQ=255)
              atatCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtt                                                                                                                                         >  1:4213686/1‑150 (MQ=255)
                atCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtttt                                                                                                                                       <  1:1317641/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:4432828/1‑150 (MQ=255)
                            tttGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCttt                                                                                                                           >  2:2853558/1‑150 (MQ=255)
                             ttGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTg                                                                                                                          <  2:2248760/150‑1 (MQ=255)
                              tGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTgg                                                                                                                         <  1:2962153/150‑1 (MQ=255)
                              tGGGGTGGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGTAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTgg                                                                                                                         <  1:1507853/150‑1 (MQ=255)
                                 ggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTa                                                                                                                      >  2:4256609/1‑150 (MQ=255)
                                   ttGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCtgat                                                                                                                    <  1:1480879/150‑4 (MQ=255)
                                   ttGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATc                                                                                                                    <  2:3983390/150‑1 (MQ=255)
                                    tGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCa                                                                                                                   <  2:2897747/150‑1 (MQ=255)
                                       gCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGcc                                                                                                                <  1:3681790/150‑1 (MQ=255)
                                        cTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGc                                                                                                                              >  2:2152844/1‑135 (MQ=255)
                                         tGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGt                                                                                                              >  2:3976509/1‑150 (MQ=255)
                                          ggAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAGAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtg                                                                                                             >  1:224560/1‑150 (MQ=255)
                                           gagtTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAATATGCCCCGCCAGTTTGCTCTGGCGGGGGTTTTTCATTCAGCTTTGGCTATCAGCCGtgt                                                                                                            >  2:4293445/1‑149 (MQ=39)
                                           gagaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtga                                                                                                            <  2:1314047/150‑1 (MQ=255)
                                           gagaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtga                                                                                                            <  1:4006113/150‑1 (MQ=255)
                                            agaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgat                                                                                                           <  2:3482262/150‑1 (MQ=255)
                                             gaTGGAGCCTAATCTTTACTGATTTGCCGGCCGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtgatg                                                                                                          >  2:3043281/1‑150 (MQ=255)
                                                 gAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGt                                                                                                      <  2:102035/150‑1 (MQ=255)
                                                   gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                    <  1:1388879/150‑1 (MQ=255)
                                                   gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                    >  1:799829/1‑150 (MQ=255)
                                                     cTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgt                                                                                                  <  2:1108934/150‑1 (MQ=255)
                                                      tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAATAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                                 >  2:1408940/1‑150 (MQ=255)
                                                      tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                                 <  1:1180234/150‑1 (MQ=255)
                                                              aCTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcgg                                                                                         <  1:1248043/150‑1 (MQ=255)
                                                               cTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcggc                                                                                        >  1:3446018/1‑150 (MQ=255)
                                                                tGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcggcg                                                                                       <  1:610522/150‑1 (MQ=255)
                                                                  aTTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCggcggc                                                                                     >  1:4080472/1‑150 (MQ=255)
                                                                     tGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCATTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTc                                                                                  <  1:4293445/150‑1 (MQ=255)
                                                                        cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                               <  1:1932718/150‑1 (MQ=255)
                                                                         ggCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAgg                                                                              >  2:4286129/1‑150 (MQ=255)
                                                                         ggCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAgg                                                                              >  2:2028479/1‑150 (MQ=255)
                                                                         ggCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAgg                                                                              >  2:2556495/1‑150 (MQ=255)
                                                                             gtgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTg                                                                          >  1:675900/1‑150 (MQ=255)
                                                                               gAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGt                                                                        >  1:2898506/1‑150 (MQ=255)
                                                                                   ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                    >  1:2160246/1‑150 (MQ=255)
                                                                                    cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                   <  2:281788/150‑1 (MQ=255)
                                                                                     ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAATCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTgg                                                                  <  1:461523/150‑1 (MQ=255)
                                                                                     ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTgg                                                                  <  1:1103290/150‑1 (MQ=255)
                                                                                            aGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTg                                                           <  1:3201657/150‑1 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTcc                                                                     >  2:599535/1‑135 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTcc                                                                     <  1:599535/135‑1 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                      >  2:730383/1‑150 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                      <  1:4457387/150‑1 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                      >  2:2996996/1‑150 (MQ=255)
                                                                                                 gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                      <  1:3034847/150‑1 (MQ=255)
                                                                                                      cccGCGAAGAGGCCAGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                 <  1:2078828/150‑1 (MQ=255)
                                                                                                       ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCAGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTcc                                                >  1:3391106/1‑150 (MQ=255)
                                                                                                         gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                              >  1:4349757/1‑150 (MQ=255)
                                                                                                         gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                              >  1:2794288/1‑150 (MQ=255)
                                                                                                         gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTAGGAAGCTGCTGCAGCTTccca                                              <  2:2005706/150‑3 (MQ=255)
                                                                                                          cgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGAGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgt                                             <  2:3730370/150‑1 (MQ=255)
                                                                                                                    cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                   <  1:870417/150‑1 (MQ=255)
                                                                                                                    cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                   <  2:743691/150‑1 (MQ=255)
                                                                                                                     gTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgc                                  >  2:175868/1‑150 (MQ=255)
                                                                                                                         ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                              >  2:3236641/1‑150 (MQ=255)
                                                                                                                           cAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCagcagc                            <  2:2679309/150‑1 (MQ=255)
                                                                                                                            aaaCCAAAAAATGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGcc                           >  2:3488912/1‑150 (MQ=255)
                                                                                                                            aaaCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGcc                           >  2:3743461/1‑150 (MQ=255)
                                                                                                                                cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                       >  1:1775680/1‑150 (MQ=255)
                                                                                                                                cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                       >  1:2524244/1‑150 (MQ=255)
                                                                                                                                cAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCgg                       <  2:3514340/150‑1 (MQ=255)
                                                                                                                                 aaaaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGc                      >  1:244659/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGccctggaag                    <  2:692682/150‑8 (MQ=255)
                                                                                                                                   aaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCt                    >  1:3147504/1‑150 (MQ=255)
                                                                                                                                    aaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTGTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTc                   <  2:4213686/150‑1 (MQ=255)
                                                                                                                                     aaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTcc                  <  1:3211372/150‑1 (MQ=255)
                                                                                                                                     aaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTcc                  <  1:1861256/150‑1 (MQ=255)
                                                                                                                                     aaaGCCCCGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTcc                  >  2:158570/1‑150 (MQ=255)
                                                                                                                                      aaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                 <  2:1561813/150‑1 (MQ=255)
                                                                                                                                      aaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                 <  1:2996080/150‑1 (MQ=255)
                                                                                                                                          cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt             >  2:450608/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt             >  1:167382/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt             >  1:394171/1‑150 (MQ=255)
                                                                                                                                          cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt             >  2:1929290/1‑150 (MQ=255)
                                                                                                                                           ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg            <  2:2108842/150‑1 (MQ=255)
                                                                                                                                           ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGAGCTCATCGGCGGCTTCCAGGGCAGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtag            >  2:2258044/1‑148 (MQ=255)
                                                                                                                                            cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc           <  1:84318/150‑1 (MQ=255)
                                                                                                                                              ccAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcgg         >  1:3858677/1‑150 (MQ=255)
                                                                                                                                                     gCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCa  >  2:3678274/1‑150 (MQ=255)
                                                                                                                                                     gCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCa  <  2:2747192/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCA  >  NC_002947/2149026‑2149324

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