Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,411,689 G→T E13* (GAA→TAA)  yhdJ → DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,411,6890GT93.9% 107.4 / 2.5 51E13* (GAA→TAA) yhdJDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (23/23);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

TAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTAT  >  NC_000913/3411443‑3411933
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGCCCCGGTTTGGTAATTAAgg                                                                                                                                                                                                                                                  >  2:336912/1‑250 (MQ=255)
          ccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACACGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTAGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCatt                                                                                                                                                                                                                                        >  2:781042/1‑251 (MQ=255)
                ggACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGAACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTACGACCATTATTCAc                                                                                                                                                                                                                                  >  1:530871/1‑251 (MQ=255)
                ggACATGGTGATGCAATACACCCGTGGTAACAAGAGCCGTACTGCGCTGGAGGTGGGCATCAACCGTGGTACGCTGCGTAAAACATTGAAAAAATACGGAATGAACTAATTCAGATTATCGACATGATTGAATAATAAGGCGCTCCTCGGTGTGGTGAAGAGCCTTATTTATAGCTGTCACAAAGGGGGTGACAATTGGATAAGGATGAGCGCCGAGCCGGTTTGGTAATTTAGCTAAGAACATAATTCAc                                                                                                                                                                                                                                  >  2:679735/1‑251 (MQ=255)
                            gccaaccccccGTTGTAACCATACCCGTTCTGCGCTTATGATGGGCAGCAACCGTGGTACGTGGCGTAAATAATTGAAAAAAGACGGCATGAACTAATTCATGTTAGCTAAATGCTTGATTAAAAAGGCGCTATTCGGCAGGGGGAAGCGCCTTTTGTATAGGCGTCACAAAGGGAGTGACCATGAGACCAGGAGGTGAACCGACCCGGTTTGGTAATTAAGCTCAGACCATTATTCACGGTGATGCGCtt                                                                                                                                                                                                                      <  2:564071/243‑1 (MQ=255)
                                tACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGGGATGCGCTTGCCg                                                                                                                                                                                                                  >  2:517699/1‑251 (MQ=255)
                                   acCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAAc                                                                                                                                                                                                               >  1:22667/1‑251 (MQ=255)
                                             aCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCTCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaaa                                                                                                                                                                                                         >  1:378264/1‑247 (MQ=255)
                                              ccAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAAATGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATc                                                                                                                                                                                                    <  1:181537/251‑1 (MQ=255)
                                                aGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTTACGGTGATGCGCTTGCCGAACTTAAAAAGATccc                                                                                                                                                                                                  >  2:1991/1‑251 (MQ=255)
                                                 gACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATcccc                                                                                                                                                                                                 <  1:262333/251‑1 (MQ=255)
                                                  aCCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCccg                                                                                                                                                                                                <  2:571493/251‑1 (MQ=255)
                                                  aCCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCccg                                                                                                                                                                                                >  2:79356/1‑251 (MQ=255)
                                                   cccGTGATGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACACAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCccgc                                                                                                                                                                                               >  1:102796/1‑251 (MQ=255)
                                                                 atgatgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTg                                                                                                                                                                                                                     >  1:424233/1‑215 (MQ=255)
                                                                 atgatgGGCATCAACCGTGGTACGCTGCGTAAAAAAGTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTg                                                                                                                                                                                                                     <  2:424233/215‑1 (MQ=255)
                                                                      gggaaccAACCGTGGTTCGTGGCTTAAAAAATGGAAAAAATACGGAATGAACAAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTCCTCGGCATGGGGAAGCGCCTTTTTTTTAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGGGAACCGCCCCGGTTTGGAAATTAAGCTAAGTCCATTTTTCACGGTGATGCGCTTGCCGAACTTAAAAAGAGCCCCGCCGAAAGTGTCgatctgatc                                                                                                                                                                            <  2:763782/245‑1 (MQ=255)
                                                                        gCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGGCGATCTGATCtt                                                                                                                                                                          <  1:662009/251‑1 (MQ=255)
                                                                         cATCAACCGGGGTCCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTAAGCTAAATGTTGGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTTAAAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATGCACGGTGATGCGCTTGCCTAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCttt                                                                                                                                                                         <  2:50594/251‑1 (MQ=255)
                                                                          aTCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGACGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGTGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGTCCATTAATCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTg                                                                                                                                                                        >  1:445153/1‑251 (MQ=255)
                                                                              aCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCg                                                                                                                                                                     >  1:322490/1‑250 (MQ=255)
                                                                                      aCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTATCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGATTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACccacca                                                                                                                                                             <  1:781042/250‑1 (MQ=255)
                                                                                      aCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGAGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCat                                                                                                                                                            >  1:566003/1‑251 (MQ=255)
                                                                                                  aaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTATCGATCTGATCTTTGCCGACCCACCATATAACATCGGt                                                                                                                                                 >  1:809784/1‑250 (MQ=255)
                                                                                                   aaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTa                                                                                                                                                >  1:365628/1‑250 (MQ=255)
                                                                                                        aaaaaaTACGGCAGGAACGAATTCAGGTTATCTAAATGCTTGATTAAAAAGTCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGAGGCGCTTGCCGACCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAAtt                                                                                                                                          <  1:577321/251‑1 (MQ=255)
                                                                                                                       aactaatgaaggtaagtgaaaTGCTGGAGTAAAAAGGCGCGACTCGGAATGGGGAAGCGTCTTTGTTATAGGGGGCACAAAGGGAGGAACCATGAGAAGAGGATGTGAACCGACCCGGTTTGGTAATAAAGTTAAGACCATTATTCGGGGTGATGCGCTTGCCGATCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCAGATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                                           <  2:101156/233‑1 (MQ=255)
                                                                                                                               aaggtaagatAAATGTTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCAGGAGACCAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGCCCATTATTCACGTGGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTg                                                                                                                    <  1:459258/241‑1 (MQ=255)
                                                                                                                                  ggtatcTAAATGCTTGATTAAAAAGGCTCTCCTGGGCATGGGGAAGCGCTTTTTTTATAGGTGCCACAAAGGGATTGACCATGAGAACAGGATGTGAACCGCCCCGGTTTGGTAATTAAGTGAAGACCATTATTCACGGTGATGCTCTTGCCGAACTTAAAAAGATCCCCGCCTAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaaa                                                                                                                <  2:378264/246‑1 (MQ=255)
                                                                                                                                         aaaTGTTTGATAAAAAAGGCGGTACGCGGCAGGGGGACGCGCCTTTTTTAGAGGTGTCAAAAAGGGAGTGACCATGAGAACAGGAGGGGAACCGCCCCGGTTTGGTAATTAAGCTAAGCCCATTATTCACGTGGAGGCGCTTGCAGAATTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGTTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATc                                                                                                         <  2:530756/251‑1 (MQ=255)
                                                                                                                                           gtgcttgctTAAAAAGGCGCTACTCGGCCTGGGGACCCGCCTTTCTTCTAGGAGTTACAAAGGGAGTGACCATGAGAACAGGCTGAGAACCGCCCCGGTTTGGTAATGCACCTAAGAACCTTATTCACGGTGATGCGCTTGCCGAAATTAAAAAGATCCCCGCCGAAAGTGTCGATCTGTTCTTTGCCGACCCACACTATAACACCGGTAACAATTCTGATTGTCTGATAGAAGCCTGGAATGTAGATCTg                                                                                                       >  2:412704/2‑251 (MQ=255)
                                                                                                                                                 gATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTGGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATc                                                                                                 <  1:517699/251‑1 (MQ=255)
                                                                                                                                                  aTTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCg                                                                                                >  2:676878/1‑251 (MQ=255)
                                                                                                                                                     aaaaaGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGAc                                                                                              <  1:336912/250‑1 (MQ=255)
                                                                                                                                                            cgcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGtt                                                                                      >  1:446931/1‑251 (MQ=255)
                                                                                                                                                                   cGGCAGGGGGAAGCGCGTTGTTTATAGGGCTCACAATGGGAGTTACCATGAGAACAGGTTGTAAACCTCCCCGGTTCGTTAATTAAGCAAAGCCCATTATCCACGGGTATGCTCTTGCGTAACTAAAAAAGTTCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCGGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGaa                                                                                  <  2:610200/248‑1 (MQ=255)
                                                                                                                                                                    ggCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTg                                                                               >  1:6911/1‑250 (MQ=255)
                                                                                                                                                                      cATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGAtt                                                                            <  1:82194/251‑1 (MQ=255)
                                                                                                                                                                       aTGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTg                                                                           >  1:226154/1‑251 (MQ=255)
                                                                                                                                                                            gAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCt                                                                                                                                <  1:584159/193‑1 (MQ=255)
                                                                                                                                                                            gAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCt                                                                                                                                >  2:584159/1‑193 (MQ=255)
                                                                                                                                                                                cgcCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGAAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGcc                                                                  <  1:5457/251‑1 (MQ=255)
                                                                                                                                                                                           taGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTg                                                       >  1:291775/1‑251 (MQ=255)
                                                                                                                                                                                           aagggttctcaAATGGAGAGACCAGGAGAACAGGAGGGAACCCGACCGGGTTTGGTAATTAATATAAGCCCTGTATTCACGTTGATGCTCTTGCCGAACTAAAAAATATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAATAAGATCTGTTTATCGACGGGCTGTTTGAAGTGATGGCAGAGTGCCATCGCGTTCTg                                                       <  2:125676/242‑1 (MQ=255)
                                                                                                                                                                                              gtgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGaa                                                     >  1:517778/1‑250 (MQ=255)
                                                                                                                                                                                                 tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGaaaa                                                   >  1:185128/1‑249 (MQ=255)
                                                                                                                                                                                                        gggagggaCAATGAGAAGAGGAGGGGAACCGACCCGGTTTTGTAATTAAGCTAAGACCTTTATTCACGGTGATGCGCTTGCCGAACTAAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGgca                                          <  2:388064/251‑1 (MQ=255)
                                                                                                                                                                                                         ggAGTGACCATGAGAACAGGATTTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGACCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGgcag                                         <  1:52649/251‑1 (MQ=255)
                                                                                                                                                                                                                    gagaACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGGAcatc                               >  1:20836/1‑250 (MQ=255)
                                                                                                                                                                                                                      aaaaaggaggTGAACCGCCCCGGTTGGGTAATTAAGCTAAGACCATTATTCACGGGGAGGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATg                            <  2:517778/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggTAATAAACCTAAGACCATTATGCACGGTGATGCTCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTAt  <  1:113403/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTAT  >  NC_000913/3411443‑3411933

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