Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,403,596 G→A 14.1% intergenic (‑56/+25) 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,5960GA14.1% 272.3 / 32.3 114intergenic (‑56/+25)PP_1227/PP_1228membrane protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (48/49);  new base A (9/7);  total (57/57)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 117 total aligned reads displayed.

TGGCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCA  >  NC_002947/1403450‑1403745
                                                                                                                                                  |                                                                                                                                                     
tggCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGc                                                                                                                                                    <  2:1353509/150‑1 (MQ=255)
    gATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgc                                                                                                                                                >  1:1400021/1‑150 (MQ=255)
      tGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGa                                                                                                                                              <  1:692206/150‑1 (MQ=255)
       gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGACACGCCCGCGaa                                                                                                                                             >  1:617856/1‑150 (MQ=255)
       gAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGACACGCCCGCGaa                                                                                                                                             <  2:617856/150‑1 (MQ=255)
         aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGTCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                           <  2:671098/150‑1 (MQ=255)
          cAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAagag                                                                                                                                          >  1:761164/1‑150 (MQ=255)
            gggTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                             <  1:2152399/145‑1 (MQ=255)
            gggTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                             >  2:2152399/1‑145 (MQ=255)
             ggTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGcc                                                                                                                                       <  2:1785882/150‑1 (MQ=255)
               tGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAagag                                                                                                                                          <  1:2044819/145‑1 (MQ=255)
               tGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAagag                                                                                                                                          >  2:2044819/1‑145 (MQ=255)
                    gcgATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAg                                                                                                                                <  2:292348/150‑1 (MQ=255)
                                  tAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAAc                                                                                                                  <  1:1635976/150‑1 (MQ=255)
                                       gCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTc                                                                                                             >  1:2136525/1‑150 (MQ=255)
                                         ccTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCa                                                                                                           >  2:1316312/1‑150 (MQ=255)
                                              ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                      <  1:227415/150‑1 (MQ=255)
                                              ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                      <  1:2154722/150‑1 (MQ=255)
                                              ttCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                      >  2:2154722/1‑150 (MQ=255)
                                                   tGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                           <  2:493390/108‑1 (MQ=255)
                                                   tGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAaga                                                                                                                                           >  1:493390/1‑108 (MQ=255)
                                                    ggTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCAt                                                                                                          >  1:1459707/1‑140 (MQ=255)
                                                       gATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                             <  2:1400021/150‑1 (MQ=255)
                                                       gATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                             >  1:618709/1‑150 (MQ=255)
                                                          gTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTgg                                                                                          <  2:1920131/150‑1 (MQ=255)
                                                           tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGaa                                                                                                                                             >  2:1620466/1‑98 (MQ=255)
                                                           tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACa                                                                                                                                 >  1:1076187/1‑110 (MQ=255)
                                                           tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACa                                                                                                                                 <  2:1076187/110‑1 (MQ=255)
                                                            cGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTgg                                                                                          <  1:496676/148‑1 (MQ=255)
                                                            cGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTgg                                                                                          >  2:496676/1‑148 (MQ=255)
                                                            cGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                        <  2:1303219/150‑1 (MQ=255)
                                                            cGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                        >  2:1634671/1‑150 (MQ=255)
                                                             gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGaaa                                                                                                                   >  2:1921008/1‑122 (MQ=255)
                                                             gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGaaa                                                                                                                   <  1:1921008/122‑1 (MQ=255)
                                                             gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTg                                                                                       >  1:2146462/1‑150 (MQ=255)
                                                             gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                           >  2:1063271/1‑98 (MQ=255)
                                                             gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                           <  1:1063271/98‑1 (MQ=255)
                                                              aagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGAc                                                                                                                                  <  1:722261/106‑1 (MQ=255)
                                                              aagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGAc                                                                                                                                  >  2:722261/1‑106 (MQ=255)
                                                               agaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGc                                                                                     <  1:1706175/150‑1 (MQ=255)
                                                                 aaaGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCa                                                                                   >  2:1036367/1‑150 (MQ=255)
                                                                       gATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgct                                                                                            >  1:36493/1‑135 (MQ=255)
                                                                       gATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgct                                                                                            <  2:36493/135‑1 (MQ=255)
                                                                           gggTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAg                                                                         <  1:1634671/150‑1 (MQ=255)
                                                                              ttGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGc                                                                      <  2:1106928/150‑1 (MQ=255)
                                                                              ttGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTCCTGCTGGTTGGCCAGCCGGTTCAGGGc                                                                      >  1:1708806/1‑150 (MQ=255)
                                                                                   ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                 <  2:761164/150‑1 (MQ=255)
                                                                                   ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                 <  1:1316312/150‑1 (MQ=255)
                                                                                    tAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGaaa                                                                                                                   <  2:393577/99‑1 (MQ=255)
                                                                                    tAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCt                                                                <  2:2146462/150‑1 (MQ=255)
                                                                                      ggCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGa                                                              >  2:898243/1‑150 (MQ=255)
                                                                                       gCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGAt                                                             >  2:690601/1‑150 (MQ=255)
                                                                                         aTAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTg                                                                   >  2:1719263/1‑142 (MQ=255)
                                                                                              ttGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                        >  2:2091124/1‑116 (MQ=255)
                                                                                              ttGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGtt                                                                                        <  1:2091124/116‑1 (MQ=255)
                                                                                                    tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGc                                                                                 <  1:777687/117‑1 (MQ=255)
                                                                                                    tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGc                                                                                 >  2:777687/1‑117 (MQ=255)
                                                                                                    tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                 <  2:1169868/133‑1 (MQ=255)
                                                                                                    tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGc                                                                 >  1:1169868/1‑133 (MQ=255)
                                                                                                    tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCa                                                >  1:362671/1‑150 (MQ=255)
                                                                                                       aGAATAACGATGGTGCGCTGCCCGTCACGCCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                             <  2:1708806/150‑1 (MQ=255)
                                                                                                         aaTAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCAGc                                           >  1:251327/1‑150 (MQ=255)
                                                                                                         aaTAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCCCATGGTTGGCGACACGGCGCAGGTCGGGGCTGATTCGTGCCGGTTCATCGGc                                           >  2:612759/1‑150 (MQ=37)
                                                                                                            aaCGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCa                                                                                   <  1:1564853/107‑1 (MQ=255)
                                                                                                            aaCGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTGGGCCa                                                                                   >  2:1564853/1‑107 (MQ=255)
                                                                                                               gATGGTGCGGTGCCTGTCAGGACCTCTTCGCGGGCGTGCCCGCGCAGAAGCGCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGccc                                     >  2:1495493/1‑150 (MQ=255)
                                                                                                                    tgcggtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCa                                                                                                        <  1:962059/78‑1 (MQ=255)
                                                                                                                     gcggtgcCTGGCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagccag                               <  2:1637125/150‑1 (MQ=255)
                                                                                                                            cTGTCAGGGCCTCTTCGCTGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagcc                                 <  2:2198406/141‑1 (MQ=255)
                                                                                                                            cTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCt                        >  1:292277/1‑150 (MQ=255)
                                                                                                                             tggcAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGAt                                                   >  1:2331407/4‑122 (MQ=255)
                                                                                                                             tGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGGTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCGGTGATTCTg                       >  1:1278701/1‑150 (MQ=255)
                                                                                                                             tGTCAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGAt                                                   <  2:2331407/122‑1 (MQ=255)
                                                                                                                              tctaGGGCCTCCGCGCGGGCGTGCCCGCGCAGTGGCACTGACAGGCTAGACGCGCATGCGCGCCATCAGCGCCTGCTGCTGGTGGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGt                      <  1:1780162/147‑1 (MQ=255)
                                                                                                                                cAGGGCCGCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCa                    <  1:1036367/150‑1 (MQ=255)
                                                                                                                                 aGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGc                                                       >  1:465946/1‑114 (MQ=255)
                                                                                                                                 aGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGc                                                       <  2:465946/114‑1 (MQ=255)
                                                                                                                                  gggCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACg                  >  2:1253362/1‑150 (MQ=255)
                                                                                                                                  gggCCCCTTCGCGGGCGTGCCCAAGAAGAGAAACTGACAGACCAGACGCGAAACTGCTGCATGAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTCGCTGAGCCGTGCCGGTTCAACGGCGTGGCCAGCCAGTGATTCTGTccga                  >  2:1780162/1‑147 (MQ=255)
                                                                                                                                   ggCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTc                     >  1:1631703/1‑146 (MQ=255)
                                                                                                                                   ggCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTc                     <  2:1631703/146‑1 (MQ=255)
                                                                                                                                   ggCCTCTTCGCGGGAGTGCCAGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCCGCCCGTGATTCTGTCACGt                 >  2:855337/1‑150 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTcc                                                                                                            >  1:827129/1‑57 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTcc                                                                                                            <  2:827129/57‑1 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGACAAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtct               <  1:612759/150‑1 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATc                                              >  1:1750467/1‑119 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATc                                              <  2:1750467/119‑1 (MQ=255)
                                                                                                                                      ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCa                                                                          <  1:1851318/90‑1 (MQ=255)
                                                                                                                                      ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCa                                                                          >  2:1851318/1‑90 (MQ=255)
                                                                                                                                      ctctTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAACGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGTCCTACCCAGCCAGTGATTCTGTCACGtctc              >  2:205013/1‑150 (MQ=255)
                                                                                                                                         ttCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCCGCCAGTGATTCTGTCACGGCTCGaa           >  1:1864374/1‑150 (MQ=255)
                                                                                                                                         ttCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGaa           >  2:553995/1‑150 (MQ=255)
                                                                                                                                             cgGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtc                >  1:1322126/1‑141 (MQ=255)
                                                                                                                                             cgGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtc                <  2:1322126/141‑1 (MQ=255)
                                                                                                                                              gggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTgg                                                                  >  1:1012621/1‑90 (MQ=255)
                                                                                                                                              gggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCg      <  2:618709/150‑1 (MQ=255)
                                                                                                                                              gggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCg      <  1:1523923/150‑1 (MQ=255)
                                                                                                                                                  gTGCCCGCGAGGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCa  <  1:1710132/150‑1 (MQ=255)
                                                                                                                                                  gTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGAtt                                                  <  1:1919302/102‑1 (MQ=255)
                                                                                                                                                  gTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGAtt                                                  >  2:1919302/1‑102 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGGCGATGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCCGGCCA  >  NC_002947/1403450‑1403745

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