Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,739,489 G→A P104S (CCA→TCA)  punR ← DNA‑binding transcriptional activator PunR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,739,4890GA98.2% 190.3 / NA 57P104S (CCA→TCA) punRDNA‑binding transcriptional activator PunR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/33);  minor base C (0/1);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACGGGCCATCCATCAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGGCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCC  >  NC_000913/1739243‑1739733
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
tGAACGGGCCATCCATCAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTg                                                                                                                                                                                                                                                   >  1:421987/1‑251 (MQ=255)
 gAACGGGCCATCCATCAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTg                                                                                                                                                                                                                                                   >  1:393522/1‑250 (MQ=255)
  aaCGGGCCATCCATCAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTTTACTTTATGACCTGAc                                                                                                                                                                                                                                                 >  2:75927/1‑251 (MQ=255)
     cgtccatccatCAACGCCAGCGGGTGGGGACGGGCAACAACGCAACTCCAGCTGAGCAGCCCCAGAGCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCCGTGCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGCCAAGAAGTTCGACATCAGCAAAATGGCGAGAAAAATCAACGCTCACCTGCAGTGTACGTTCTGACCGGACAat                                                                                                                                                                                                                                              <  2:322193/248‑1 (MQ=255)
          atccatcAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAACAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAAtattat                                                                                                                                                                                                                                         >  2:267061/1‑251 (MQ=255)
              atcAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCAGCTGCCGTGTACGTTCTGACCTGACAATATTATCCAc                                                                                                                                                                                                                                     <  1:461698/251‑1 (MQ=255)
                cAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACt                                                                                                                                                                                                                                    <  2:367950/250‑1 (MQ=255)
                 aaCGCCAGCGGGTGGTGACAGGCAACAACGCAACTCCAGCGTAGCATCCCCATAACCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCAGTCGGAAAGCGCATCCCAGACACCGCTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTGCTGACCTGACAATATTATCCACTGc                                                                                                                                                                                                                                  <  2:168574/251‑1 (MQ=255)
                     ccAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATa                                                                                                                                                                                                                              <  1:404200/251‑1 (MQ=255)
                        gCGGGGGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCt                                                                                                                                                                                                                           <  2:406979/251‑1 (MQ=255)
                            ggggtgACGGGCAACAACGCAACTCCAGCTTAGCATCCCCATACCCCGGAAGGCATAACGACCGCCGACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACt                                                                                                                                                                                                                       <  2:132041/249‑1 (MQ=255)
                               gtgACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGAc                                                                                                                                                                                                                    >  2:192398/1‑251 (MQ=255)
                                  aCTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACcgc                                                                                                                                                                                                                 <  1:95392/251‑1 (MQ=255)
                                     ggCAACAACGCAACTCCAGCTTAGCCTCCCCATATCCCGGAGGGCATCACGACCGCCTACCGGAATCGCGCGAGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGACGTTCGACATCATCAAAATGGCGATAAAACTCAACGATCATCTGCCGTGTCCGTTCTGACCTGACAATATTATCCACTGCGATCGCTAACTGACCGCGcc                                                                                                                                                                                                              >  2:31670/1‑251 (MQ=255)
                                         acaacGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGcc                                                                                                                                                                                                          >  2:157539/1‑251 (MQ=255)
                                         acaacGCAACTCCAGCTTAGCATCCCCATAACCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGcc                                                                                                                                                                                                          <  2:241856/251‑1 (MQ=255)
                                          caacGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCg                                                                                                                                                                                                         <  1:159373/251‑1 (MQ=255)
                                             cGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTc                                                                                                                                                                                                      <  1:307181/251‑1 (MQ=255)
                                               cAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGc                                                                                                                                                                                                    <  1:243417/251‑1 (MQ=255)
                                                  caccAGCCTAGCATACCCATAGCCCGGAAGGCATAACCACCGCCTACCGGCATCGCGCGTGTGGCGCCAATCGCCAGTTCCACGCGCCCGTCCGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCAACATCTTCAAAGTGGCGAAAAAAATCACCGATCAGCTGCCGTGTACGTTCGGACTTGACAATATTATCCACCGCGATAGCTACCTGAGCGTTCCTGCCGTTCGcaat                                                                                                                                                                                                 <  2:462700/249‑4 (MQ=255)
                                                      aGCTTAGCATCCCCATATCCCGGAAGGCATAACGCCCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTc                                                                                                                                                                                                      >  2:454201/1‑242 (MQ=255)
                                                      aGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTc                                                                                                                                                                                                      <  1:454201/242‑1 (MQ=255)
                                                                ccccATAGCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACGCCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAAGCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACGGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCt                                                                                                                                                                                   <  1:220781/251‑1 (MQ=255)
                                                                    atatCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCg                                                                                                                                                                               >  1:119425/1‑251 (MQ=255)
                                                                    atatCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCg                                                                                                                                                                               <  1:98180/251‑1 (MQ=255)
                                                                      atCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCgg                                                                                                                                                                              <  1:173007/250‑1 (MQ=255)
                                                                        cccGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCGATCTGCTGACATTGCTGGCGGGt                                                                                                                                                                            <  1:109172/250‑1 (MQ=255)
                                                                            gAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGGTCt                                                                                                                                                                       >  1:437537/1‑251 (MQ=255)
                                                                               ggCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTCTCTGCTGACATTGCTGGCGGGGGTTCTGcc                                                                                                                                                                    >  2:11950/1‑250 (MQ=255)
                                                                                 cATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCAtt                                                                                                                                                                  <  1:157539/251‑1 (MQ=255)
                                                                                 cATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCAtt                                                                                                                                                                  <  1:186748/251‑1 (MQ=255)
                                                                                  aTAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCAttt                                                                                                                                                                 <  1:434198/251‑1 (MQ=255)
                                                                                     aCGACCGCCGACCGGAAGCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGTCATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCAtttttt                                                                                                                                                              <  1:27505/251‑1 (MQ=255)
                                                                                           gCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCg                                                                                                                                                                                                     <  1:372877/206‑1 (MQ=255)
                                                                                           gCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCg                                                                                                                                                                                                     >  2:372877/1‑206 (MQ=255)
                                                                                             cTACCGGAAGCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTTATCTGCCTTTTTTTGtcaac                                                                                                                                                       >  2:257690/1‑250 (MQ=255)
                                                                                                       cgcgcgTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTg                                                                                                                                                                                     >  2:279653/1‑210 (MQ=255)
                                                                                                       cgcgcgTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTg                                                                                                                                                                                     <  1:279653/210‑1 (MQ=255)
                                                                                                       cgcgcgTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCAGCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCGGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCctt                                                                                                                                            <  2:464236/251‑1 (MQ=255)
                                                                                                         cgcgTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCcttct                                                                                                                                          <  1:264186/251‑1 (MQ=255)
                                                                                                         cgcgTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTGCCTGAAAGACAAGAAGTGCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATGGCTGGCTGGTGATCTGCATTTTTTTGACAACAGAGCGCCcttct                                                                                                                                          <  1:458913/251‑1 (MQ=255)
                                                                                                                 gcgcCAATCGCCAGTGCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTTCATTTTTTTGACAACAGAGCGCCCGTCTTTGAGaa                                                                                                                                   <  2:421987/250‑1 (MQ=255)
                                                                                                                    ccAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCa                                                                                                                               <  2:324265/251‑1 (MQ=255)
                                                                                                                    ccAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACGTCCTGAAAGACAAGAAGGTCGTCATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTATCTGACCGCGCCATCCGTTCGCTATCTGCTGACATTGCTGGCTGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCa                                                                                                                               <  1:226057/251‑1 (MQ=255)
                                                                                                                     cAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCAc                                                                                                                              <  1:365824/251‑1 (MQ=255)
                                                                                                                                                    cgcATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGt                                                                                               <  2:119425/251‑1 (MQ=255)
                                                                                                                                                      cATCCCAGACCCCGTTGAACCCTTCCTGAAAGACAAGAAGTTCGACATCAGCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAAGAGGATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTGCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTGTGACAACAGAGCGCCCTTCTTTGAGAAACCGCGCTCCAGCAGCGGTCAGCTCCACATCACGGTGc                                                                                             <  2:437537/251‑1 (MQ=255)
                                                                                                                                                                      gAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCATTTTTTTGTCAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCCCATCACGGTGCCGTCGTTCCAagaacg                                                                             >  2:146515/1‑248 (MQ=255)
                                                                                                                                                                       aaCACTTCCTGAAAGACAAGAAGTTCGACATCCTCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGGTCTGCATTTTTTTGACAACAGAGCGCCCGTCTTTGAGAAACAACGCCCCCGCAGCGGTCAGCTCCCCATCACGGTGCCGTCGttacaagggcgg                                                                            >  2:482641/1‑241 (MQ=255)
                                                                                                                                                                                tGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGATCTGCCTTTTTTTGACAACAGAACGCCCTTCTTTGAGAAACAACGCTCCAGCAGCGGGCCGCGCCACATCACGGTGCCGCCGTTCAAAGAGCGGCACCgcc                                                                     >  1:314713/1‑249 (MQ=255)
                                                                                                                                                                                          aaagTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGCCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATGTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCAGCGCCAGCCACCCTTcc                                                        <  1:332970/250‑1 (MQ=255)
                                                                                                                                                                                                  aCATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGGGGTCTGCATTTTTTTGCCAACAGCGCGCCCCTCTTTGAGAAACCACGCCCCCGCAGCGGGCCGCTCCACATCACGGGGCCGCCGTTCAAAAAGCGGGACCcccagccccccctcccaccgccg                                                 >  2:396187/1‑235 (MQ=255)
                                                                                                                                                                                                   catcatCAAAATGGCGATAACAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGTGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGTTCCGTCGTTCAAAGAGCGGCCCCGCCCGCCACCCTTTCCGCTGACGc                                               >  1:92799/1‑251 (MQ=255)
                                                                                                                                                                                                   catcatCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGGTCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCCCGGGGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGCCGc                                               >  1:393228/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccGTGTACGTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAACGCGCCCTTCTTTGAGAAACCACGCTCCAGCCGCGGTCAGCCACACCTCACGGTGCCGCCGGTCAAAGAACGGCACCGCCAGCCCACTCTCCAGCTGACGCCCGCTATAGCTGACCCCAGAAGGGACGCCATACAGCt          >  2:324295/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGcc  >  1:120830/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gTTCTGACCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCCCCTGTGcc  >  1:19057/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
TGAACGGGCCATCCATCAACGCCAGCGGGTGGTGACTGGCAACAACGCAACTCCAGCTTAGCATCCCCATATCCCGGAAGGCATAACGACCGCCTACCGGAATCGCGCGTGTTGCGCCAATCGCCAGTTCCACGCGCCCGTCGGAAAGCGCATCCCAGACACCGTTGAACACTTCCTGAAAGACAAGAAGTTCGACATCATCAAAATGGCGATAAAAATCAACGATCATCTGCCGTGTACGTTCTGGCCTGACAATATTATCCACTGCGATAGCTAACTGACCGCGCCAGCCGTTCGCTATCTGCTGACATTGCTGGCGGGTGATCTGCATTTTTTTGACAACAGAGCGCCCTTCTTTGAGAAACCACGCTCCAGCAGCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCC  >  NC_000913/1739243‑1739733

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