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,6890GT98.0% 162.0 / NA 50E13* (GAA→TAA) yhdJDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (23/26);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTA  >  NC_000913/3411444‑3411932
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
aGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGCCCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGc                                                                                                                                                                                                                                                 <  2:482998/250‑1 (MQ=255)
   aCAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTa                                                                                                                                                                                                                                               >  1:57663/1‑249 (MQ=255)
      gCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAcc                                                                                                                                                                                                                                          <  2:443586/251‑1 (MQ=255)
        cccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCat                                                                                                                                                                                                                                        >  1:123796/1‑251 (MQ=255)
             tttgACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATTGGCATCAACCGTGGTACGCTTCGTAAAAAATTGAAAAAATACGGCATGAACGAATTCAGGTTATCTAAATGCTTGATTAAAAAGGCGCGTCTCGGCATGGGGAAGTGCCTTTTTTATAGGTTTTGCAATGGGAGTGTCCATAATAACAGTATGTGATCCGGCCCGGTTTGGTAATTAAGCTAAGACCattatt                                                                                                                                                                                                                                    <  1:27698/247‑1 (MQ=255)
               ggACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATGAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCAc                                                                                                                                                                                                                                 <  2:437941/251‑1 (MQ=255)
                gACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACg                                                                                                                                                                                                                                <  2:207421/251‑1 (MQ=255)
                gACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACg                                                                                                                                                                                                                                <  1:451203/251‑1 (MQ=255)
                  cATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTAGTCACGGt                                                                                                                                                                                                                              <  1:513943/251‑1 (MQ=255)
                     ggTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGa                                                                                                                                                                                                                            <  2:154902/250‑1 (MQ=255)
                             aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATTAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGc                                                                                                                                                                                                                   >  1:80800/1‑251 (MQ=255)
                                      gTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCTCTTGCCGAACTTaa                                                                                                                                                                                                          <  1:185142/251‑1 (MQ=255)
                                        ggTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaa                                                                                                                                                                                                         >  2:555036/1‑250 (MQ=255)
                                                   ccGTGCTGCGCTGATGATGGGCATCAACCGTGGTCCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCGTGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCccgcc                                                                                                                                                                                             <  1:18356/251‑1 (MQ=255)
                                                             gtgatgaTGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCg                                                                                                                                                                                   <  2:57663/250‑1 (MQ=255)
                                                                aggatgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACAAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCgatc                                                                                                                                                                                <  1:19040/249‑1 (MQ=255)
                                                                        catcagcggTGGTACGCGGCGTAAAAAATTGAAAAAATACGGAAGAAACTAATCAAGGTTAGCTAAATGTTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGCCCATGAAAACAGGATGTGATCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGTTTGCCGAACTTAAAAAGACCCCCGCCGAAAGTGTCGATCTGATCttt                                                                                                                                                                        <  2:210939/243‑1 (MQ=255)
                                                                         aTCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGCACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTg                                                                                                                                                                       >  1:83884/1‑251 (MQ=255)
                                                                            aaCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCg                                                                                                                                                                    >  2:43231/1‑251 (MQ=255)
                                                                                     aCGCTGCGTAAAAAATTGAAAAAATACGGCATGAACAAATTCAGGTTAGCTAAATGTTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGGGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCAGTATTCACGGTGATGCGCTTGCCGAACTAAAAAAGATCCCCGCCGAAAGTGGCGATCTGATCTTGGCCGACccacca                                                                                                                                                            <  2:287656/250‑1 (MQ=255)
                                                                                            gTAAAAAAGTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCAAAGCCAATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACAt                                                                                                                                                    <  2:177275/251‑1 (MQ=255)
                                                                                               aaaaaGTGAAAAAATACGGCAGAAACTAATTCAGGTTAGCTAAATGGTTGATTAAAAGGGCGGTCCGCGGAATGGGGAAGCGCGTTTTTGATAGGGGACAGAAGGGGGGTGACAAGGAGAACAGGGTGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCTCTTGCCGAACTTAAAAATACCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCgg                                                                                                                                                 <  2:536259/251‑1 (MQ=255)
                                                                                                      gAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAg                                                                                                                                                                                        >  2:472954/1‑205 (MQ=255)
                                                                                                      aaaaaaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAg                                                                                                                                                                                        <  1:472954/204‑1 (MQ=255)
                                                                                                          aaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTg                                                                                                                                      >  2:543301/1‑251 (MQ=255)
                                                                                                             tACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCTGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATg                                                                                                                                   >  2:408061/1‑251 (MQ=255)
                                                                                                              aCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATg                                                                                                                                   <  1:243042/250‑1 (MQ=255)
                                                                                                                   aTGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTg                                                                                                                              >  1:93941/1‑250 (MQ=255)
                                                                                                                          aaTTCAGGTTAGCTAAATGCTTGATTAAAACGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGCGACCATGAGCACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAACATTATTCACGGTGATGAGCTTGCCGAACTTAAAAAGATCCCCGCCGAACGGGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTCACAATTTTGATGCTCTGATCGAAGc                                                                                                                      >  2:167199/1‑251 (MQ=255)
                                                                                                                          aaTTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAg                                                                                                                       <  1:439173/250‑1 (MQ=255)
                                                                                                                           aTTCAGGACAGCTAAATGCTTGATTAAAAAGCAGCTCCTCGGCATGGGGAAGCGCATTTTTTATAGGTGTCACAAAGGCAGTGACAATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCAGACCCACCATATAACATCGGTAAAAATTTTGATGGTATGATCGAAGcc                                                                                                                     >  2:352811/1‑251 (MQ=255)
                                                                                                                              aaGGTTAGCTAAAGGCTGGATTAAAAGGGGGCGACTCGGCAGGGGGAAGCGCATTTTTTATGGGTGGCACAAGGGGAGTGACCATGAGAAAAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTGTGATGGTCTGATCGAAGCCTgg                                                                                                                  <  2:429215/250‑1 (MQ=255)
                                                                                                                               aGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTgg                                                                                                                  >  2:48012/1‑250 (MQ=255)
                                                                                                                               aGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGa                                                                                                                 <  1:337318/251‑1 (MQ=255)
                                                                                                                                  ttAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaa                                                                                                                >  2:485091/1‑249 (MQ=255)
                                                                                                                                   tAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaa                                                                                                                >  2:414790/1‑248 (MQ=255)
                                                                                                                                        aaaTGTTTGATAAAAAAGGCGCGACTCGTCATGGGGAAGCGCCTTTTTTATGGGTGTTACAAAGGGAGTGCCCATGAGAACAGGATGGGAAGCTACCCGTTTTGGTAATTAAGCTAAGACCATTATTCACGGTGAGGCGCTGGCCGAACTTAAAAAGATCCCCTCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAATAACGTGGATGTTCTGATCGAATCCTGGAAAGGAGATc                                                                                                        <  2:239218/251‑1 (MQ=255)
                                                                                                                                         aaTGTTTGATTAAAAAGGCGCTACGCGGCATGGGGACGCGCTTTTTTTATAGGTGTCACAAAGGGAGTGACCAGGAGAACAGAATGTGAACCAACCGGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTGTGATGGTCTGATCGAAGCCTGGAAAGAAGATCt                                                                                                       <  2:262171/251‑1 (MQ=255)
                                                                                                                                                gATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATc                                                                                                >  2:44664/1‑251 (MQ=255)
                                                                                                                                                   tAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGAc                                                                                             >  2:275276/1‑251 (MQ=255)
                                                                                                                                                      aaaGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTgg                                                                                          <  1:234957/251‑1 (MQ=255)
                                                                                                                                                               ccTCGGAATGGGGAAGCGCTTTTTTTATAGGGGTCACAAAGGAAGTGACCATGAGACCAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTGGCATAACTTAAAAATATCCCCGCCGAAAGTGTCGATCTGATCGTTGCCTACCCACCATATAACATCGGTAAAAATTTTGATGGTCTTATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGaa                                                                                 <  1:324617/250‑1 (MQ=255)
                                                                                                                                                                            aaGCGCCTTTTTTCTAGGTGTCACACAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGCCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGt                                                                    >  2:342802/1‑251 (MQ=255)
                                                                                                                                                                                                             gACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGt                                   <  1:555036/251‑1 (MQ=255)
                                                                                                                                                                                                                cATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTAc                                 >  2:16513/1‑250 (MQ=255)
                                                                                                                                                                                                                          ggATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAAc                        >  1:455427/1‑249 (MQ=255)
                                                                                                                                                                                                                              gtgAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACg                  >  1:280195/1‑251 (MQ=255)
                                                                                                                                                                                                                               tgAACCGACCCGGGTTGGTAATTAAGCGAAGACCATTATTCAAGGGGATGCGATTGCAGAACTTAAAACGATCCCCGACGAAAGAGTCGATCTGATATTTGCCGAACCCCCATATACCATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGGAGATCTGTTTATCGAATGGCTGTTTGAAGTGATTGCCGAGTGCCACCGCGTTCTGACAAAGAAGGGCAGCAAGAACATCATGAACAGTACgg                 >  1:542252/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGACCCGGTTTGGTAATTACGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAAc            >  2:65508/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCAAACTTAAAAAGAACCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTa  <  2:75464/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTA  >  NC_000913/3411444‑3411932

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