Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,434 C→T N16N (AAC→AAT)  araC → DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,4340CT97.6% 118.2 / ‑4.1 42N16N (AAC→AAT) araCDNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/20);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

GGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGT  >  NC_000913/70187‑70673
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gttctatgcTGCCGCGATGATAGTCTCTTTTGTTACGCTCTTTTGTCATGGCTTTGGTCCCGCGTTGTTACTGAATGCTTTAAATAATCGGGGTCACCGGTTTGGTGAGCGAGAAGAGCCAGGAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAGTGGTTTCTTCTCTGAATTGCGGGAGTATGAAAAGTATGTCTGAAGCGCAAAATGATCCCCTGCTGTCGGGATACTCGATTAATGcc                                                                                                                                                                                                                                                  <  2:516151/243‑1 (MQ=255)
                       actcTTTTGTTACGCGGTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGGTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                                                                                                                                                                                                           <  1:434164/248‑1 (MQ=255)
                              tGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTACTCGTTTAATGCCCATCTGGGGGCGGGGTTAAACCCGATTggg                                                                                                                                                                                                                    >  1:46907/1‑249 (MQ=255)
                                ttACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGAATGGCTGGAGCGCCAAATGATTCCCTGCTGGCGGGGTACTCGTTTTATGCCCCTCTGGTGGCGGGttttgagccgcttgaggg                                                                                                                                                                                                                  >  2:98028/1‑236 (MQ=255)
                                   cgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGTTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGccc                                                                                                                                                                                                                >  1:199882/1‑249 (MQ=255)
                                                       ggTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGATTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGCTTGTGGCCAACGGGTATTTCGAtttttttt                                                                                                                                                                                           >  1:624616/1‑250 (MQ=255)
                                                           acGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCATGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGACGATTGGGGCCAACGGTTTTCTCGATTTTTTTATcg                                                                                                                                                                                        >  1:24829/2‑250 (MQ=255)
                                                             gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc                                                                                                                                                                                     <  2:521364/251‑1 (MQ=255)
                                                             gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGGCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc                                                                                                                                                                                                                              <  1:407069/210‑1 (MQ=255)
                                                              cTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTTACGCCGCTTGAGGCCCACGGTTTTTTCGATTTTTTTTTcgaccc                                                                                                                                                                                    >  1:389190/1‑250 (MQ=255)
                                                               tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGGAAAGTATGGGTGAAGAGCCAAATGATCCCCTGCTGCCGGGATAATCGTTTAATGCCCATCTGGTGGCGGGTTTTACGCCCATTTAGGCCCAAGGGTATCTCGATTTTTTTTTcgcacgg                                                                                                                                                                                   >  1:697393/1‑246 (MQ=255)
                                                                ttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGCCCACGGGTATCTCGGTTTTTTTATcgaccgac                                                                                                                                                                                  >  1:604919/1‑251 (MQ=255)
                                                                      cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTgg                                                                                                                                                                            <  1:267106/251‑1 (MQ=255)
                                                                          aTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAAt                                                                                                                                                                        <  1:98028/251‑1 (MQ=255)
                                                                                   ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTTAAGGCGATTGGGGGCCACGGGTTTCTCGGTTTTTTTTTCGACCGACCGCTgggggatcggaaa                                                                                                                                                                 >  1:486778/1‑239 (MQ=255)
                                                                                   aaaaGCGTGGTTACCTGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTggg                                                                                                                                                                           <  2:486778/237‑1 (MQ=255)
                                                                                         ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                                                                                                                                                                                                                                  >  1:227290/1‑178 (MQ=255)
                                                                                         ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                                                                                                                                                                                                                                  <  2:227290/178‑1 (MQ=255)
                                                                                            gttttcGGTTTGGTTAGCGGGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCACTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAAGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑aa                                                                                                                                                   <  2:429565/246‑1 (MQ=255)
                                                                                               actggTTTGGTTAGAGAGAAGAGCCATTAAAAGACGCAGTTACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAAGGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGTCAACGGTTATCTCGATTTTTT‑ATCGACCGACCGCTGTGAATGAAA‑GGTTATATTCT‑C‑AAtctc                                                                                                                                               <  2:81429/248‑1 (MQ=255)
                                                                                               accgtttttgTTAGCGAGAATAGCCAGTAAAAGCCGCAGGGACGGTAATGTCTGATGCAATATGGACAATTGGTTTCGTCTCTGAATGGCGGGAGGATGAAAAGTATTGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTGTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AAtct                                                                                                                                                <  2:383132/242‑1 (MQ=255)
                                                                                               aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AAtct                                                                                                                                                <  2:199882/251‑1 (MQ=255)
                                                                                                 cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGGGGCCAACGGTTAt                                                                                                                                                                                                        >  1:189118/1‑196 (MQ=255)
                                                                                                        gTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTAATCGTTTAATGCCCATCTGGGGGCGGGGTTTACGCCGGTTTGGGCCAACGGTTCTCTCGATTTTTTTATTGACCGCACGCGGGGGATGAAA‑GGGTATTTTTT‑C‑AATCTCACCATTcg                                                                                                                                       >  2:192308/1‑251 (MQ=255)
                                                                                                                    agCCAGTTAAAGTCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGggt                                                                                                                           <  2:397649/251‑1 (MQ=255)
                                                                                                                           agaaGACGCAGTGAAGGTAATGGCTGATGTAATAGGGACGATGTGTTTCTTCTCTTAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAATATGATCCCCTGCTGCCTGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGAGTGAGGCCAACGGTTATCGCGATTTTTTTATCGAGCGACTGCTGGGGAGGAAA‑GGTTATATTCG‑C‑AATCTCCCCATTCGCGGTCAGGGGGTGG‑TGaaa                                                                                                                   <  2:51452/249‑1 (MQ=255)
                                                                                                                           aaaaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGTCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGaaa                                                                                                                   <  1:321355/251‑1 (MQ=255)
                                                                                                                                    aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATg                                                                                                                                                                       >  1:373188/1‑194 (MQ=255)
                                                                                                                                    aGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATg                                                                                                                                                                       <  2:373188/194‑1 (MQ=255)
                                                                                                                                               tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGGATGAAAAGTATGGCTGAAGCGCAAAAAGATCCCCTGCTGCCGGGATAATCGTTTAATGCACATTTGGTGGAGGGGTTAAAGCCGATTTAGGACCACGGGTCTCTCGAATTTTTTTTCGGCCGTCAGGTGGGAATGAAAAGGTT‑TATTCT‑CCAAT‑TAACCAATACAGGTCAGGGGGTGGGTGATAATTCAGG‑ACGAGAAttttt                                                                                               >  1:86873/1‑249 (MQ=255)
                                                                                                                                                 tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGAATGGAAAGTATGGCTGCAGCGCAAAATCATCCCCTGCTGCCGGGATACTCGATTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGAGGCCACAGGTTTTCTCGTTTTTTTTATTGCACGCCCGCTGGGACTTAAA‑CGTTATTTTCT‑C‑CATCTTCCCATTCTCGCGCAGGGGGGGG‑GGAAAAATCAGGGatcaggatttgttt                                                                                             >  1:40661/1‑238 (MQ=255)
                                                                                                                                                    gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGttt                                                                                             <  1:192308/248‑3 (MQ=255)
                                                                                                                                                    gATGCAATATGGACAATTGGTTTCTTCTATGAATGGCGGGAGTATGAAAAGTATGGCTGCAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACGGTTATCTCCGTTTTTTTATCGACCGACCCCTGGGGATGAAA‑GGTTTTATTCT‑C‑AATCTCACCCTTCGCGGTCAGGGGGGGG‑GGAAAAATTAGGGACGAGAATTTGTTTGcc                                                                                          >  1:136549/1‑251 (MQ=255)
                                                                                                                                                                aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                                                                                                                                                                                                    <  2:593818/105‑1 (MQ=255)
                                                                                                                                                                aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                                                                                                                                                                                                    >  1:593818/1‑105 (MQ=255)
                                                                                                                                                                        tttCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCCAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAAGCCGATTGAAGCCAAAGGTTATCTCGGTTTTTTTTTCGAACGACCGCTGGGAATGAAA‑GATTATATTCTTC‑AAT‑TTACCAATTGCGGGCATGGGGTGG‑TGAAAAATCCGGGGGGAGAATTTTTTTTCCGACCGGGGGATATTTTGCTg                                                                      >  1:191411/1‑251 (MQ=255)
                                                                                                                                                                                ctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc                                                                 <  2:604919/248‑1 (MQ=255)
                                                                                                                                                                                  gAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                                              <  1:11023/249‑1 (MQ=255)
                                                                                                                                                                                                                             ccccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTGAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATg                 <  1:163747/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGGTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCCGGGGGTGG‑TGAAAAATCAGGGCCGCGAATTTTTTTGCCCACCGGGTGATATTTTACTGTTCCCGCCAGGAGAGTTTACTAACTACGGTCGTCCTCCGGAGGCTCGCGAATGGTAtacccagggggt  >  2:613841/1‑240 (MQ=255)
                                                                                                                                                                                                                                             aCTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑a                                                                                                                                                    <  1:264464/105‑1 (MQ=255)
                                                                                                                                                                                                                                             aCTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑a                                                                                                                                                    >  2:264464/1‑105 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAA‑GGTTATATTCT‑C‑AATCTCACCATTCGCGGTCAGGGGGTGG‑TGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGT  >  NC_000913/70187‑70673

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