Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,152,201 T→C 7.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,152,2010TC 7.2% 1024.9 / 44.5 332intergenic (+82/+35)murB/rneUDP‑N‑acetylenolpyruvoylglucosamine reductase/ribonuclease E
Reads supporting (aligned to +/- strand):  ref base T (129/178);  new base C (12/12);  total (141/191)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 340 total aligned reads displayed.

GCAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGAC  >  NC_002947/2152052‑2152348
                                                                                                                                                     |                                                                                                                                                   
  aaGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCCGTTTGc                                                                                                                                                   >  2:4201138/1‑150 (MQ=255)
     aTTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACt                                                                                                                                                <  2:1756843/150‑1 (MQ=255)
     aTTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACt                                                                                                                                                <  2:1949730/150‑1 (MQ=255)
          ggCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCgg                                                                                                                                           <  2:3259794/150‑1 (MQ=255)
           gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                          <  2:1586640/150‑1 (MQ=255)
           gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                          <  1:6356881/150‑1 (MQ=255)
           gCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTAAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCggg                                                                                                                                          <  2:1420045/150‑1 (MQ=255)
              gATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCt                                                                                                                                       >  1:1878021/1‑150 (MQ=255)
                 atCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtttt                                                                                                                                    <  1:3090536/150‑1 (MQ=255)
                   cTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCACAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCATTttc                                                                                                                                  <  2:6735808/150‑1 (MQ=255)
                   cTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttc                                                                                                                                  <  1:4111111/150‑1 (MQ=255)
                   cTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttc                                                                                                                                  <  1:5517589/150‑1 (MQ=255)
                    ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                 >  1:6707013/1‑150 (MQ=255)
                    ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                 >  2:2076119/1‑150 (MQ=255)
                    ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                 >  2:4520786/1‑150 (MQ=255)
                    ttCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttca                                                                                                                                 >  2:4592330/1‑150 (MQ=255)
                       gAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTttcattc                                                                                                                              <  2:289753/150‑1 (MQ=255)
                         gCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAg                                                                                                                            >  1:732770/1‑150 (MQ=255)
                         gCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGCAAACTGGCGGGGCtttttttggttt                                                                                                                            <  2:3592087/150‑8 (MQ=38)
                         gCGGTTTGGGGTAGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCATTATCATTCAg                                                                                                                            >  2:3794455/1‑150 (MQ=255)
                          cGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGc                                                                                                                           <  2:255438/150‑1 (MQ=255)
                               tGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTgg                                                                                                                      <  1:4832336/150‑1 (MQ=255)
                                ggggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGc                                                                                                                     <  2:1189717/150‑1 (MQ=255)
                                ggggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGc                                                                                                                     <  2:3179243/150‑1 (MQ=255)
                                 gggTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCt                                                                                                                    >  2:4703046/1‑150 (MQ=255)
                                   gTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTAt                                                                                                                  <  1:4964249/150‑1 (MQ=255)
                                           ggAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtg                                                                                                          >  2:1058806/1‑150 (MQ=255)
                                            gagaTGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGtga                                                                                                         <  1:3480071/150‑1 (MQ=255)
                                                 ggAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAgg                                                                                                    >  1:5143134/1‑150 (MQ=255)
                                                    gCCTAATCTTTACTGATTTGCCGGCTGTTAACCCCCGTAGGAGCGGTTTTACCCGCGATGCGGCCCTTCTGGAAAACAAAAAAAAGCCCGGCCTGTTTGCACTGGGGGGGATTTTTCATTCAGCTTTGGCTATCAGCGGTGTTTTGGttt                                                                                                 <  2:5843153/150‑1 (MQ=255)
                                                    gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                 <  1:1016666/150‑1 (MQ=255)
                                                    gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                 <  1:3753097/150‑1 (MQ=255)
                                                    gCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGttt                                                                                                 <  2:6248833/150‑1 (MQ=255)
                                                     ccTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtg                                                                                                >  1:6470260/1‑150 (MQ=255)
                                                       tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCtttttttgg                                                                                                                               <  2:6796357/117‑5 (MQ=255)
                                                       tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                              >  2:2481099/1‑150 (MQ=255)
                                                       tAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCACCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTtgtg                                                                                              >  1:2060362/1‑150 (MQ=255)
                                                         aTCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCt                                                                                            >  1:4949637/1‑150 (MQ=255)
                                                            tttACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCAt                                                                                         >  2:3659322/1‑150 (MQ=255)
                                                                cTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAg                                                                                                     <  1:6139474/134‑1 (MQ=255)
                                                                 tGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcggcg                                                                                    >  1:1052604/1‑150 (MQ=255)
                                                                 tGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATcggcg                                                                                    >  1:5770748/1‑150 (MQ=255)
                                                                     ttGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCtt                                                                                >  2:6117287/1‑150 (MQ=255)
                                                                       gCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTcc                                                                              <  1:21880/150‑1 (MQ=255)
                                                                        ccGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCa                                                                             <  2:932043/150‑1 (MQ=255)
                                                                         cGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAg                                                                            <  1:6624213/150‑1 (MQ=255)
                                                                             tgtgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCt                                                                        <  2:495839/150‑1 (MQ=255)
                                                                              gtgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTg                                                                       <  2:2740946/150‑1 (MQ=255)
                                                                               tgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTgg                                                                      >  1:2612955/1‑150 (MQ=255)
                                                                               tgAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGGTTGTGCTCATCGGCGGCTTCCAGGGCTgg                                                                      >  2:2025452/1‑150 (MQ=255)
                                                                                 aaCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGtt                                                                    >  2:3229291/1‑150 (MQ=255)
                                                                                 aaCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGtt                                                                    >  2:5192186/1‑150 (MQ=255)
                                                                                  aCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTc                                                                   <  2:1805998/150‑1 (MQ=255)
                                                                                   cccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTcc                                                                  >  1:6397314/1‑150 (MQ=255)
                                                                                    ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                 <  2:1258689/150‑1 (MQ=255)
                                                                                    ccccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCt                                                                 <  1:3229291/150‑1 (MQ=255)
                                                                                     cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                <  1:2338693/150‑1 (MQ=255)
                                                                                     cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTg                                                                <  2:726620/150‑1 (MQ=255)
                                                                                         ttggAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGaa                                                            <  2:2219485/148‑1 (MQ=255)
                                                                                         tAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGaa                                                            <  1:2797428/150‑1 (MQ=255)
                                                                                             aGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTg                                                        <  2:3339633/150‑1 (MQ=255)
                                                                                              gCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgc                                                       >  1:6579502/1‑150 (MQ=255)
                                                                                               cGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAg                                                           <  1:1786817/145‑1 (MQ=255)
                                                                                                  gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                   >  2:6833633/1‑150 (MQ=255)
                                                                                                     tACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGTTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCt                                                >  1:1171429/1‑150 (MQ=255)
                                                                                                       cccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                              <  1:6068835/150‑1 (MQ=255)
                                                                                                        ccGCGAAGAGGCCCGTCAGGCAAACCAAAAATAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGACGCTGCAGCAGCTTcc                                             >  2:3566156/1‑150 (MQ=255)
                                                                                                         cgcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTGTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCt                                            <  1:3786833/150‑1 (MQ=255)
                                                                                                          gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                           >  1:337183/1‑150 (MQ=255)
                                                                                                          gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                           >  2:550462/1‑150 (MQ=255)
                                                                                                            gAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtgtg                                         <  1:5810282/150‑1 (MQ=255)
                                                                                                             aaGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGc                                        <  2:2670950/150‑1 (MQ=255)
                                                                                                                  gCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCtt                                   >  2:468029/1‑150 (MQ=255)
                                                                                                                   cccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCttt                                  <  1:2405743/150‑1 (MQ=255)
                                                                                                                    ccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTg                                 >  2:2814239/1‑150 (MQ=255)
                                                                                                                    ccGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTg                                 <  1:6764203/150‑1 (MQ=255)
                                                                                                                     cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                <  1:3640445/150‑1 (MQ=255)
                                                                                                                     cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                <  2:5387926/150‑1 (MQ=255)
                                                                                                                     cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCAGTCAGCTTTGGCTATCAGCCGTGATGAGGGTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                >  2:148295/1‑150 (MQ=255)
                                                                                                                      gTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGTGGGTTGTGCTCATCGGCGGCTTCCAGGGCGGGTTCCTGGGAAGCTGCGGCAGCTTCCTGTGCCGCGTTGgc                               >  2:5114521/1‑150 (MQ=255)
                                                                                                                       tCAGGCAAACCAACAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgca                              <  1:1686802/150‑1 (MQ=255)
                                                                                                                          ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                           <  1:6863542/150‑1 (MQ=255)
                                                                                                                          ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                           <  2:6507011/150‑1 (MQ=255)
                                                                                                                          ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCATTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                           <  2:563511/150‑1 (MQ=255)
                                                                                                                           gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGcagcag                          >  2:336766/1‑150 (MQ=255)
                                                                                                                           gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTGTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGcagcag                          >  1:1527953/1‑150 (MQ=255)
                                                                                                                                ccAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCg                     >  2:6095894/1‑150 (MQ=255)
                                                                                                                                     aaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTc                <  1:4576068/150‑1 (MQ=255)
                                                                                                                                         gCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGc            <  1:5387038/150‑1 (MQ=255)
                                                                                                                                          ccccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCt           >  2:1922447/1‑150 (MQ=255)
                                                                                                                                           cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt          >  2:2444464/1‑150 (MQ=255)
                                                                                                                                           cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt          >  1:3949855/1‑150 (MQ=255)
                                                                                                                                           cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt          >  2:5403061/1‑150 (MQ=255)
                                                                                                                                           cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt          >  1:4193975/1‑150 (MQ=255)
                                                                                                                                            ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg         <  1:4703046/150‑1 (MQ=255)
                                                                                                                                             cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc        >  1:1895510/1‑150 (MQ=255)
                                                                                                                                              gCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcg       >  1:2718592/1‑150 (MQ=255)
                                                                                                                                                  ttttGCACTGGCGGTGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGATGTTTGTGCTCATCGGAGGCTTACAGTGCTGGTTCCTGGGAAGCTTCAGCATCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGa   <  2:3455361/149‑1 (MQ=255)
                                                                                                                                                  gTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGa   <  2:5723607/150‑1 (MQ=255)
                                                                                                                                                   tttGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGAc  <  1:52443/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCAAGATTCAGGCGGATATCTTCGAGCGGTTTGGGGTTGAGCTGGAGATGGAGCCTAATCTTTACTGATTTGCCGGCTGTGAACCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGAC  >  NC_002947/2152052‑2152348

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: