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

AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATC  >  NC_000913/3411444‑3411934
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCt                                                                                                                                                                                                                                                  >  1:483513/1‑251 (MQ=255)
  aaCAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTaa                                                                                                                                                                                                                                                <  2:481798/251‑1 (MQ=255)
                 aCATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACgg                                                                                                                                                                                                                                 <  1:133872/251‑1 (MQ=255)
                  aaTGGGGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAAGTAAGCTAAGACCATTATTCACGGt                                                                                                                                                                                                                                <  2:319301/250‑1 (MQ=255)
                      gTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCGGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATg                                                                                                                                                                                                                            <  2:511044/251‑1 (MQ=255)
                        gATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATgcg                                                                                                                                                                                                                          <  1:24526/251‑1 (MQ=255)
                                acacCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATTAACGAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAAGTAAGTTAAGACCATTATTCTCGGTGATgagttttccg                                                                                                                                                                                                                   <  2:436424/250‑10 (MQ=255)
                                      gTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                                            <  1:174209/251‑1 (MQ=255)
                                        ggTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaa                                                                                                                                                                                                           <  2:177467/250‑1 (MQ=255)
                                            aCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTAGTCTCGGTGATGCGCTTGCCGAACTTAAAAAGAt                                                                                                                                                                                                      <  1:491567/251‑1 (MQ=255)
                                               aGACCCGTGCTGCGCTGATGATTGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATccc                                                                                                                                                                                                   <  1:434468/251‑1 (MQ=255)
                                                            gcTGATGCTGGTCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGAATGGGGAAGCGCCTTTTTTATAGGTTTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTTCCGAACTTAAAAAGATCCCCGCCGaacgggtc                                                                                                                                                                                      >  2:322132/1‑245 (MQ=255)
                                                                         aTCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTg                                                                                                                                                                         >  1:146363/1‑251 (MQ=255)
                                                                            aaCCGTGGTACGCTGCGTAAAAAAGTGAAAAAATTCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGGGACCATGAGAACAGGATGTTAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCg                                                                                                                                                                      <  1:409121/251‑1 (MQ=255)
                                                                                           cGTAAAAAATTGCAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACGAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAACCTTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGAGTCGATCTGATCTTTGCCGACCCACCATATAAc                                                                                                                                                        >  2:185958/1‑250 (MQ=255)
                                                                                                  aaTTGAAAAAATACGGCATGAACTAATTAAGGTTAGCTAAATGTTGAATTAAAAAGGAGCTACTAGGCAGGGGGAAGGGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCGAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTaa                                                                                                                                                <  2:146363/251‑1 (MQ=255)
                                                                                                              aCGGAATGAACTAATTCAGGTTAGCAAAATGTTTGATTAAAAAGGCGCAACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATg                                                                                                                                     <  2:201657/250‑1 (MQ=255)
                                                                                                                      aaCTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                                            >  2:165931/1‑251 (MQ=255)
                                                                                                                          atttCAGGTTAGCTAAATTATTGATTAAAAAGGCGCTACGCGTCATGGGGACGCCGCTTTTTGATAGGGTTCACAAAGTGATGTACCATGAAAACAGGAGGTGAACCGACCCGGTTTGGTAATTTAGCTAACACCATTATTCCCGGTGATGCGCTTGCCAAACTTAAAAAGATCCCCGCCGAAATTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATAGTCTGTTCGaa                                                                                                                          <  2:534518/247‑1 (MQ=255)
                                                                                                                          aagttagtttAGCTAAATGCTGGATTAAAAAGGGGCTACTCGGCAGGGGTAAGTGCCTTTTTTAGAGGTGTCACAAAGGGAGTGCCCATGAGAACGGGATGTGACCCGACCGGGTTTGGTAATTAAGCTAAGACCATTATTTACGGTGATGCGCTTGCCGATCTTAAAAAGATCCCCTCCTAAAGTGTCTATCTGATCTTTGCCGACCCACCATATAACAGCGTGAAAAATTGTGATTGTCGGATCGAAGc                                                                                                                        <  2:253296/243‑1 (MQ=255)
                                                                                                                                   tAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaa                                                                                                                  >  1:24653/1‑248 (MQ=255)
                                                                                                                                         aaTGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGGTTCGATCTGCTCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCt                                                                                                         >  2:240236/1‑251 (MQ=255)
                                                                                                                                              ttGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGttt                                                                                                     <  2:177236/250‑1 (MQ=255)
                                                                                                                                               tGATAAAAAAGGCGCTACTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAATGGAGTGACCATGAGAACAGGATGTGCACCGACCCGGTTTGGAAATTAAGCTAAGACAATTATTCACGGGGATGCGCTTGCCTAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCAGATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTa                                                                                                    <  2:516674/250‑1 (MQ=255)
                                                                                                                                                    tAAAAAGGCGCTACTAGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGAc                                                                                               >  2:373150/1‑251 (MQ=255)
                                                                                                                                                          ggCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTg                                                                                         >  2:225857/1‑251 (MQ=255)
                                                                                                                                                            cgcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAACGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGCCTGGCTGtt                                                                                       >  2:338234/1‑251 (MQ=255)
                                                                                                                                                            cgcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAACGGGAGTGACCATGAGAACAGGATCTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGTTGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGCTCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGCAGCCTGGAAAGACGATCTGTTTATCGACTGGCTGtt                                                                                       >  2:232862/1‑251 (MQ=255)
                                                                                                                                                             gcTACGCGGCATGGGGAAGCGGCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAAAACAGTATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATGATTCACGGTGATGCGTTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTGTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCTACTGGCTGttt                                                                                      <  2:52753/251‑1 (MQ=255)
                                                                                                                                                                                 gcCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGcc                                                                   >  2:351418/1‑250 (MQ=255)
                                                                                                                                                                                      ttttATAGGTGTCACAAAGGGAGTGACCAGGAGAACAGGATGTGAACCGACCCGGTTTGTTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACcgc                                                              <  1:516157/250‑1 (MQ=255)
                                                                                                                                                                                              gtgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGCTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATACCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGACGCCTGGAAAGAAGATCTGTTTATCGCCTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGaa                                                      >  1:94292/1‑250 (MQ=255)
                                                                                                                                                                                                      aaGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAggg                                             <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTTATC  >  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: