Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,854,796 A→C 19.0% V20G (GTG→GGG)  ydjF ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,854,7960AC19.0% 154.9 / 2.4 77V20G (GTG→GGG) ydjFputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/24);  new base C (3/14);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 4.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGGTTTTATTGTCGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTCCCTCATCAATGT  >  NC_000913/1854566‑1855029
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
tggtggTTTTATTGTCAATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACTATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTa                                                                                                                                                                                                                                   <  1:832716/239‑1 (MQ=255)
               cGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTGTAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCAAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCCACCAAATTTGAGACGTTCACCTTTTTATCGTTGGCAACCAtt                                                                                                                                                                                                                    >  2:117423/1‑239 (MQ=255)
               cGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAACTTTGTGACTGTACCCTTTTTATCGTTGGCACCCCtt                                                                                                                                                                                                                    >  1:380879/1‑239 (MQ=255)
                       ggTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTgg                                                                                                                                                                                                                            <  1:41404/223‑1 (MQ=255)
                        gtggTGCCTTGGGTGAAATAAGCTGCTTCTCTTCAAAAAACGATGAAGCGCGATTAGAAAAATGGATATTCCCCGTCAACAACGCTGCATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTGCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAAt                                                                                                                                                                                                            <  2:243888/235‑1 (MQ=255)
                         tAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTAAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTAACTGAAAAATCCCACTAAAATTTGCGCATGTAAACATTTTATCTTTGGCAACCATTTGCGTAAtt                                                                                                                                                                                                           >  2:828607/1‑238 (MQ=255)
                            tGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAACAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTACAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCt                                                                                                                                                                                                       >  1:245941/1‑239 (MQ=255)
                            tGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCTTTGGCAACCATTTGCTTAATTg                                                                                                                                                                                                          >  1:680590/1‑236 (MQ=255)
                            tGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCAACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTCATTCCCt                                                                                                                                                                                                       >  1:69490/1‑239 (MQ=255)
                                ttGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCCTTGGCAACCATTTGCTTAATTGCCTGaa                                                                                                                                                                                                    >  2:505497/1‑238 (MQ=255)
                                 tGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACCTCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAACTTCTCAAAATCGCGGCGAATGGTTTCTTCGGTTACCTTAAAAATCCCCATCCAATTTTAGAATGGCACCTTTTTTTCGTTGGCAAACATTTGcccttctcttataacc                                                                                                                                                                                                  >  2:41404/1‑224 (MQ=255)
                                 tGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTTTCGTTGGCAACCATTTGCTTAATTGccgtaat                                                                                                                                                                                                   >  2:753985/1‑233 (MQ=255)
                                 ggcttGCAATAAGCCGCTTCTCTTCATAAAACGATGAAGCCCGCTTATAAAAATGGATATTCGCCGTCAATATCGCTGTATTTAAAACCGCACCACCAGAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATc                                                                                                                                                                                                  <  1:191870/235‑1 (MQ=255)
                                       cAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCt                                                                                                                                                                                            >  1:579462/1‑239 (MQ=255)
                                             gctgctTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCACCTTCTCAAGATCACGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGACACCTTTTTTACATTGACAACCAATTCCTTAATTGCCTGCATCCTGCCTTTTgcc                                                                                                                                                                                      >  2:191870/1‑239 (MQ=255)
                                                tgttttctTCATAAAATGATGAACTGCGCTTAAAAAAATGGATATTCTCCCGCAACACCGCTTTATTTAAAACAGCACCACTATAGGTTCGGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACTTGAAAAATCCCCCTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCGTGAATCCTGTCTTTTgccgcc                                                                                                                                                                                   <  2:874819/234‑1 (MQ=255)
                                                gctTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTgccgcc                                                                                                                                                                                   <  1:193847/239‑1 (MQ=255)
                                                 ctTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTTACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCa                                                                                                                                                                                  >  1:56378/1‑239 (MQ=255)
                                                  ttCTCTTCATAAAACGACGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTTAGACTGTCACCTTTTTATCTTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTATTTTTCCCCCAc                                                                                                                                                                                 >  2:315318/1‑239 (MQ=255)
                                                    ctctTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCCCGa                                                                                                                                                                               >  1:345548/1‑239 (MQ=255)
                                                       ttCATAAAACGTTGAAGCGCCTTTATAAAAATGGATTTTCTCCGTCAACATCGCTGTATTTAAAACAGCCCCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTa                                                                                                                                                                            <  2:319756/239‑1 (MQ=255)
                                                       ttCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCCCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATTGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTa                                                                                                                                                                            <  1:828607/239‑1 (MQ=255)
                                                               aCGTTGAAGCGCGCTTAAAAAAATGGATTTTTCCGGTCAACATCGTTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAAGGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTg                                                                                                                                                                    <  2:570072/239‑1 (MQ=255)
                                                                  tttaaGCGCGCTTATAAAAAGGGATATTCCCCGTCAACCTCGCTTTATTTAAAACAGCCCCCCCATAGGTTCTGTTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTAt                                                                                                                                                                 <  2:855218/236‑1 (MQ=255)
                                                                   tGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCGGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTAt                                                                                                                                                                 <  1:505497/238‑1 (MQ=255)
                                                                           gcTTATAAAAATGGATTTTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACa                                                                                                                                                        <  1:753985/239‑1 (MQ=255)
                                                                           gcTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACa                                                                                                                                                        >  1:77053/1‑239 (MQ=255)
                                                                           ccTTAAAAAAATGGATTTTCCCCGTAAAACCCCCTTTATTTAAAACAGCCCCCCCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACa                                                                                                                                                        <  2:680590/238‑1 (MQ=255)
                                                                                         atatTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGa                                                                                                                                          >  1:546875/1‑239 (MQ=255)
                                                                                          attttcccgtaaactaccgtttATTTAAAACAGCCCCACCATAGGTTCTGGTCAAAAAGCCTTCACCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGAt                                                                                                                                         <  2:579462/218‑1 (MQ=255)
                                                                                               tccttaaacatgcttttATTTAAAACAGCCCCCCCATAGGTTTTGGTCAAAAAGCCTTCATCTTCCAGCCTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCg                                                                                                                                    <  2:77053/223‑1 (MQ=255)
                                                                                                ccctCAACATCTCTGTATTTAAAACACCACCCCCAATGTTTCTGGTCAAAAAGCCTTCATCCTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAACCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGc                                                                                                                                   <  2:99536/236‑1 (MQ=255)
                                                                                                  tttttaagaaacaattttAAAACCTCCCCCCCATTGTTTTTGGTAAAAAAGCCTTCCTCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAATTCCCCCTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCAt                                                                                                                                 <  2:253252/224‑1 (MQ=255)
                                                                                                             tGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGttt                                                                                                                      >  2:685995/1‑239 (MQ=255)
                                                                                                                         cgCACCACCATAGGTTTGGGTCAAAAAGCCTTCATCTTCCAGCTTCCCAAGATCGGGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATa                                                                                                          <  1:117423/238‑1 (MQ=255)
                                                                                                                         aGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCCCATTAGCGCGATTATCGCATAACCGAGGTTTACTTTCAAAATa                                                                                                          >  1:873733/1‑239 (MQ=255)
                                                                                                                         aGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTAAATCACATTAACCCGATTAACCCAAAACCGAGGTTTACTTTCaaaaaa                                                                                                          >  2:423916/1‑237 (MQ=255)
                                                                                                                           cccccccaTAGGTTCGGTGAAAAAGGCTTTCATCTTCCAGATTTCAAAGATCGCGGGTAATGGTTTTTTGGGTTACCGAAAAAATCCCACTTAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGATTAATTTCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAAc                                                                                                        <  2:380879/234‑1 (MQ=255)
                                                                                                                           caccacCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATTACATTAgc                                                                                                                                              >  1:327150/1‑201 (MQ=255)
                                                                                                                           caccacCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGtt                                                                                                                       >  1:253252/1‑224 (MQ=255)
                                                                                                                                     ggTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCTCTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTaca                                                                                                                                                                          <  2:405109/163‑1 (MQ=255)
                                                                                                                                     ggTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTaca                                                                                                                                                                          >  1:405109/1‑163 (MQ=255)
                                                                                                                                       ttCTGGTCAAAAAGCCTTCATCTTCCAGCTTCCCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAAt                                                                                              <  1:88047/237‑1 (MQ=255)
                                                                                                                                        tCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCa                                                                                                                                                                                                                                           <  1:203132/95‑1 (MQ=255)
                                                                                                                                        tCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCa                                                                                                                                                                                                                                           >  2:203132/1‑95 (MQ=255)
                                                                                                                                        tCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCCCGATTATCGCATAACCGGTGTTTACTTTCAAAATAACCTGTTTGAAtaaa                                                                                           >  2:475852/1‑236 (MQ=255)
                                                                                                                                        tCTGGGCAAAAAGCCTTATTTCTCCAGCTTCTCAAGATCGCTGGAAATGGTTTCCTCGGTTCCCTGAAAAATCCCCCTCAAATTTGAGACTGTCACCTTTTTAGCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATcac                                                                                           <  2:69490/239‑1 (MQ=255)
                                                                                                                                          tGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCCCATTATCGCGATTATCCCATAACCGATGTTTACTTTCAAAAAAACCTGTTTGAATccaag                                                                                         >  2:411814/1‑235 (MQ=255)
                                                                                                                                          tGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCCCATTAGCGCGATTATCGCATAACCGCTGTTTACTTTCAAAAAAACCTGTTTGAATCAAAg                                                                                         >  2:719095/1‑239 (MQ=255)
                                                                                                                                          tGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCAt                                                                                                                                 >  2:860116/1‑199 (MQ=255)
                                                                                                                                          tGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAg                                                                                         >  2:676885/1‑239 (MQ=255)
                                                                                                                                          tGGTCAAAAAGCCTTCATCTTCCAGCTTCCCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCAt                                                                                                                                 <  1:860116/199‑1 (MQ=255)
                                                                                                                                          tGGTCAAAAAGCCTTCATCTTCCAGCT‑‑TCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTCCCGCCACGATTACACCCCCGTATCTTTTTACATCCCATTATCGCGATTATCGCATAAACGCTGTTTACTTTCAAAAGAACCTGTTTGAATCACAAAt                                                                                       >  2:173905/1‑239 (MQ=255)
                                                                                                                                           ggTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGCAACATCCCACTCACATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTATTTGCCTGAATCCTGTCTTTTGCCTCCACGATTACACCCCTGTTTCTTTTTCCCTCCCATTAACGCGCTTTTCGCATAACCCCTTTTTACTTTCAAAAAAAACTTTTTTAATccacga                                                                                        >  2:122935/1‑234 (MQ=255)
                                                                                                                                            gTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTAAATCACATTAGCGCGATTATCGCTTAACCGATGTTTACTTTCAAAAAAACACGTTTGatacaaacat                                                                                       >  2:579602/1‑230 (MQ=255)
                                                                                                                                                       cTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATcaca                                                                             >  1:669927/1‑239 (MQ=255)
                                                                                                                                                       cTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATcaca                                                                             >  1:616317/1‑239 (MQ=255)
                                                                                                                                                        ttCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCTACGATTACACCCCTGTATCTTTTTACATCCCATTAGCCCGATTCACGCATAACCGATGTTTACTTTCAAAATAACCCGTTTGAATcacaacattttaacacag                                                                            >  2:854235/1‑225 (MQ=255)
                                                                                                                                                               ccaagtttCTCAAGATCGCGGCGAATGGTTTCTTCTGTTACCTGAAAATTCCCACTCAAATGTGAGACTGTCGCCTTTTGATGGTTGGCAACCATTTGGTTAATGGCCTGAATCCTGTTTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATAACATTAGCGCGATTATCGCAGAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTcac                                                                     <  1:251126/233‑1 (MQ=255)
                                                                                                                                                                     ttCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGGCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGaaa                                                                <  2:108557/238‑1 (MQ=255)
                                                                                                                                                                       ctcAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAAt                                                                                              >  1:308130/1‑206 (MQ=255)
                                                                                                                                                                               cgcgtgtattgttttttttgtttACTTGAAAAACCCCCCCAAAAATGGAACTTGTCCCCTTTTGTTGGTGGAAAACTATTTGTTTATTTGCCTGAACCCTGTTTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGGTTGAATCACAGATTTTCATTACAGTTTTCACAGAAACAGAGGTGAAt                                                     <  2:669927/218‑1 (MQ=255)
                                                                                                                                                                                  ggcgAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCCCAAAAACAGAGGTaaacagt                                                  >  2:194690/1‑235 (MQ=255)
                                                                                                                                                                                  ggcgAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAAAAGAGGAGAATCgg                                                  >  2:716840/1‑238 (MQ=255)
                                                                                                                                                                                     gattgtttttttttgTTTCCTCAAAAACCCCCCTCAAATTTGAGACTTCCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTGGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGtt                                               <  2:546875/225‑1 (MQ=255)
                                                                                                                                                                                      gatttttttttTGGTTCCCTGAAAAATCCCCCTCAATTTGGGGACTTGCCCCTTTTTAACGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTGGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGtt                                               <  2:873733/233‑1 (MQ=255)
                                                                                                                                                                                      aaTGGTTTCTTCGGTTTCCTGAAAAATCCCCCCCAAATTTGGAACTTTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGtt                                               <  1:423916/238‑1 (MQ=255)
                                                                                                                                                                                      aaTGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGtt                                               <  1:676885/238‑1 (MQ=255)
                                                                                                                                                                                      aaTGGTTTCTTCGGTTAACGTAAAAATCCCACTCAAATTTGAGTCTTTCACCTTTTTTCCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGGCTTTTGCCGCCACGATTACGCCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGtt                                               <  1:845189/238‑1 (MQ=255)
                                                                                                                                                                                       aTGGTTTCTTCGGTTACCTGAAAAATCCCACGCAAATTTGAGCCTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGa                                             <  1:411814/239‑1 (MQ=255)
                                                                                                                                                                                               tttgttttccTGAAAAATCCCACCCAAATTTGGGACTGTCCCCTTTTTGCCGGTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAACAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTaa                                     <  2:345548/231‑1 (MQ=255)
                                                                                                                                                                                               ttcGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTaa                                     <  1:122935/239‑1 (MQ=255)
                                                                                                                                                                                               ttcGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTaa                                     <  1:719095/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTTCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTATTTTCTAAATAACCTGTTTGAATCACCGATTTTCATCACATTTTTCACGTAAACAGAGGTGACACGAGTTGTGTTTTTAAc                                    >  2:705794/1‑238 (MQ=255)
                                                                                                                                                                                                             aaaTCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCctc                       <  1:315318/239‑1 (MQ=255)
                                                                                                                                                                                                                       cAAATTTGGGTCTTTCCCCTTTTTTGCGTTTGCAACCATTTGCTTAATTCCCTTAAACCTGTCTTTTGCAGCCACGATTACACCCCTTTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTcc             <  2:741724/239‑1 (MQ=255)
                                                                                                                                                                                                                            ttGGGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTCCCtcat        <  1:579602/239‑1 (MQ=255)
                                                                                                                                                                                                                             tgggctttgcACCTTTTTACCGTTGGCACCCTTTGTCCTAATTGCCTGAATCCTTGTTTTTGCCGCCACGATTACCCCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTCCCtcatc       <  1:303848/230‑1 (MQ=255)
                                                                                                                                                                                                                                  atGTCAACTTTTTATCGCTGCCAACCATTTGATTAATTGACTGGATCCTGTCTTTTGCTGCCACGATTACACCCCTGTATCTTTTTACATCACCTTTGCGCGATTATCCCATTATCGATGGTTACTTTCATAATAACCTGTTTGAATCACAGCTTTTACTCACAGTTTTCACCGAAAAAGAGGTGGATCGCCTTGCGACTTTAAAGATGCAGGGCCTCTACCATATCCCTCAACAAttt  >  2:34708/2‑237 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
TGGTGGTTTTATTGTCGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAAGCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTCCCTCATCAATGT  >  NC_000913/1854566‑1855029

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