Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,408,546 T→C 21.3% 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,408,5460TC21.3% 494.6 / 133.6 264intergenic (‑57/+24)PP_1227/PP_1228membrane protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (88/105);  new base C (27/27);  total (119/145)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-02

Only 100 of 268 total aligned reads displayed.

GTGGCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAC  >  NC_002947/1408398‑1408695
                                                                                                                                                    |                                                                                                                                                     
    cGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcg                                                                                                                                                  >  1:6842231/1‑150 (MQ=255)
     gATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgc                                                                                                                                                 <  2:6966726/150‑1 (MQ=255)
     gATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgc                                                                                                                                                 >  2:2132519/1‑150 (MQ=255)
       tGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGa                                                                                                                                               <  1:3570173/150‑1 (MQ=255)
       tGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGa                                                                                                                                               <  1:1087899/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAACAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              <  1:788818/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGCTCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              <  1:1013357/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                              <  2:1667970/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                              <  2:777551/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                              <  1:6010082/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                              <  1:4758452/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              <  2:6520187/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              <  1:2731628/150‑1 (MQ=255)
        gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGAGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                              <  1:1196421/150‑1 (MQ=255)
          ccAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCCTAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                            <  2:7002170/149‑1 (MQ=255)
          aCAGGGTGGCGGCGATCGTGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCTAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                            <  2:365572/150‑1 (MQ=255)
          aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                            <  1:2951290/150‑1 (MQ=255)
          aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                            <  1:4385412/150‑1 (MQ=255)
          aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                            <  2:1077946/150‑1 (MQ=255)
          aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                            <  1:734151/150‑1 (MQ=255)
          aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                            <  1:6059681/150‑1 (MQ=255)
           cAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTGCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCGCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAagag                                                                                                                                           <  1:2367391/150‑1 (MQ=255)
            aGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAgg                                                                                                                                          <  1:4999854/150‑1 (MQ=255)
            aGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAgg                                                                                                                                          <  1:495693/150‑1 (MQ=255)
               gTGGGGGCGATCGAGGTGCTTAAACGCCCTGGTTCGCTGGTGATGTCGCTGACAGTGGATCGGGTTGTATTAGGCATATGTTGTGCCCCAGCATAACGATGTTGCGGTGCCTGTCCGGGCGTATTCGCGGGCGTGCCCGCGAAGAGGccc                                                                                                                                       <  2:5702940/150‑1 (MQ=255)
                       gATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGc                                                                                                                               >  1:4229508/1‑150 (MQ=255)
                         tCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTa                                                                                                                             <  2:4591898/150‑1 (MQ=255)
                              gTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGAcgc                                                                                                                        >  2:3689272/1‑150 (MQ=255)
                                gCGTAAACGCCCTTGTTCGCTGGTGATGTCGAAGAATGTGTATCGGGTTGGATTCGGCATAAGTTGTGCTCAATAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGa                                                                                                                      <  1:7133193/150‑1 (MQ=38)
                                 cGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGaa                                                                                                                     >  1:5200667/1‑150 (MQ=255)
                                    aaaaGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACt                                                                                                                  >  2:1286822/1‑150 (MQ=255)
                                            ttGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATc                                                                                                          >  2:2256557/1‑150 (MQ=255)
                                             tGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCa                                                                                                         <  2:4866493/150‑1 (MQ=255)
                                               ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                       >  2:4625157/1‑150 (MQ=255)
                                                   cTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCct                                                                                                   <  2:507042/150‑1 (MQ=255)
                                                       tgaggccGAAGAACGTGGAGCGGGTTGGATCCGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCGGCCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctg                                                                                               <  1:5964571/144‑1 (MQ=38)
                                                        gATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTCCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                              <  2:1866336/150‑1 (MQ=255)
                                                            tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGt                                                                                          >  2:3814772/1‑150 (MQ=255)
                                                             ctaagaaTGTGGATGGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGGGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                         <  1:4666969/148‑1 (MQ=255)
                                                             cGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGTTGGTGCGGTGCCTGTCCGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                         <  1:3916206/150‑1 (MQ=255)
                                                             cGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                         <  1:5805004/150‑1 (MQ=255)
                                                              gaagacAGTGGCTCGGGTGGGATTCGGCATAAGTGGTGCTCAAGAAGACCGAGGGTGCGGTGCCGGTCTGGGCCTCTTCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTg                                                                                        <  1:4045872/150‑1 (MQ=37)
                                                              gaagaaAGTGTATCGGGTTGGATTAGGAATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTg                                                                                        <  1:6376773/150‑1 (MQ=255)
                                                              gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCAGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTg                                                                                        >  1:3203412/1‑150 (MQ=255)
                                                                  aaaGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGGTCCCCCGAAGGGGCCCTGACAGGCTAGACGCGAAACAGCGCCAACAGCGCCAGCGGCTGGTTGGccc                                                                                    >  1:858414/1‑149 (MQ=255)
                                                                     ggggATCGGGTTGGATGAGGCATCCGGTGTGCTCAAGAATACCGAGGGTGCGGTGCCTGTCAGGGCCTCGTCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGcc                                                                                 <  2:7299902/148‑1 (MQ=38)
                                                                     gTGGATCGGGTCGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGcc                                                                                 >  2:910046/1‑150 (MQ=255)
                                                                        gATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGt                                                                              <  2:4997971/150‑1 (MQ=255)
                                                                                   aTTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTgg                                                                   <  1:3610217/150‑1 (MQ=255)
                                                                                    ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                  <  2:2267187/150‑1 (MQ=255)
                                                                                    ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                  <  2:1991512/150‑1 (MQ=255)
                                                                                       ggCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGa                                                               >  1:4333614/1‑150 (MQ=255)
                                                                                        gCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGAt                                                              <  2:3419027/150‑1 (MQ=255)
                                                                                          aTAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCGTGGCTGATGc                                                            >  2:6497889/1‑150 (MQ=255)
                                                                                           tAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCg                                                           <  2:4167393/150‑1 (MQ=255)
                                                                                            aaGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGt                                                          <  1:5258618/150‑1 (MQ=255)
                                                                                                tgtgCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCg                                                      <  2:4283938/150‑1 (MQ=255)
                                                                                                     tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCa                                                 >  2:7044860/1‑150 (MQ=255)
                                                                                                     tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCa                                                 >  2:5422362/1‑150 (MQ=255)
                                                                                                     tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCa                                                 >  2:5972041/1‑150 (MQ=255)
                                                                                                      cAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCAt                                                <  2:5939587/150‑1 (MQ=255)
                                                                                                        aGAATACCGTTGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGGGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                              <  1:2132519/150‑1 (MQ=255)
                                                                                                        aGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                              <  1:5079391/150‑1 (MQ=255)
                                                                                                            tAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCt                                          <  2:2755000/150‑1 (MQ=255)
                                                                                                            aaaCGATGCTGCGGTGCCTGTCAGGGCCTCTTCGCGGCCCTGCCCCCCAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCAACAGCGCCTGCTGCTGCTTGGCCACCCGCTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCt                                          <  2:5222250/149‑1 (MQ=255)
                                                                                                                   ggtgcggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagc                                   <  1:243044/150‑1 (MQ=255)
                                                                                                                        ggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCCGTg                              >  2:2652159/1‑150 (MQ=255)
                                                                                                                             cTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCt                         <  1:6729064/150‑1 (MQ=255)
                                                                                                                              tGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTg                        >  2:3535924/1‑150 (MQ=255)
                                                                                                                               gTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGt                       <  1:2186484/150‑1 (MQ=255)
                                                                                                                               gTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGt                       <  1:7044860/150‑1 (MQ=255)
                                                                                                                                tctgggCCTCTTCGCGGGCGGGCCCGCTAAGATGCCCTGACAGGCTAGTCGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTc                      <  1:6866261/147‑1 (MQ=255)
                                                                                                                                 cAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCa                     <  2:7414522/150‑1 (MQ=255)
                                                                                                                                  aGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCAc                    >  2:3213829/1‑150 (MQ=255)
                                                                                                                                    ggCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATAGGCCTGCCCAGCCAGTGATTCTGTCACGt                  >  1:4589346/1‑150 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtc                 >  2:1225770/1‑150 (MQ=255)
                                                                                                                                       ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtctc               <  1:4782088/150‑1 (MQ=255)
                                                                                                                                       ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtctc               <  2:6193376/150‑1 (MQ=255)
                                                                                                                                       ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtctc               <  2:5200667/150‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:1946867/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGa                             <  1:1714332/135‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGa                             >  2:5485753/1‑135 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGa                             <  1:5485753/135‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:5833613/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:1343285/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg              >  2:3040438/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGaa            >  1:3738388/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCa                                                                                    >  2:4362358/1‑78 (MQ=255)
                                                                                                                                          ttCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGaa            >  1:4833678/1‑150 (MQ=255)
                                                                                                                                           tCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAAt           >  2:7110494/1‑150 (MQ=255)
                                                                                                                                           tCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAAt           >  2:3281190/1‑150 (MQ=255)
                                                                                                                                            cgcgGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATg          >  1:206034/1‑150 (MQ=255)
                                                                                                                                             gcgGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGc         >  1:3092889/1‑150 (MQ=255)
                                                                                                                                               gggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCg       <  1:1905259/150‑1 (MQ=255)
                                                                                                                                                ggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCgg      >  1:262226/1‑150 (MQ=255)
                                                                                                                                                ggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCgg      <  2:2868582/150‑1 (MQ=255)
                                                                                                                                                 gCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGc     >  1:2857384/1‑150 (MQ=255)
                                                                                                                                                  cGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGcc    <  1:3687132/150‑1 (MQ=255)
                                                                                                                                                   gTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGt                                                          >  2:4572530/1‑95 (MQ=255)
                                                                                                                                                    tGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAc  >  1:2494065/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCAC  >  NC_002947/1408398‑1408695

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: