Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 526,982 A→T 6.5% N1241I (AAT→ATT)  rhsD → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913526,9820AT 6.5% 218.1 / 4.5 61N1241I (AAT→ATT) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (23/34);  new base T (2/2);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGAATACTCTGGTATTGTCTGTAAAGATACTAATGGAAAA  >  NC_000913/526696‑527229
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGACTGGGGCAACCAGCTTAATGAGGACAACCCGCATCATGTGTATCAGCCCTACCGTCTGCCCGGGCCGCCGCATGATGAGGAATCAGGGCTGTACTATAACGTCACCGGTACTACGATCCGTTCCAGGGTCGGTATATTACTCAGGACCCGATGGGTTTGCAAGGGGGATGGAAttt                                                                                                                                                                                                                                                      >  2:192269/1‑289 (MQ=255)
 cGGAAGACACGCCTGCGCGAAAAGCTCATCTTTATCAGTGCGAGCACCGGGGACTTCCGCGGGGGCTTATCAGCTAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGaagt                                                                                                                                                                                                                                                       <  2:618859/288‑3 (MQ=17)
         acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGt                                                                                                                                                                                                                                              <  1:627906/289‑1 (MQ=255)
          cacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTa                                                                                                                                                                                                                                             >  1:731799/1‑289 (MQ=255)
              ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCt                                                                                                                                                                                                                                         >  1:609467/1‑289 (MQ=255)
                   gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTaaa                                                                                                                                                                                                                                    <  1:710402/289‑1 (MQ=255)
                     aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATc                                                                                                                                                                                                                                  <  1:529901/289‑1 (MQ=255)
                             tCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTacaa                                                                                                                                                                                                                          <  1:587929/289‑1 (MQ=255)
                                     aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGa                                                                                                                                                                                                                  <  1:457909/289‑1 (MQ=255)
                                      ccgCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGAc                                                                                                                                                                                                                 <  1:258068/287‑1 (MQ=255)
                                       tGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGAc                                                                                                                                                                                                                 >  2:557078/1‑288 (MQ=255)
                                        gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGAc                                                                                                                                                                                                                 >  1:520991/1‑287 (MQ=255)
                                               cctggggCTGCAGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAACGGTACTACGAGCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGa                                                                                                                                                                                                        <  1:376964/286‑1 (MQ=255)
                                                      cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTg                                                                                                                                                                                                  <  2:71373/288‑1 (MQ=255)
                                                                      tCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGc                                                                                                                                                                                 <  1:397522/289‑1 (MQ=255)
                                                                                             tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACgggg                                                                                                                                                           <  2:609467/288‑1 (MQ=255)
                                                                                                               gaATGGGGCAGCCGGCTTAATGAGGGAAACCCGAACCATGTGTATCAGCCGTACAGTCTGCGAGGGCGGCAGAATGATCAGGAATCAGGGCGGTACAATAACCGCGACCGTGACGACGTTGCGTTGCAGGGGCGGTATATTACTCAGTACCCTATGGGGTTGAAAGGGTGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGGTTATTGCAGGCTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCgtt                                                                                                                                        <  2:67738/289‑1 (MQ=255)
                                                                                                                 aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCtt                                                                                                                                         >  1:745139/1‑286 (MQ=255)
                                                                                                                 aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCtt                                                                                                                                         <  2:745139/286‑1 (MQ=255)
                                                                                                                 aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTc                                                                                                                                       >  1:90543/1‑288 (MQ=255)
                                                                                                                    gggCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    <  2:90543/288‑1 (MQ=255)
                                                                                                                    gggCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    <  2:731799/288‑1 (MQ=255)
                                                                                                                    gggCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGGCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    <  2:531693/288‑1 (MQ=255)
                                                                                                                     ggCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGta                                                                                                                                  >  1:676475/1‑289 (MQ=255)
                                                                                                                                                          aTCAGGCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtg                                                                                              <  2:40743/288‑1 (MQ=255)
                                                                                                                                                              gCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    <  1:712729/246‑1 (MQ=255)
                                                                                                                                                              gCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    >  2:712729/1‑246 (MQ=255)
                                                                                                                                                               ccGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGATTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCt                                                                                                                                                                          <  1:185525/207‑1 (MQ=255)
                                                                                                                                                               ccGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGATTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCt                                                                                                                                                                          >  2:185525/1‑207 (MQ=255)
                                                                                                                                                                                     agcaTGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACg                                                                    <  2:520991/287‑1 (MQ=255)
                                                                                                                                                                                     agcaTGATGAGGAATCAGGGCTGTACTATAACGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACTCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCa                                                                  <  1:192269/288‑1 (MQ=255)
                                                                                                                                                                                         tgatgaGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAAt                                                              >  1:331186/1‑289 (MQ=255)
                                                                                                                                                                                          gatgaGGAAACGAGGCTGTACTATAACCGTCACCGGGACTACGATCGGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATa                                                             <  2:462745/289‑1 (MQ=255)
                                                                                                                                                                                                  aTCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTa                                                     <  2:331186/289‑1 (MQ=255)
                                                                                                                                                                                                        gCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                                                                                                    <  1:274298/220‑1 (MQ=255)
                                                                                                                                                                                                        gCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                                                                                                    >  2:274298/1‑220 (MQ=255)
                                                                                                                                                                                                        gCTGTACTATAAACGTCAACGGTACTAAGATCCGTTGCAGGGGCGGTATATTACTCAGGAACCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTaa                                                <  2:149524/288‑1 (MQ=255)
                                                                                                                                                                                                            tACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    >  1:463812/1‑200 (MQ=255)
                                                                                                                                                                                                            tACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    <  2:463812/200‑1 (MQ=255)
                                                                                                                                                                                                             aCTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAAtt                                                                                                                                                                                                                    <  2:582258/119‑1 (MQ=255)
                                                                                                                                                                                                             aCTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAAtt                                                                                                                                                                                                                    >  1:582258/1‑119 (MQ=255)
                                                                                                                                                                                                                ataACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGa                                       <  1:557078/289‑1 (MQ=255)
                                                                                                                                                                                                                 taACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCa                                                                  <  2:207521/261‑1 (MQ=255)
                                                                                                                                                                                                                 taACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCa                                                                  >  1:207521/1‑261 (MQ=255)
                                                                                                                                                                                                                 taACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGaa                                      <  2:676475/289‑1 (MQ=255)
                                                                                                                                                                                                                    ccGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTGTATGTAATGAAATGGAAt                                     >  1:442944/1‑287 (MQ=255)
                                                                                                                                                                                                                      gTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGATTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGAATAc                                   >  2:77271/1‑287 (MQ=255)
                                                                                                                                                                                                                      gTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    >  2:595822/1‑190 (MQ=255)
                                                                                                                                                                                                                      gTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACg                                                                                                                                    <  1:595822/190‑1 (MQ=255)
                                                                                                                                                                                                                         aCCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtg                                                                                              <  2:605839/225‑1 (MQ=255)
                                                                                                                                                                                                                         aCCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtg                                                                                              >  1:605839/1‑225 (MQ=255)
                                                                                                                                                                                                                         aCCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACg                                                                    <  2:40870/251‑1 (MQ=255)
                                                                                                                                                                                                                         aCCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACg                                                                    >  1:40870/1‑251 (MQ=255)
                                                                                                                                                                                                                                ctacGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGATTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGAATACTCTGGTATTGTc                       <  1:77271/289‑1 (MQ=255)
                                                                                                                                                                                                                                        cGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGAATACTCTGGTATTGTCTGTAAAGa               >  1:72199/1‑289 (MQ=255)
                                                                                                                                                                                                                                                     gTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGAATACTCTGGTATTGTCTGTAAAGATACTAATGGaaaa  >  2:273093/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                aGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGc                                                                                                                                                                                 >  1:64243/1‑103 (MQ=255)
                                                                                                                                                                                                                                                                aGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGc                                                                                                                                                                                 <  2:64243/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACg                                                                    <  2:7416/212‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACg                                                                    >  1:7416/1‑212 (MQ=255)
                                                                                                                                                                                                                                                                  gACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGaa                                      >  1:147861/1‑240 (MQ=255)
                                                                                                                                                                                                                                                                  gACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGaa                                      <  2:147861/240‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
CCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATATTACTCAGGACCCGATGGGGTTGAAAGGGGGATGGAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAACGCAGAATAGATCTCTATGTAATGATATGGAATACTCTGGTATTGTCTGTAAAGATACTAATGGAAAA  >  NC_000913/526696‑527229

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