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

AACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATGAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTT  >  NC_000913/3411446‑3411931
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aaCAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTa                                                                                                                                                                                                                                               >  1:474279/1‑250 (MQ=255)
     ccccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGAcc                                                                                                                                                                                                                                          >  1:668542/1‑250 (MQ=255)
       ccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGAGATACTCGGCAAGGGGAAGAGACTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTTCACCGACCCGGTTTGGTA‑ATTAAGCTAAgcacatt                                                                                                                                                                                                                                       >  1:124479/1‑245 (MQ=255)
          gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGATGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCattatt                                                                                                                                                                                                                                    <  2:159858/251‑1 (MQ=255)
              gACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑Cg                                                                                                                                                                                                                               <  2:356823/251‑1 (MQ=255)
                 aTGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTg                                                                                                                                                                                                                            <  1:546505/251‑1 (MQ=255)
                   ggTGATGCAATAAACCCGAGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAACTACGGCATGAAATAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGAGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACCAAGGGAGTGAGCATGAGAACAGGATGTGAACCGGCCCAGTTTGGTC‑ATTAAGCTAAGACCATTATTAA‑Cgggga                                                                                                                                                                                                                           >  2:541367/1‑247 (MQ=255)
                    gTGATGCAAGTCACCCGTGGTAACCAGACCCGTGCTGCGCTGAGGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGGATGAACTAATTCAGGTTAGTTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTGTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACGGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATg                                                                                                                                                                                                                         <  2:588162/251‑1 (MQ=255)
                      gATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATgcg                                                                                                                                                                                                                       <  1:508039/251‑1 (MQ=255)
                       aTGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATgcgc                                                                                                                                                                                                                      >  1:23451/1‑251 (MQ=255)
                          cAATACACCCGTGGTAATCAGACCCGTGCTGCGCTGATGATGGGCATCTACCGTGGTACGCAGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGATTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCtt                                                                                                                                                                                                                    <  2:474279/250‑1 (MQ=255)
                                         aaCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAg                                                                                                                                                                                                     <  1:406241/250‑1 (MQ=255)
                                             aGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTAAGCTGCGTAAAAAATTGAAAAAATACGGCATGAACGAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGGCCATGAGAACAGGATGTGAACCGACCCGGTGTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATccc                                                                                                                                                                                                <  2:343506/251‑1 (MQ=255)
                                             aGACACGTGCTGCGATGATGATGGGCATCAACCGTTGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGAcccc                                                                                                                                                                                                <  2:110467/251‑1 (MQ=255)
                                                       tacgatgatgaTGGGCATCAACCGTGGTAAGCTGCGTAAAAAATTGGAAGAAGACGGCATGGACTAATTCAGGTTAGCTAAATGCAGGATTAAAAAGGGGCTACTCGGACTGGGGAATCGCCTTTTTTATAGGTGTCACATAAGGAGTGAACAAGAGAAAAGGGTATTAACCGGCCAGGTTTGGTA‑ATGAAGCAAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAGAGATTCCCGCCGAAAgg                                                                                                                                                                                      >  2:355812/6‑250 (MQ=255)
                                                                          aaCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGAATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTTGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTc                                                                                                                                                                                   >  1:207377/1‑235 (MQ=255)
                                                                             cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                                                                                                                                >  1:217307/1‑251 (MQ=255)
                                                                                ggTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAACGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTCCCGACccac                                                                                                                                                             >  2:69697/1‑251 (MQ=255)
                                                                                     gCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                                         <  1:550444/202‑1 (MQ=255)
                                                                                     gCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                                         >  2:550444/1‑202 (MQ=255)
                                                                                         cGTAAAACATTCAAAAAATACGGCATGAACTAATTAAGGTTCGCTCAATGCTTGGTTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAAGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACAGACCCGGTTTGGAC‑ATTAAGCTAAGACCATTATTAC‑CGGTGGAGCGCTTGCCGAACTAAAAAAGATACCCGCCGAAATGGGCGATCTGATCTTTGcaaacaccacataaacc                                                                                                                                                     >  2:619152/1‑234 (MQ=255)
                                                                                               aaaTTGAAAAAATACGGCATGAACTAATTCAGGATAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTAA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACCTCGGt                                                                                                                                               >  1:482192/1‑250 (MQ=255)
                                                                                                        aaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGGCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAAttttt                                                                                                                                     >  2:246210/1‑250 (MQ=255)
                                                                                                          aTACGGCATGAACTATTTCAGTTTAGCTAAATGCTTTATTAAAAGGGCGCTATTCGGCAGGGGGATGCGCCTTTTTAATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATTTGAACCGACCCGGTTTGGTA‑ATTAACTTAAGACCATTTTTCA‑CGTTGATTCGTTTGCCGACCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGAt                                                                                                                                   <  1:448111/251‑1 (MQ=255)
                                                                                                                    aaCTAATTAAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATCCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                                         <  1:561209/251‑1 (MQ=255)
                                                                                                                                                    aaaGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTgg                                                                                         >  1:175754/1‑251 (MQ=255)
                                                                                                                                                      aGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCt                                                                                       >  1:501897/1‑251 (MQ=255)
                                                                                                                                                                   cagggggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGAtt                                                                          <  1:246210/248‑1 (MQ=255)
                                                                                                                                                                      ggggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAACAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACCTCGGTACAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGCAGTGATTGc                                                                        >  1:100320/1‑250 (MQ=255)
                                                                                                                                                                      ggcgaaccgcCTTTTTTATAGGTGTCCCAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTCAGCTAAGACCATAATTCA‑CGGTTATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGAACTTAGCCGACCCACCATATAACATCGGTAAAACTTTAGAAGGTATGTTCGAATCCAGGAAAGCAGATCTGTTTATCGCCTGGTCGTTTGGAGTgctttc                                                                        >  2:233130/8‑245 (MQ=255)
                                                                                                                                                                        ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAATACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAAc                                                                                                                                                     <  1:391251/171‑1 (MQ=255)
                                                                                                                                                                        ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAAc                                                                                                                                                     >  2:391251/1‑171 (MQ=255)
                                                                                                                                                                                  tttttATAGTTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATCCA‑CGGTGATGCGTTTGCCGAACTAAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACcg                                                            <  1:262362/250‑1 (MQ=255)
                                                                                                                                                                                       atatgggtCAAAAAGGGAGTGACTATGACAACGGGTTGGGAACCGACCTGTTTGGGGA‑ATAAAGCGAAGACTATTATTCA‑CGTTGATGCGTTTGCCGAACGTAAAAAGATCACCGCCGAAAGGGTCGATCCGATCATTGCCGACCCACCATATAACATCGGTCAAAACTTTGATGGTCTGATCGAAGCCTGGAAAGAATATCTGTTTATCGACTGGCTGTTTGAAGTGATTTCAGAGTGCCAACGCGTTCt                                                      <  2:173462/245‑1 (MQ=255)
                                                                                                                                                                                           gtgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGCCCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGCTCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACAGGCTGTTTGAAGGGATTGCAGAGTGACAACGCGTTCttaa                                                   >  2:287782/1‑247 (MQ=255)
                                                                                                                                                                                            gcttccaAAGGGAGTGACCATGATAACAGGATGTGAACCGACCCTGTTTG‑TATATAAAGCTAAGCCCACTAT‑CAGCGGTGATTCGCTTTCTGAATTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGTCCCACCATATAACATCGTTAAAAATGTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCGGTTTATCGACTGGCTGTTTGAAGTGATTGTAGAGGGCCACCGCGTTCTTaaaa                                                 <  2:73231/246‑1 (MQ=255)
                                                                                                                                                                                                             ccATGAAAACAGGATGTGAACCGCCCCGGTTTGGTA‑ATAAAGCTAAGTCCATTATTCA‑CGGTGATGCGCTTGCCGAACTAAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTAc                                <  1:541367/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaCAGGATGTGAACCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGACCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGa                         <  1:180418/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCCGACCCGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGaa              >  1:53900/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gccccGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGTTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACa          <  1:69697/249‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGTTTGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCAAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGcc      <  1:131433/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTTGGAA‑ATTAAGCTAAGACCATTATTCA‑CGGTGCTGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAACGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGAGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGCACAGAACGGAAAACATgacc     >  2:421695/1‑248 (MQ=255)
                                                                                                                                                                                                                                            tGGTA‑ATTAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCttt  <  1:158497/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
AACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTA‑ATGAAGCTAAGACCATTATTCA‑CGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTTT  >  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: ‑