Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,956 T→C L117L (CTA→CTG arcA ← DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,9560TC91.0% 171.6 / 6.9 68L117L (CTA→CTGarcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (40/21);  total (43/25)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

CGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCG  >  NC_000913/4639711‑4640199
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGtt                                                                                                                                                                                                                                                  >  2:222959/1‑251 (MQ=255)
              cAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATc                                                                                                                                                                                                                                     <  2:45983/250‑1 (MQ=255)
                      ggcaggTCATTTTATTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGAACTGCTCGCCATCAGGGCCGATCAACAAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTAAACGCTACGAAGTTCTTCGGTGACCGTACCCAGATTAATGGTCCGGGACAGCAGGTTGCGTGCACGAATCGTCCGTTc                                                                                                                                                                                                                             >  1:355555/5‑250 (MQ=255)
                             cATTTTCTTCAACCGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCTCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCCACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGCTTCATGGTACGGGACAGACGGTTGCGAGCACGAATTGTCAATTCAcgcaggg                                                                                                                                                                                                                     >  2:71947/1‑247 (MQ=255)
                              aTTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGtt                                                                                                                                                                                                                    <  1:222959/251‑1 (MQ=255)
                                      cagcagTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCATGCGGGTTGAACGGtt                                                                                                                                                                                                            <  2:77227/251‑1 (MQ=255)
                                      cagcagTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGtt                                                                                                                                                                                                            >  1:470228/1‑251 (MQ=255)
                                           gTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTg                                                                                                                                                                                                          >  1:42445/1‑248 (MQ=255)
                                           gTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTg                                                                                                                                                                                                          <  2:42445/248‑1 (MQ=255)
                                                   cgagACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTc                                                                                                                                                                                               >  1:32213/4‑251 (MQ=255)
                                                      gACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACGTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATc                                                                                                                                                                                            <  2:455006/251‑1 (MQ=255)
                                                      gACTGAATTGTGCCTGGGTTTTCACAGTAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATc                                                                                                                                                                                            <  2:321264/251‑1 (MQ=255)
                                                              ttggcgggggTTTTTACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTGGTACTGCTCGCCATAAGGGCCGATCAACGAACGGCTGTTGATGGCCAGTTCCCAACCATGGAACTTGTAGCTTGCAACGCTAAGACGTTCGTCGCTGACAGTACCCAGATTCATGTTACGGGACAGCAGGGTGCGTGAACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCg                                                                                                                                                                                     <  2:336149/243‑1 (MQ=255)
                                                                ggCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCACCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATTGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGAt                                                                                                                                                                                   <  1:71947/249‑1 (MQ=255)
                                                                    tGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGCTTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGCGATGTAGTCATCTGCACCGCTTTCgg                                                                                                                                                                              >  1:482190/1‑250 (MQ=255)
                                                                         tttCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGcc                                                                                                                                                                          <  2:344047/250‑1 (MQ=255)
                                                                             acaGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCgaga                                                                                                                                                                      >  1:107873/1‑250 (MQ=255)
                                                                                gAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTTATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAAttt                                                                                                                                                                  <  2:406190/251‑1 (MQ=255)
                                                                                  aGTGAAGCATGGCGCGGAACTCGCTGAGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGCTGTCCAGTTCCCAACCATTGAACTTGTAGCTATCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGGAAGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGGTTGGTGCTGTAGTCATCTGCACCGATTTCGATGCCGAGAAttttt                                                                                                                                                                >  2:232361/1‑250 (MQ=255)
                                                                                   gTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTAt                                                                                                                                                               >  1:67848/1‑251 (MQ=255)
                                                                                         cATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACtt                                                                                                                                                         <  2:253364/251‑1 (MQ=255)
                                                                                         cATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCAAGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACtt                                                                                                                                                         >  1:188730/1‑251 (MQ=255)
                                                                                                  gAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCAc                                                                                                                                                >  1:153051/1‑251 (MQ=255)
                                                                                                     cccGCTGCTCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGGTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTGCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGc                                                                                                                                             <  2:32213/249‑1 (MQ=255)
                                                                                                     cTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCTGGGCAGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGGGTGCACGAATCGTCAGTTAACGGGGGTTGAACGGTTTGGTGGTGTAGTCATCTGCACCGATTTCGAGGCCGAGACTTTTATCGACTTCGTTGTCACGGc                                                                                                                                             >  2:205841/1‑251 (MQ=255)
                                                                                                     cTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGc                                                                                                                                             >  2:64493/1‑251 (MQ=255)
                                                                                                         cTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGCTCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGAAGCTTTCAACGCTACGACGTTCTTCCCTGAAAGAACCCAGATTCAAGGTACGGGACAGCCGGTTGCGTGCACGAAACGTCAGTTCACGCGGGTTAAACGGTTTGGTGCTGGATGCATCTACAACGATTTCGAGGCCGAGATTTTGATCGACTTCGTTGTCACGGccc                                                                                                                                           >  2:410323/1‑248 (MQ=255)
                                                                                                         cTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGACCGGCTGGTGATGTCCAGTTACCAACCATTGAACTTGTAGCTTTCAACGCTACGACGATCTTCGCTGACAGTACCCAGATTCATGGGACGGGACAGCAGGATGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGTGAATTTTATAGACTTCGATGTCACGGcc                                                                                                                                            >  2:410330/1‑248 (MQ=255)
                                                                                                            cgcgGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTc                                                                                                                                        >  1:86749/1‑249 (MQ=255)
                                                                                                            cgcgGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGAAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGACACGGCCAGTc                                                                                                                                        >  1:479728/1‑249 (MQ=255)
                                                                                                              cgGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACTGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGa                                                                                                                                    <  2:213409/251‑1 (MQ=255)
                                                                                                                gCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAAc                                                                                                                                  >  1:5719/1‑251 (MQ=255)
                                                                                                                       gTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACg                                                                                                                           >  1:123632/1‑251 (MQ=255)
                                                                                                                        tACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACCTCAACGc                                                                                                                          >  2:117801/1‑251 (MQ=255)
                                                                                                                               cGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTc                                                                                                                  <  1:41848/251‑1 (MQ=255)
                                                                                                                                  cATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCa                                                                                                                                                                                                                            >  2:155160/1‑143 (MQ=255)
                                                                                                                                  cATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCa                                                                                                                                                                                                                            <  1:155160/143‑1 (MQ=255)
                                                                                                                                   aTCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCt                                                                                                              >  1:273628/1‑251 (MQ=255)
                                                                                                                                   aTCAGGGCCGATCAACGAACGGCTGGTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGCACAGCAACGccaacttcgccc                                                                                                               >  2:326205/1‑240 (MQ=255)
                                                                                                                                            gATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCgcgcc                                                                                                     >  1:457184/1‑250 (MQ=255)
                                                                                                                                             aTCAACGAACGGCTGTGGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACAGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAg                                                                                                    <  1:457157/251‑1 (MQ=255)
                                                                                                                                                   gaacgGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTCCGACGTTCTTCGCTGCCAGTACCCAGCTTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGCTTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCCCCTTCTCGCGCCGTTCAcc                                                                                              >  2:470177/1‑250 (MQ=255)
                                                                                                                                                         cTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTa                                                                                         >  1:552/1‑250 (MQ=255)
                                                                                                                                                                     aGTTCCCAACCATTGTACTTGTGGCTTTCAACGCTACGACATTCTTCGCTTACATTAACCAGATTCATGGTACGGGACAGCCGGTTGCGTGCACGAATCGTCCGTTCACGCGGGTTGAACGTTTTGGTGCTGTAGACATCTGCACCGATTTCGAGTCCGCGAATTTTATCGCCTTCGTTGTCACGGCCAGTCAGGAACATCAACCACACATTCGCCCTCTCTCCCAGTTCACTCCCTAACAGAAGACCGtt                                                                            >  2:267826/1‑251 (MQ=255)
                                                                                                                                                                         cccAACCTATGAACTTGTAGCTTTCAACGCTTCGACGTTTGTCGCTGACAGTACCCAGGTGCATGGTACGTGACAGCAGGTTGCGTGCACGAATCGGCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGGACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCATGCAGGAACGTCAACGCAACATTCGCCTGCTCGCGCAGTTCTCGCGCTAACAGAAGACCGTTCTTa                                                                        <  2:456014/251‑1 (MQ=255)
                                                                                                                                                                             aCCATTGAACTTGTAGCTTTCAACGCTCCGACGTTCTTCGCTGACCGTACCCAGATTCATGGTACGGGACAGCAGGGTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCCTCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCa                                                                                                                                                 >  1:160409/1‑175 (MQ=255)
                                                                                                                                                                             aCCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCa                                                                                                                                                 <  2:160409/175‑1 (MQ=255)
                                                                                                                                                                             aCCATTGAACAAGTAGCGTTCAACGCAACGACGGTCTTAGCTGACAGTACCCAGATTCCTGGTACGGGTCAGCAGGATGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGACTGCGAGGCCGAGAGTTTTATCGACTTCGTTTTCACGGCCATTCAGGAACACCAACGCAAGATACGCCTGCTCTCGCAGTTCACGCGTTCACCGAAGACCGGTCTTACCCg                                                                    >  2:286985/1‑251 (MQ=255)
                                                                                                                                                                                  tGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACt                                                                                                                                                          >  1:364023/1‑161 (MQ=255)
                                                                                                                                                                                  tGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGGTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACt                                                                                                                                                          <  2:364023/161‑1 (MQ=255)
                                                                                                                                                                                  tGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGAACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAg                                                                <  1:412443/250‑1 (MQ=255)
                                                                                                                                                                                          tAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGGTATc                                                       >  2:92889/1‑251 (MQ=255)
                                                                                                                                                                                            gCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATcc                                                      >  1:6660/1‑250 (MQ=255)
                                                                                                                                                                                               ttCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATg                                                   >  1:409768/1‑250 (MQ=255)
                                                                                                                                                                                                tCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGa                                                  >  1:328864/1‑250 (MQ=255)
                                                                                                                                                                                                tCAACGCCACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGa                                                  <  2:328864/250‑1 (MQ=255)
                                                                                                                                                                                                  aaCGCGACGACGCTCTTCGCTGACAGCACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTTGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCa                                               <  2:457184/251‑1 (MQ=255)
                                                                                                                                                                                                         cgacgTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACa                                                                                                                                 <  2:84161/163‑1 (MQ=255)
                                                                                                                                                                                                         cgacgTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACa                                                                                                                                 >  1:84161/1‑163 (MQ=255)
                                                                                                                                                                                                               tcttcGCTGACAGTACCCAGATTCATGGTACGGGCCAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGAAATCCATGATCACCAGGTTGATGTc                                  >  1:369119/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttcGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCat                                >  1:42076/1‑251 (MQ=255)
                                                                                                                                                                                                                       gACAGTACCCAGATTCATGGTACGGGACAGCAGGGTGCGTGCAAGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTGGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTc                            <  2:17192/249‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGCTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGAt             >  1:470578/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGc                                                                                                                                             >  2:353875/1‑123 (MQ=255)
                                                                                                                                                                                                                                     tCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGc                                                                                                                                             <  1:353875/123‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCAttt       >  1:91821/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATTATCCCCAGGTTGATGTCATATTCAAAGAGGCTCTGAAGCCTTTc      >  1:322389/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gggACAGCAGGATGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATAGACTTCGTTGTCACGGCCAGTCAGGTACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAAAAGCAGACCGTTCTTACCCGGCAGATTGATATACATGATCACCAGGTTGATGTCATATTCAGAGAGATACAGATGCCTTTCcgcg  >  2:31616/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGTGATGTAGTCATCTGCACCGATTTCGAGGCCGAGAATTTTATCGACTTCGTTGTCACGGCCAGTCAGGAACATCAACGCAACATTCGCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCG  >  NC_000913/4639711‑4640199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: