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,6890GT92.9% 116.6 / ‑3.2 43E13* (GAA→TAA) yhdJDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base T (20/19);  minor base A (0/3);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTT  >  NC_000913/3411446‑3411931
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aaCAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTa                                                                                                                                                                                                                                               >  1:474279/1‑250 (MQ=255)
     ccccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAcc                                                                                                                                                                                                                                          >  1:668542/1‑250 (MQ=255)
       ccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGAGATACTCGGCAAGGGGAAGAGACTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTTCACCGACCCGGTTTGGTAATTAAGCTAAgcacatt                                                                                                                                                                                                                                       >  1:124479/1‑245 (MQ=255)
          gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGATGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCattatt                                                                                                                                                                                                                                    <  2:159858/251‑1 (MQ=255)
              gACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACg                                                                                                                                                                                                                               <  2:356823/251‑1 (MQ=255)
                 aTGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTg                                                                                                                                                                                                                            <  1:546505/251‑1 (MQ=255)
                   ggTGATGCAATAAACCCGAGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAACTACGGCATGAAATAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGAGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACCAAGGGAGTGAGCATGAGAACAGGATGTGAACCGGCCCAGTTTGGTCATTAAGCTAAGACCATTATTAACgggga                                                                                                                                                                                                                           >  2:541367/1‑247 (MQ=255)
                    gTGATGCAAGTCACCCGTGGTAACCAGACCCGTGCTGCGCTGAGGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGGATGAACTAATTCAGGTTAGTTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTGTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACGGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATg                                                                                                                                                                                                                         <  2:588162/251‑1 (MQ=255)
                      gATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATgcg                                                                                                                                                                                                                       <  1:508039/251‑1 (MQ=255)
                       aTGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATgcgc                                                                                                                                                                                                                      >  1:23451/1‑251 (MQ=255)
                          cAATACACCCGTGGTAATCAGACCCGTGCTGCGCTGATGATGGGCATCTACCGTGGTACGCAGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGATTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCtt                                                                                                                                                                                                                    <  2:474279/250‑1 (MQ=255)
                                         aaCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAg                                                                                                                                                                                                     <  1:406241/250‑1 (MQ=255)
                                             aGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTAAGCTGCGTAAAAAATTGAAAAAATACGGCATGAACGAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGGCCATGAGAACAGGATGTGAACCGACCCGGTGTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATccc                                                                                                                                                                                                <  2:343506/251‑1 (MQ=255)
                                             aGACACGTGCTGCGATGATGATGGGCATCAACCGTTGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGAcccc                                                                                                                                                                                                <  2:110467/251‑1 (MQ=255)
                                                       tacgatgatgaTGGGCATCAACCGTGGTAAGCTGCGTAAAAAATTGGAAGAAGACGGCATGGACTAATTCAGGTTAGCTAAATGCAGGATTAAAAAGGGGCTACTCGGACTGGGGAATCGCCTTTTTTATAGGTGTCACATAAGGAGTGAACAAGAGAAAAGGGTATTAACCGGCCAGGTTTGGTAATGAAGCAAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAGAGATTCCCGCCGAAAgg                                                                                                                                                                                      >  2:355812/6‑250 (MQ=255)
                                                                          aaCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGAATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTTGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTc                                                                                                                                                                                   >  1:207377/1‑235 (MQ=255)
                                                                             cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                                                                                                                                >  1:217307/1‑251 (MQ=255)
                                                                                ggTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAACGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTCCCGACccac                                                                                                                                                             >  2:69697/1‑251 (MQ=255)
                                                                                     gCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                                         <  1:550444/202‑1 (MQ=255)
                                                                                     gCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                                         >  2:550444/1‑202 (MQ=255)
                                                                                         cGTAAAACATTCAAAAAATACGGCATGAACTAATTAAGGTTCGCTCAATGCTTGGTTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAAGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACAGACCCGGTTTGGACATTAAGCTAAGACCATTATTACCGGTGGAGCGCTTGCCGAACTAAAAAAGATACCCGCCGAAATGGGCGATCTGATCTTTGcaaacaccacataaacc                                                                                                                                                     >  2:619152/1‑234 (MQ=255)
                                                                                               aaaTTGAAAAAATACGGCATGAACTAATTCAGGATAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTAACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACCTCGGt                                                                                                                                               >  1:482192/1‑250 (MQ=255)
                                                                                                        aaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGGCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAAttttt                                                                                                                                     >  2:246210/1‑250 (MQ=255)
                                                                                                          aTACGGCATGAACTATTTCAGTTTAGCTAAATGCTTTATTAAAAGGGCGCTATTCGGCAGGGGGATGCGCCTTTTTAATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATTTGAACCGACCCGGTTTGGTAATTAACTTAAGACCATTTTTCACGTTGATTCGTTTGCCGACCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGAt                                                                                                                                   <  1:448111/251‑1 (MQ=255)
                                                                                                                    aaCTAATTAAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATCCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                                         <  1:561209/251‑1 (MQ=255)
                                                                                                                                                    aaaGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTgg                                                                                         >  1:175754/1‑251 (MQ=255)
                                                                                                                                                      aGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCt                                                                                       >  1:501897/1‑251 (MQ=255)
                                                                                                                                                                   cagggggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGAtt                                                                          <  1:246210/248‑1 (MQ=255)
                                                                                                                                                                      ggggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAACAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACCTCGGTACAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGCAGTGATTGc                                                                        >  1:100320/1‑250 (MQ=255)
                                                                                                                                                                      ggcgaaccgcCTTTTTTATAGGTGTCCCAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTCAGCTAAGACCATAATTCACGGTTATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGAACTTAGCCGACCCACCATATAACATCGGTAAAACTTTAGAAGGTATGTTCGAATCCAGGAAAGCAGATCTGTTTATCGCCTGGTCGTTTGGAGTgctttc                                                                        >  2:233130/8‑245 (MQ=255)
                                                                                                                                                                        ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAATACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAAc                                                                                                                                                     <  1:391251/171‑1 (MQ=255)
                                                                                                                                                                        ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAAc                                                                                                                                                     >  2:391251/1‑171 (MQ=255)
                                                                                                                                                                                  tttttATAGTTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATCCACGGTGATGCGTTTGCCGAACTAAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACcg                                                            <  1:262362/250‑1 (MQ=255)
                                                                                                                                                                                       atatgggtCAAAAAGGGAGTGACTATGACAACGGGTTGGGAACCGACCTGTTTGGGGAATAAAGCGAAGACTATTATTCACGTTGATGCGTTTGCCGAACGTAAAAAGATCACCGCCGAAAGGGTCGATCCGATCATTGCCGACCCACCATATAACATCGGTCAAAACTTTGATGGTCTGATCGAAGCCTGGAAAGAATATCTGTTTATCGACTGGCTGTTTGAAGTGATTTCAGAGTGCCAACGCGTTCt                                                      <  2:173462/245‑1 (MQ=255)
                                                                                                                                                                                           gtgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGCCCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGCTCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACAGGCTGTTTGAAGGGATTGCAGAGTGACAACGCGTTCttaa                                                   >  2:287782/1‑247 (MQ=255)
                                                                                                                                                                                            gcttccaAAGGGAGTGACCATGATAACAGGATGTGAACCGACCCTGTTTGTATATAAAGCTAAGCCCACTATCAGCGGTGATTCGCTTTCTGAATTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGTCCCACCATATAACATCGTTAAAAATGTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCGGTTTATCGACTGGCTGTTTGAAGTGATTGTAGAGGGCCACCGCGTTCTTaaaa                                                 <  2:73231/246‑1 (MQ=255)
                                                                                                                                                                                                             ccATGAAAACAGGATGTGAACCGCCCCGGTTTGGTAATAAAGCTAAGTCCATTATTCACGGTGATGCGCTTGCCGAACTAAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTAc                                <  1:541367/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaCAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGACCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGa                         <  1:180418/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGaa              >  1:53900/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gccccGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGTTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACa          <  1:69697/249‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCAAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGcc      <  1:131433/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTTGGAAATTAAGCTAAGACCATTATTCACGGTGCTGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAACGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGAGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGCACAGAACGGAAAACATgacc     >  2:421695/1‑248 (MQ=255)
                                                                                                                                                                                                                                            tGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCttt  <  1:158497/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
AACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTT  >  NC_000913/3411446‑3411931

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