Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,403,597 T→C 21.1% intergenic (‑57/+24) PP_1227 ← / ← PP_1228 membrane protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,403,5970TC21.1% 449.5 / 122.5 243intergenic (‑57/+24)PP_1227/PP_1228membrane protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (81/103);  new base C (26/24);  total (108/135)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-05

Only 100 of 246 total aligned reads displayed.

GGTGGCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAC  >  NC_002947/1403448‑1403746
                                                                                                                                                     |                                                                                                                                                     
ggtggCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGt                                                                                                                                                       <  2:3009080/150‑1 (MQ=255)
    gCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGccc                                                                                                                                                   <  2:7060277/150‑1 (MQ=255)
      gATGAAGAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgc                                                                                                                                                 <  2:7249926/150‑1 (MQ=255)
      gATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgg                                                                                                                                                 >  1:6215158/1‑149 (MQ=255)
      gATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgc                                                                                                                                                 <  1:660028/150‑1 (MQ=255)
         gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                              <  2:7809244/150‑1 (MQ=255)
         gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              <  1:7009971/150‑1 (MQ=255)
           aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                            >  2:6164976/1‑150 (MQ=255)
           aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                            <  2:787418/150‑1 (MQ=255)
           aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                            <  2:4122486/150‑1 (MQ=255)
           aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                            <  2:3638980/150‑1 (MQ=255)
            cAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAagag                                                                                                                                           <  1:2916847/150‑1 (MQ=255)
                 tGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCt                                                                                                                                      <  2:6725467/150‑1 (MQ=255)
                   gcggcgATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGa                                                                                                                                    <  2:7105750/150‑1 (MQ=255)
                     ggcgATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACa                                                                                                                                  >  2:5277272/1‑150 (MQ=255)
                        gATCGAGGTGCGTAAAAGCCCTTGTTCTCTGGTGATGTCGAAGAAATTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTAGCGGGCGTGCCCGCGAAGGGGCCCTGACAGGc                                                                                                                               >  1:2723501/1‑150 (MQ=255)
                        gATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGc                                                                                                                               >  2:6832527/1‑150 (MQ=255)
                            gAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGa                                                                                                                           >  1:8342613/1‑150 (MQ=255)
                            gAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCCTGACAGGCTAGa                                                                                                                           >  2:182209/1‑150 (MQ=255)
                             cggTGCGTACAAGCCCTTGTTCGCGGGTGATGTCGAAGAAAGTGGATCGGGTTGGCTGCGGCATAAGGTGTGCTCAAGAATACCGCTGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGAc                                                                                                                          <  2:9161238/149‑1 (MQ=38)
                                tGCGTAAAAGCCCTGGTTCGCTGGTGATGTCGAAGAACGTGGCCCGGGTTGGATTAGGCATATGTTGTGCACAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGAcgcg                                                                                                                       <  1:561302/150‑1 (MQ=38)
                                 gCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGTTGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGa                                                                                                                      <  1:3602279/150‑1 (MQ=255)
                                 gCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGa                                                                                                                      >  2:3119674/1‑150 (MQ=255)
                                  cGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGaa                                                                                                                     >  1:9158943/1‑150 (MQ=255)
                                    tAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAAc                                                                                                                   >  2:5119343/1‑150 (MQ=255)
                                      aaaGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTg                                                                                                                 >  1:2112805/1‑150 (MQ=255)
                                                ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                       <  2:5164083/150‑1 (MQ=255)
                                                ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                       <  2:8569259/150‑1 (MQ=255)
                                                ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                       <  1:5119343/150‑1 (MQ=255)
                                                ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                       >  2:8090535/1‑150 (MQ=255)
                                                    cTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCct                                                                                                   <  2:4288832/150‑1 (MQ=255)
                                                     tGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctg                                                                                                  >  2:4570550/1‑150 (MQ=255)
                                                       gTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              >  1:6991832/1‑104 (MQ=255)
                                                        tGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATACGTTGTGCTCAAGAATAACGATTGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctg                                                                                               <  2:2398030/150‑1 (MQ=255)
                                                        tGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctg                                                                                               >  2:4412478/1‑150 (MQ=255)
                                                         tttGTCTAAGAATGTGGATCGGGTTGGATTAGGCATATGTTGTGCTAAATAATACCGTTGGTGCGGTGCCTGTCAGGGCCTCTGCGGGGGCGGGCCCGCCAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                              <  1:790806/148‑1 (MQ=255)
                                                         gATGTCTACGAATGTGGTGGGGGTTGTATTAGGCATATGTTGGGCTCACGAATACCGAGGGGGCGGTGCCTGCCCGGGCCTCTTCGGGGGGGTGCCCGCGACGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                              <  1:1752322/150‑1 (MQ=255)
                                                               gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTg                                                                                        >  2:2551828/1‑150 (MQ=255)
                                                                 agaaAGTGTTTCGGGTTGGATTCGGCATAAGTTGTGCTCAAGAATAACGTTGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGc                                                                                      <  1:4963623/150‑1 (MQ=255)
                                                                       tGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCg                                                                                <  1:1449210/150‑1 (MQ=255)
                                                                         gATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACCCGAAAAAGCTCCAACAGCGCCAGCTGCTGGTTGGCCCGCCGGt                                                                              >  1:6527471/1‑150 (MQ=38)
                                                                           tCGGGTTGGTTTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTc                                                                            <  1:7272050/150‑1 (MQ=255)
                                                                               gTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCGGCTGCTGGTTGGCCAGCCGGTTCAggg                                                                        >  2:4117271/1‑150 (MQ=255)
                                                                                 tGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCt                                                                      <  2:2112805/150‑1 (MQ=255)
                                                                                  ggATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCtt                                                                     <  1:8585712/150‑1 (MQ=255)
                                                                                  ggATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTAGCGGGCGTGCCCCCGAAGAGGCCCAGACAGGCTAGAAGCGAAACAGCACCAAAACCGCCAGCAGCTGGTTTGCCACCCCGTTCAGGGCtt                                                                     >  1:6182176/1‑150 (MQ=255)
                                                                                     ttAGGCATACGTTGTGGTCACGTATAACTATGGTGCGGTGCCTGTCTGGGCCGGTTGGGGGGGGTGCCCGCGAAGAGGGCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCTGCCGGTTCAGGGCTTGGc                                                                  <  1:5896198/150‑1 (MQ=255)
                                                                                     ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGCAGAGGCCCAGACAGGCTAGACGCGAAAATGCTCCAACAGCGCCGGCTGCTGGGTGGCCACCCCGTTCAGGGCGTGGc                                                                  >  1:3930112/1‑150 (MQ=38)
                                                                                     ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                  >  2:8131034/1‑150 (MQ=255)
                                                                                     ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                  <  2:3442053/150‑1 (MQ=255)
                                                                                     ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGGTCAGGGCTTGGc                                                                  >  2:7266110/1‑150 (MQ=255)
                                                                                         gCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGAt                                                              >  1:3563034/1‑150 (MQ=255)
                                                                                         gCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTAGCGGGCGTGGCCGGGCAGATGCCCAGACAGGCTAGACGCGCAAAAGCTCCAACAGCGCCGGCTGCGGGTTGGCCCGCCGGGTCCGGGCGTGGCTGAt                                                              >  1:3277669/1‑150 (MQ=255)
                                                                                          cATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCGTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATg                                                             >  2:1641736/1‑150 (MQ=255)
                                                                                                ttGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGcc                                                       >  2:6344211/1‑150 (MQ=255)
                                                                                                  gtccTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGa                                                     <  2:2288520/147‑1 (MQ=255)
                                                                                                     cTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTc                                                  <  2:1220999/150‑1 (MQ=255)
                                                                                                        aaGAATAACGATGGTGCGGTGCCTGTCAGTGCCTCTTCGCGGGGGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAATCTGCTCGATCATCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGTCTGATGCGTGCCGATTCATc                                               <  2:8732759/150‑1 (MQ=39)
                                                                                                         aGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                              <  1:2186488/150‑1 (MQ=255)
                                                                                                         aGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTGCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                              <  1:8171548/150‑1 (MQ=255)
                                                                                                          gAATAACGATGGTGCGGTGCCTGTCAGGGCCTATTGGCGGGCGTGCCCGCGATGAGGCCCTGACAGGGTAGTCGCGAAAGTTGTCCATGAGCGGCGGCTGCTGGTTGGCCAGCCGGTTGAGGGCTTGGCTGATGCGTGCCGATTCATCgg                                             >  2:3257573/1‑150 (MQ=255)
                                                                                                              aaCGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAGCTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTg                                         >  2:1742590/1‑150 (MQ=255)
                                                                                                                 gATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGccc                                      <  1:487012/150‑1 (MQ=255)
                                                                                                                  aTGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCcca                                     <  2:5364697/150‑1 (MQ=255)
                                                                                                                  aTGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCcca                                     <  2:6553587/150‑1 (MQ=255)
                                                                                                                   tGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccag                                    <  1:9058412/150‑1 (MQ=255)
                                                                                                                     gtgcggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagcc                                  >  2:8215988/1‑150 (MQ=255)
                                                                                                                      tgcggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGTGGCCCTGACAGGGTAGACGCGAAACTGGTCCATGAGCGCCTGGTGCTGGTTGGCCAGCCGGTTGAGGGCTTGGCTGATGCGTGGCGTTTCATCGGCCTGCccagcca                                 >  2:83464/1‑150 (MQ=38)
                                                                                                                      tgcggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagcca                                 <  1:4365036/150‑1 (MQ=255)
                                                                                                                      ggcggggcCTGTCTGGGCCTCGTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagcca                                 <  1:162944/149‑1 (MQ=255)
                                                                                                                       gcggtgcCTGTCAGGGCCTCTTGGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGCAACAGCTCCATCAGCGCCGGCTGCGGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCGGCccaccccc                                >  1:5700173/1‑148 (MQ=38)
                                                                                                                       gcggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagccag                                <  1:1472441/150‑1 (MQ=255)
                                                                                                                         ggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGGCGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTg                              <  1:6398845/150‑1 (MQ=255)
                                                                                                                          gtgcCTGTCAGGGCCTCTTCTCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGtta                             <  1:745734/150‑3 (MQ=255)
                                                                                                                              cTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCt                         <  2:463702/150‑1 (MQ=255)
                                                                                                                                 tCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTc                      <  2:5746136/150‑1 (MQ=255)
                                                                                                                                 tCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTc                      >  2:1433602/1‑150 (MQ=255)
                                                                                                                                 tCAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCAGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTc                      >  2:2155395/1‑150 (MQ=255)
                                                                                                                                   tgggCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCAc                    <  1:30368/149‑1 (MQ=255)
                                                                                                                                   aGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCAc                    <  1:6832527/150‑1 (MQ=255)
                                                                                                                                        ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtctc               <  1:5277272/150‑1 (MQ=255)
                                                                                                                                         tctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:566923/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  1:5539144/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  1:8755180/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:7463403/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  1:2226129/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCACGCCCGCGAAGAGGCCCTCACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:693250/1‑150 (MQ=255)
                                                                                                                                          ctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGa             >  1:6596998/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGaa            >  2:905649/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGaa            >  2:292761/1‑150 (MQ=255)
                                                                                                                                             cgcgGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATg          >  1:1261233/1‑150 (MQ=255)
                                                                                                                                              gcgGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGc         >  1:1866069/1‑150 (MQ=255)
                                                                                                                                              gcgGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGc         >  1:3411315/1‑150 (MQ=255)
                                                                                                                                              gcgGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGc         >  1:483336/1‑150 (MQ=255)
                                                                                                                                              gcgGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGc         >  1:5340789/1‑150 (MQ=255)
                                                                                                                                                  gCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGc     <  1:8215988/150‑1 (MQ=255)
                                                                                                                                                  gCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGc     <  2:5346669/150‑1 (MQ=255)
                                                                                                                                                    gTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCa   <  1:2551828/150‑1 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAc  <  1:1036585/150‑1 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAc  <  1:4257472/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGGCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAC  >  NC_002947/1403448‑1403746

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: