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,6890GT97.6% 118.9 / ‑6.3 42E13* (GAA→TAA) yhdJDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (19/22);  minor base A (0/1);  total (19/23)
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.99e-01

AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATC  >  NC_000913/3411444‑3411934
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCt                                                                                                                                                                                                                                                  >  1:483513/1‑251 (MQ=255)
  aaCAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTaa                                                                                                                                                                                                                                                <  2:481798/251‑1 (MQ=255)
                 aCATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACgg                                                                                                                                                                                                                                 <  1:133872/251‑1 (MQ=255)
                  aaTGGGGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAAGTAAGCTAAGACCATTATTCACGGt                                                                                                                                                                                                                                <  2:319301/250‑1 (MQ=255)
                      gTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCGGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATg                                                                                                                                                                                                                            <  2:511044/251‑1 (MQ=255)
                        gATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATgcg                                                                                                                                                                                                                          <  1:24526/251‑1 (MQ=255)
                                acacCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATTAACGAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAAGTAAGTTAAGACCATTATTCTCGGTGATgagttttccg                                                                                                                                                                                                                   <  2:436424/250‑10 (MQ=255)
                                      gTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                                            <  1:174209/251‑1 (MQ=255)
                                        ggTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaa                                                                                                                                                                                                           <  2:177467/250‑1 (MQ=255)
                                            aCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTAGTCTCGGTGATGCGCTTGCCGAACTTAAAAAGAt                                                                                                                                                                                                      <  1:491567/251‑1 (MQ=255)
                                               aGACCCGTGCTGCGCTGATGATTGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATccc                                                                                                                                                                                                   <  1:434468/251‑1 (MQ=255)
                                                            gcTGATGCTGGTCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGAATGGGGAAGCGCCTTTTTTATAGGTTTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTTCCGAACTTAAAAAGATCCCCGCCGaacgggtc                                                                                                                                                                                      >  2:322132/1‑245 (MQ=255)
                                                                         aTCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTg                                                                                                                                                                         >  1:146363/1‑251 (MQ=255)
                                                                            aaCCGTGGTACGCTGCGTAAAAAAGTGAAAAAATTCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGGG‑ACCATGAGAACAGGATGTTAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCg                                                                                                                                                                      <  1:409121/251‑1 (MQ=255)
                                                                                           cGTAAAAAATTGCAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACGAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAACCTTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGAGTCGATCTGATCTTTGCCGACCCACCATATAAc                                                                                                                                                        >  2:185958/1‑250 (MQ=255)
                                                                                                  aaTTGAAAAAATACGGCATGAACTAATTAAGGTTAGCTAAATG‑TTGAATTAAAAAGGAGCTACTAGGCAGGGGGAAGGGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCGAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTaa                                                                                                                                                <  2:146363/251‑1 (MQ=255)
                                                                                                              aCGGAATGAACTAATTCAGGTTAGCAAAATGTTTGA‑TTAAAAAGGCGCAACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATg                                                                                                                                     <  2:201657/250‑1 (MQ=255)
                                                                                                                      aaCTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                                            >  2:165931/1‑251 (MQ=255)
                                                                                                                          atttCAGGTTAGCTAAATTATTGA‑TTAAAAAGGCGCTACGCGTCATGGGGACGCCGCTTTTTGATAGGGTTCACAAAGTGA‑TGTACCATGAAAACAGGAGGTGAACCGACCCGGTTTGGTAATTTAGCTAACACCATTATTCCCGGTGATGCGCTTGCCAAACTTAAAAAGATCCCCGCCGAAATTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATAGTCTGTTCGaa                                                                                                                          <  2:534518/247‑1 (MQ=255)
                                                                                                                          aagttagtttAGCTAAATGCTGGA‑TTAAAAAGGGGCTACTCGGCAGGGGTAAGTGCCTTTTTTAGAGGTGTCACAAAGGGAGTG‑CCCATGAGAACGGGATGTGACCCGACCGGGTTTGGTAATTAAGCTAAGACCATTATTTACGGTGATGCGCTTGCCGATCTTAAAAAGATCCCCTCCTAAAGTGTCTATCTGATCTTTGCCGACCCACCATATAACAGCGTGAAAAATTGTGATTGTCGGATCGAAGc                                                                                                                        <  2:253296/243‑1 (MQ=255)
                                                                                                                                   tAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaa                                                                                                                  >  1:24653/1‑248 (MQ=255)
                                                                                                                                         aaTGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGGTTCGATCTGCTCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCt                                                                                                         >  2:240236/1‑251 (MQ=255)
                                                                                                                                              ttGA‑TTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGttt                                                                                                     <  2:177236/250‑1 (MQ=255)
                                                                                                                                               tGA‑TAAAAAAGGCGCTACTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAATGGAGTG‑ACCATGAGAACAGGATGTGCACCGACCCGGTTTGGAAATTAAGCTAAGACAATTATTCACGGGGATGCGCTTGCCTAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCAGATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTa                                                                                                    <  2:516674/250‑1 (MQ=255)
                                                                                                                                                    tAAAAAGGCGCTACTAGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGAc                                                                                               >  2:373150/1‑251 (MQ=255)
                                                                                                                                                          ggCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTg                                                                                         >  2:225857/1‑251 (MQ=255)
                                                                                                                                                            cgcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAACGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGCCTGGCTGtt                                                                                       >  2:338234/1‑251 (MQ=255)
                                                                                                                                                            cgcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAACGGGAGTG‑ACCATGAGAACAGGATCTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGTTGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGCTCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGCAGCCTGGAAAGACGATCTGTTTATCGACTGGCTGtt                                                                                       >  2:232862/1‑251 (MQ=255)
                                                                                                                                                             gcTACGCGGCATGGGGAAGCGGCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAAAACAGTATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATGATTCACGGTGATGCGTTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTGTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCTACTGGCTGttt                                                                                      <  2:52753/251‑1 (MQ=255)
                                                                                                                                                                                 gcCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGcc                                                                   >  2:351418/1‑250 (MQ=255)
                                                                                                                                                                                      ttttATAGGTGTCACAAAGGGAGTG‑ACCAGGAGAACAGGATGTGAACCGACCCGGTTTGTTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACcgc                                                              <  1:516157/250‑1 (MQ=255)
                                                                                                                                                                                              gtgtCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGCTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATACCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGACGCCTGGAAAGAAGATCTGTTTATCGCCTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGaa                                                      >  1:94292/1‑250 (MQ=255)
                                                                                                                                                                                                      aaGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAggg                                             <  1:501203/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTAc                                   >  1:601103/1‑251 (MQ=255)
                                                                                                                                                                                                                     gagaACAGGATGTGAACCGACCCGGTTTGGTAATAAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTAcatca                               <  1:300576/251‑1 (MQ=255)
                                                                                                                                                                                                                       gaACAGGATGTGAACCGACCCGGGTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCAGACCCACAATATAACATCGGGAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCAGGTAcatcagg                             >  2:147104/1‑249 (MQ=255)
                                                                                                                                                                                                                           aGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAAcc                         >  2:593245/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGa                  >  1:366248/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGa                  >  1:65131/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aCCGACCCGTCTTGGTAATTACGCTAAGAACATTATTCACGGTGATGCGCTTGCCGTACTTCAAAAGAACCCCGCCGAAAGTGTCGATCTGATCTTAGCCGACCCCCCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGATGTGCTAGCAGAGAGCCACCGCGTTCTGAAACAGCAGAGCAGCATGTACATCATGAACAGTACGGaa                 >  2:429755/1‑250 (MQ=255)
                                                                                                                                                                                                                                                tGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCttt     <  1:552590/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAATTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGTAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATc  <  1:104738/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGA‑TTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTG‑ACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATC  >  NC_000913/3411444‑3411934

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