Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,559,766 A→C Q500P (CAG→CCG)  rtcR → sigma 54‑dependent transcriptional regulator of rtcBA expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,559,7660AC98.4% 227.8 / ‑3.2 62Q500P (CAG→CCG) rtcRsigma 54‑dependent transcriptional regulator of rtcBA expression
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (28/33);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAAC  >  NC_000913/3559539‑3559994
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
cTTTGCCACTAGCGGACGCAGCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAg                                                                                                                                                                                                                               <  1:263499/236‑1 (MQ=255)
  ttGCAAATAGCGAACGCAGACCTCTGGACGTGGTTGAAGATGAGATCAACCGTCTGAGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGcc                                                                                                                                                                                                                             <  2:256138/236‑1 (MQ=255)
             gaacgcagcgatctGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACGGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGa                                                                                                                                                                                                                  <  2:766113/225‑1 (MQ=255)
                     cACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCCCGCCCGGGCACAGCCACCGTCAATGACGCGGa                                                                                                                                                                                                            >  1:220947/1‑234 (MQ=255)
                         cggtaCGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGAt                                                                                                                                                                                                           <  2:676318/227‑1 (MQ=255)
                           ggACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTAc                                                                                                                                                                                                    >  1:545677/1‑236 (MQ=255)
                            gACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACg                                                                                                                                                                                                   <  1:410863/236‑1 (MQ=255)
                                                  aCCGTCTGCGCTATAACTGGCAGGAGAGCCGGCCCGACGCGCTTACGGCGTTGCTGGGCACAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGt                                                                                                                                                                             <  2:475174/236‑1 (MQ=255)
                                                  aCCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTGTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAGATCACTGGCGCGgtt                                                                                                                                                                                 >  2:457743/1‑232 (MQ=255)
                                                  aCCGTCGGCGGTATAACGGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGt                                                                                                                                                                             <  2:56905/236‑1 (MQ=255)
                                                      tCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCTGGCAAAGCCACAGTCAACGACCCGGCTCGGCTACGCAACTACCTGGTGCGTTTTGGGCTGa                                                                                                                                                                         >  1:78532/1‑236 (MQ=255)
                                                        tGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCTGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGTACGGTTACGCAAATACCTGGCCCGGTTTGGTCttccg                                                                                                                                                                       >  2:693642/1‑232 (MQ=255)
                                                                       aGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGc                                                                                                                                                               >  1:162521/1‑229 (MQ=255)
                                                                       aGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGc                                                                                                                                                               <  2:162521/229‑1 (MQ=255)
                                                                         ggaagCCGCACCTCCGCGCTTACGGCGTTGCTGGGAGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACCCGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTGCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          <  2:545677/229‑1 (MQ=255)
                                                                         gagaGCCGCCCCTCCGCGATTACGGCGTTGCGGGGCACAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          <  2:222627/232‑1 (MQ=255)
                                                                         gagaGCCGACCCTCCGCGCTTACGGCGTGGCCGGGCGCAGAGGCAGAAAACATCGATCTCTGCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          <  1:709101/232‑1 (MQ=255)
                                                                         gacacccgccCCTCCGCCGTTACGGCGTTGCTGGGCGCAGACGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          <  2:137778/227‑1 (MQ=255)
                                                                                  cccTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGc                                                                                                                                                 <  1:457743/232‑1 (MQ=255)
                                                                                      ccGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTc                                                                                                                                          <  1:502397/235‑1 (MQ=255)
                                                                                      ccGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTc                                                                                                                                          >  2:502397/1‑235 (MQ=255)
                                                                                      ccGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCCCAGCTcc                                                                                                                                         >  1:727760/1‑236 (MQ=255)
                                                                                        gcgcTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAg                                                                                                                                       <  1:752624/236‑1 (MQ=255)
                                                                                               cGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGc                                                                                                                                                                                                       <  2:412249/165‑1 (MQ=255)
                                                                                               cGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGc                                                                                                                                                                                                       >  1:412249/1‑165 (MQ=255)
                                                                                                          ggCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCa                                                                                                                                        <  1:578321/216‑1 (MQ=255)
                                                                                                                 agGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGGTAGCaac                                                                                                              >  2:72462/1‑234 (MQ=255)
                                                                                                                 agGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGACCAGGTCCAGTTGaaa                                                                                                                                 >  2:578321/1‑215 (MQ=255)
                                                                                                                    cAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          <  1:469048/189‑1 (MQ=255)
                                                                                                                    cAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          >  2:469048/1‑189 (MQ=255)
                                                                                                                      gAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCa                                                                                                          >  1:726301/1‑235 (MQ=255)
                                                                                                                      gAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGAAGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGCGGGAAGCCGGGCAGGATCAGCCCAGCTCCAGTTGAATATGGTGGTCCGTCAGTACCTGCa                                                                                                          >  2:202973/1‑235 (MQ=255)
                                                                                                                           cATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACg                                                                                                     >  1:51354/1‑235 (MQ=255)
                                                                                                                              cGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACg                                                                                                     >  1:41543/1‑232 (MQ=255)
                                                                                                                                      tCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGACTATGGTGGTCCGTCAGCACCTGCATCACGCTTACTggcggc                                                                                         >  2:465371/1‑236 (MQ=255)
                                                                                                                                       cGACCGTATGCAACTGGAACTCGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGgcggcg                                                                                        <  1:89421/236‑1 (MQ=255)
                                                                                                                                             tATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGgc                                                                                                                                                                                       <  1:344889/135‑1 (MQ=255)
                                                                                                                                             tATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGgc                                                                                                                                                                                       >  2:344889/1‑135 (MQ=255)
                                                                                                                                             tATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCgg                                                                                  >  1:599258/1‑235 (MQ=255)
                                                                                                                                               tGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCCGAACCTGCATCCCGCTTACTGGCGACTGTGCGTct                                                                                >  2:399359/1‑235 (MQ=255)
                                                                                                                                                     tGGAACGCGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAg                                                                          <  1:631094/236‑1 (MQ=255)
                                                                                                                                                       gAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTc                                                                                                                    <  1:461498/192‑1 (MQ=255)
                                                                                                                                                       gAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTc                                                                                                                    >  2:461498/1‑192 (MQ=255)
                                                                                                                                                            cGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATAc                                                                                                                                                                                            >  2:656446/1‑115 (MQ=255)
                                                                                                                                                            cGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATAc                                                                                                                                                                                            <  1:656446/115‑1 (MQ=255)
                                                                                                                                                             gTTAGCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATc                                                                  <  1:542329/236‑1 (MQ=255)
                                                                                                                                                                 tCGCTATCGGCAGCAAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCAt                                                                   <  2:41543/231‑1 (MQ=255)
                                                                                                                                                                       tCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCatgct                                                         >  1:250126/1‑235 (MQ=255)
                                                                                                                                                                       tCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGCCGGTGTAGACGGCATCTACCatgct                                                         >  1:440552/1‑235 (MQ=255)
                                                                                                                                                                       tCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGAGGCGGGGTAGAAGGCATCTACCatgct                                                         >  1:563511/1‑235 (MQ=255)
                                                                                                                                                                             gccAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCatgctgatgctg                                                  <  2:250126/236‑1 (MQ=255)
                                                                                                                                                                                aGGCAAAGTCGCTTGCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGccc                                               <  1:588530/236‑1 (MQ=255)
                                                                                                                                                                                      aGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCatat                                           <  1:693642/234‑1 (MQ=255)
                                                                                                                                                                                            tttCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGa                                       <  1:72462/231‑1 (MQ=255)
                                                                                                                                                                                                     ccGGACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATTTACCATGCTGATTCTGTCCCTATTGACCATTGCGTTACgg                         >  2:296239/1‑236 (MQ=255)
                                                                                                                                                                                                        gACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCg                                                                                   <  2:126097/176‑1 (MQ=255)
                                                                                                                                                                                                        gACGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCg                                                                                   >  1:126097/1‑176 (MQ=255)
                                                                                                                                                                                                          cGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          >  1:266156/1‑103 (MQ=255)
                                                                                                                                                                                                          cGCCAGCTTTTTGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGc                                                                                                                                                          <  2:266156/103‑1 (MQ=255)
                                                                                                                                                                                                                  ttttGACGTTTCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGt            <  1:296239/236‑1 (MQ=255)
                                                                                                                                                                                                                          tttCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGGCGGGGTAGAAGGCATCTACCATGCTGGTGCTGCCCATCTTGCACACTCCGGTTATGCCCAAATTCCGCGGATCGACa    >  2:503252/1‑236 (MQ=255)
                                                                                                                                                                                                                            tCGCGCCCGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAAc  <  1:399359/236‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
CTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTATCGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAAC  >  NC_000913/3559539‑3559994

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