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

TGGTGGTTTTATTGTCGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTCCCTCATCAATGT  >  NC_000913/1854566‑1855029
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
tggtggTTTTATTGTCAATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACTATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTa                                                                                                                                                                                                                                   <  1:832716/239‑1 (MQ=255)
               cGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTGTAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCAAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCCACCAAATTTGAGACGTTCACCTTTTTATCGTTGGCAACCAtt                                                                                                                                                                                                                    >  2:117423/1‑239 (MQ=255)
               cGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAACTTTGTGACTGTACCCTTTTTATCGTTGGCACCCCtt                                                                                                                                                                                                                    >  1:380879/1‑239 (MQ=255)
                       ggTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTgg                                                                                                                                                                                                                            <  1:41404/223‑1 (MQ=255)
                        gtggTGCCTTGGGTGAAATAAGCTGCTTCTCTTCAAAAAACGATGAAGCGCGATTAGAAAAATGGATATTCCCCGTCAACAACGCTGCATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTGCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAAt                                                                                                                                                                                                            <  2:243888/235‑1 (MQ=255)
                         tAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTAAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTAACTGAAAAATCCCACTAAAATTTGCGCATGTAAACATTTTATCTTTGGCAACCATTTGCGTAAtt                                                                                                                                                                                                           >  2:828607/1‑238 (MQ=255)
                            tGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAACAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTACAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCt                                                                                                                                                                                                       >  1:245941/1‑239 (MQ=255)
                            tGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCTTTGGCAACCATTTGCTTAATTg                                                                                                                                                                                                          >  1:680590/1‑236 (MQ=255)
                            tGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCAACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTCATTCCCt                                                                                                                                                                                                       >  1:69490/1‑239 (MQ=255)
                                ttGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCCTTGGCAACCATTTGCTTAATTGCCTGaa                                                                                                                                                                                                    >  2:505497/1‑238 (MQ=255)
                                 tGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACCTCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAACTTCTCAAAATCGCGGCGAATGGTTTCTTCGGTTACCTTAAAAATCCCCATCCAATTTTAGAATGGCACCTTTTTTTCGTTGGCAAACATTTGcccttctcttataacc                                                                                                                                                                                                  >  2:41404/1‑224 (MQ=255)
                                 tGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTTTCGTTGGCAACCATTTGCTTAATTGccgtaat                                                                                                                                                                                                   >  2:753985/1‑233 (MQ=255)
                                 ggcttGCAATAAGCCGCTTCTCTTCATAAAACGATGAAGCCCGCTTATAAAAATGGATATTCGCCGTCAATATCGCTGTATTTAAAACCGCACCACCAGAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATc                                                                                                                                                                                                  <  1:191870/235‑1 (MQ=255)
                                       cAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCt                                                                                                                                                                                            >  1:579462/1‑239 (MQ=255)
                                             gctgctTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCACCTTCTCAAGATCACGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGACACCTTTTTTACATTGACAACCAATTCCTTAATTGCCTGCATCCTGCCTTTTgcc                                                                                                                                                                                      >  2:191870/1‑239 (MQ=255)
                                                tgttttctTCATAAAATGATGAACTGCGCTTAAAAAAATGGATATTCTCCCGCAACACCGCTTTATTTAAAACAGCACCACTATAGGTTCGGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACTTGAAAAATCCCCCTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCGTGAATCCTGTCTTTTgccgcc                                                                                                                                                                                   <  2:874819/234‑1 (MQ=255)
                                                gctTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTgccgcc                                                                                                                                                                                   <  1:193847/239‑1 (MQ=255)
                                                 ctTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTTACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCa                                                                                                                                                                                  >  1:56378/1‑239 (MQ=255)
                                                  ttCTCTTCATAAAACGACGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTTAGACTGTCACCTTTTTATCTTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTATTTTTCCCCCAc                                                                                                                                                                                 >  2:315318/1‑239 (MQ=255)
                                                    ctctTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCCCGa                                                                                                                                                                               >  1:345548/1‑239 (MQ=255)
                                                       ttCATAAAACGTTGAAGCGCCTTTATAAAAATGGATTTTCTCCGTCAACATCGCTGTATTTAAAACAGCCCCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTa                                                                                                                                                                            <  2:319756/239‑1 (MQ=255)
                                                       ttCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCCCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATTGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTa                                                                                                                                                                            <  1:828607/239‑1 (MQ=255)
                                                               aCGTTGAAGCGCGCTTAAAAAAATGGATTTTTCCGGTCAACATCGTTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAAGGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTg                                                                                                                                                                    <  2:570072/239‑1 (MQ=255)
                                                                  tttaaGCGCGCTTATAAAAAGGGATATTCCCCGTCAACCTCGCTTTATTTAAAACAGCCCCCCCATAGGTTCTGTTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTAt                                                                                                                                                                 <  2:855218/236‑1 (MQ=255)
                                                                   tGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCGGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTAt                                                                                                                                                                 <  1:505497/238‑1 (MQ=255)
                                                                           gcTTATAAAAATGGATTTTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACa                                                                                                                                                        <  1:753985/239‑1 (MQ=255)
                                                                           gcTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACa                                                                                                                                                        >  1:77053/1‑239 (MQ=255)
                                                                           ccTTAAAAAAATGGATTTTCCCCGTAAAACCCCCTTTATTTAAAACAGCCCCCCCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACa                                                                                                                                                        <  2:680590/238‑1 (MQ=255)
                                                                                         atatTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGa                                                                                                                                          >  1:546875/1‑239 (MQ=255)
                                                                                          attttcccgtaaactaccgtttATTTAAAACAGCCCCACCATAGGTTCTGGTCAAAAA‑GCCTTCACCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGAt                                                                                                                                         <  2:579462/218‑1 (MQ=255)
                                                                                               tccttaaacatgcttttATTTAAAACAGCCCCCCCATAGGTTTTGGTCAAAAA‑GCCTTCATCTTCCAGCCTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCg                                                                                                                                    <  2:77053/223‑1 (MQ=255)
                                                                                                ccctCAACATCTCTGTATTTAAAACACCACCCCCAATGTTTCTGGTCAAAAA‑GCCTTCATCCTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAACCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGc                                                                                                                                   <  2:99536/236‑1 (MQ=255)
                                                                                                  tttttaagaaacaattttAAAACCTCCCCCCCATTGTTTTTGGTAAAAAA‑GCCTTCCTCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAATTCCCCCTCAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCAt                                                                                                                                 <  2:253252/224‑1 (MQ=255)
                                                                                                             tGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGttt                                                                                                                      >  2:685995/1‑239 (MQ=255)
                                                                                                                         cgCACCACCATAGGTTTGGGTCAAAAA‑GCCTTCATCTTCCAGCTTCCCAAGATCGGGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATa                                                                                                          <  1:117423/238‑1 (MQ=255)
                                                                                                                         aGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCCCATTAGCGCGATTATCGCATAACCGAGGTTTACTTTCAAAATa                                                                                                          >  1:873733/1‑239 (MQ=255)
                                                                                                                         aGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTAAATCACATTAACCCGATTAACCCAAAACCGAGGTTTACTTTCaaaaaa                                                                                                          >  2:423916/1‑237 (MQ=255)
                                                                                                                           cccccccaTAGGTTC‑GGTGAAAAAGGCTTTCATCTTCCAGATTTCAAAGATCGCGGGTAATGGTTTTTTGGGTTACCGAAAAAATCCCACTTAAATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGATTAATTTCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAAc                                                                                                        <  2:380879/234‑1 (MQ=255)
                                                                                                                           caccacCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATTACATTAgc                                                                                                                                              >  1:327150/1‑201 (MQ=255)
                                                                                                                           caccacCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGtt                                                                                                                       >  1:253252/1‑224 (MQ=255)
                                                                                                                                     ggTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCTCTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTaca                                                                                                                                                                          <  2:405109/163‑1 (MQ=255)
                                                                                                                                     ggTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTaca                                                                                                                                                                          >  1:405109/1‑163 (MQ=255)
                                                                                                                                       ttCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCCCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAAt                                                                                              <  1:88047/237‑1 (MQ=255)
                                                                                                                                        tCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCa                                                                                                                                                                                                                                           <  1:203132/95‑1 (MQ=255)
                                                                                                                                        tCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCa                                                                                                                                                                                                                                           >  2:203132/1‑95 (MQ=255)
                                                                                                                                        tCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCCCGATTATCGCATAACCGGTGTTTACTTTCAAAATAACCTGTTTGAAtaaa                                                                                           >  2:475852/1‑236 (MQ=255)
                                                                                                                                        tCTGGGCAAAAA‑GCCTTATTTCTCCAGCTTCTCAAGATCGCTGGAAATGGTTTCCTCGGTTCCCTGAAAAATCCCCCTCAAATTTGAGACTGTCACCTTTTTAGCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATcac                                                                                           <  2:69490/239‑1 (MQ=255)
                                                                                                                                          tGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCCCATTATCGCGATTATCCCATAACCGATGTTTACTTTCAAAAAAACCTGTTTGAATccaag                                                                                         >  2:411814/1‑235 (MQ=255)
                                                                                                                                          tGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCCCATTAGCGCGATTATCGCATAACCGCTGTTTACTTTCAAAAAAACCTGTTTGAATCAAAg                                                                                         >  2:719095/1‑239 (MQ=255)
                                                                                                                                          tGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCAt                                                                                                                                 >  2:860116/1‑199 (MQ=255)
                                                                                                                                          tGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAg                                                                                         >  2:676885/1‑239 (MQ=255)
                                                                                                                                          tGGTCAAAAA‑GCCTTCATCTTCCAGCTTCCCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCAt                                                                                                                                 <  1:860116/199‑1 (MQ=255)
                                                                                                                                          tGGTCAAAAA‑GCCTTCATCTTCCAGCT‑‑TCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTCCCGCCACGATTACACCCCCGTATCTTTTTACATCCCATTATCGCGATTATCGCATAAACGCTGTTTACTTTCAAAAGAACCTGTTTGAATCACAAAt                                                                                       >  2:173905/1‑239 (MQ=255)
                                                                                                                                           ggTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGCAACATCCCACTCACATTTGAGACTGTCCCCTTTTTATCGTTGGCAACCATTTGCTTATTTGCCTGAATCCTGTCTTTTGCCTCCACGATTACACCCCTGTTTCTTTTTCCCTCCCATTAACGCGCTTTTCGCATAACCCCTTTTTACTTTCAAAAAAAACTTTTTTAATccacga                                                                                        >  2:122935/1‑234 (MQ=255)
                                                                                                                                            gTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTAAATCACATTAGCGCGATTATCGCTTAACCGATGTTTACTTTCAAAAAAACACGTTTGatacaaacat                                                                                       >  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)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
TGGTGGTTTTATTGTCGATAAAGGGTAGTGCCTTGCGTGCAATAAGCTGCTTCTCTTCATAAAACGATGAAGCGCGCTTATAAAAATGGATATTCTCCGTCAACATCGCTGTATTTAAAACAGCACCACCATAGGTTCTGGTCAAAAA‑GCCTTCATCTTCCAGCTTCTCAAGATCGCGGCGAATGGTTTCTTCGGTTACCTGAAAAATCCCACTCAAATTTGAGACTGTCACCTTTTTATCGTTGGCAACCATTTGCTTAATTGCCTGAATCCTGTCTTTTGCCGCCACGATTACACCCCTGTATCTTTTTACATCACATTAGCGCGATTATCGCATAACCGATGTTTACTTTCAAAATAACCTGTTTGAATCACAGATTTTCATCACAGTTTTCACAGAAACAGAGGTGAATCGTGTTGAGTATTTAACGCTCCAGGGCCTCTGCCATTTCCCTCATCAATGT  >  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: ‑