Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 434,375 A→G 5.7% A307A (GCA→GCG ribD → fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913434,3750AG 5.7% 215.8 / 3.0 70A307A (GCA→GCGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base G (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

GTGAGTGGCCGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATACGTCATGTA  >  NC_000913/434236‑434522
                                                                                                                                           |                                                                                                                                                   
gTGAGTTGCCGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTTCTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCa                                                                                                                                             <  2:1954/148‑1 (MQ=255)
    gTGGCCGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCgg                                                                                                                                        <  1:189540/149‑1 (MQ=255)
    gTGGCCGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCgg                                                                                                                                        <  5:166431/149‑1 (MQ=255)
         cGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATTATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGtt                                                                                                                                     >  8:310084/1‑147 (MQ=255)
         cGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGtt                                                                                                                                     <  4:48277/147‑1 (MQ=255)
         cGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGtt                                                                                                                                     >  3:48277/1‑147 (MQ=255)
         cGGAAACGGTGAGTACCTTGCTGATTCCATAGCATAAAGGTAATCTTGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTAGCTGGCGCATTGCTGCAGCCGGttt                                                                                                                                     >  1:339310/1‑147 (MQ=255)
           gAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCa                                                                                                                                                     <  7:210863/129‑1 (MQ=255)
           gAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCa                                                                                                                                                     >  8:210863/1‑129 (MQ=255)
            aaaCGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAATCAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTc                                                                                                                                >  1:264363/1‑149 (MQ=255)
              aCGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGa                                                                                                                              <  1:313491/149‑1 (MQ=255)
               cGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGAt                                                                                                                             <  3:330529/149‑1 (MQ=255)
                       aCCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCGTTGCTGCAGGCGGGTTTAGTCGATGAGCTGAt                                                                                                                     >  7:225524/1‑149 (MQ=255)
                       aCCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCGTTGCTGCAGGCGGGTTTAGTCGATGAGCTGAt                                                                                                                     >  7:230965/1‑149 (MQ=255)
                       aCCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGAt                                                                                                                     >  8:44462/1‑149 (MQ=255)
                             cTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                                                                                                               >  8:113049/1‑149 (MQ=255)
                              tGATTCCAGAGCATAAAGGTCATGTGGATCTGGTTGTACTGATGATGCAAATGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTGTAGTCGATGAGCTGATTGTCta                                                                                                               >  5:44411/1‑148 (MQ=255)
                               gATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtata                                                                                                             >  3:284682/1‑149 (MQ=255)
                               gATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGAGCATTGCTGCAGGCGTGTTTAGTCGATGAGCTGATTGTCtata                                                                                                             >  3:344296/1‑149 (MQ=255)
                                aTTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATTATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTatat                                                                                                            <  7:310084/149‑1 (MQ=255)
                                     agagCATAAAGGTCATCTGTATCTGGTTGTAATTATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCa                                                                                                        <  8:144737/148‑1 (MQ=255)
                                         cATAAACGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTaa                                                                                                   <  1:345153/149‑1 (MQ=255)
                                           tAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc                                                                                                 >  4:81756/1‑149 (MQ=255)
                                           tAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc                                                                                                 >  8:222041/1‑149 (MQ=255)
                                                 tCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAg                                                                                           >  3:264240/1‑149 (MQ=255)
                                                        gatctggTTGTACTGATGATGCAACTGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGAc                                                                                   >  8:44856/1‑149 (MQ=255)
                                                             ggTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGc                                                                                                                         <  6:230794/107‑1 (MQ=255)
                                                             ggTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGc                                                                                                                         >  5:230794/1‑107 (MQ=255)
                                                              gTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCcg                                                                              <  3:81756/149‑1 (MQ=255)
                                                              gTTGTACGTATCCTGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGAc                                                                                   <  6:44411/144‑1 (MQ=255)
                                                                    cTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCGTTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                         <  8:230965/148‑1 (MQ=255)
                                                                    cTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCGTTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAt                                                                         <  8:225524/148‑1 (MQ=255)
                                                                       atgatgCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATg                                                                     <  3:113389/149‑1 (MQ=255)
                                                                        tgatgCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGc                                                                    <  7:113049/149‑1 (MQ=255)
                                                                            gCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAtt                                                                        <  7:135275/141‑1 (MQ=255)
                                                                            gCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGAtt                                                                        >  8:135275/1‑141 (MQ=255)
                                                                            gCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGc                                                                >  1:6102/1‑149 (MQ=255)
                                                                              aaCTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGCGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTg                                                              >  3:283449/1‑149 (MQ=255)
                                                                              aaCTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTg                                                              >  2:155486/1‑149 (MQ=255)
                                                                                   ggTAAACAGCAAATTACCAGCATCTGTGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGATTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCa                                                           <  2:264363/147‑1 (MQ=255)
                                                                                   ggTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGTGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCTCGGATTATGCACGCTGCCa                                                           <  4:344296/147‑1 (MQ=255)
                                                                                   ggTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCa                                                           <  7:222041/147‑1 (MQ=255)
                                                                                          aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGa                                                                                                                              >  5:290914/1‑73 (MQ=255)
                                                                                          aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGa                                                                                                                              <  6:290914/73‑1 (MQ=255)
                                                                                           gCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTgaga                                                 <  1:144629/149‑1 (MQ=255)
                                                                                              aaTTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGTTTTAGTCGATGAGCTTATTGTCTATCTCGCACCTAAACTATTAGGCAGCGCCGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAAt                                              >  5:190087/1‑149 (MQ=255)
                                                                                              aaTTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAAt                                              >  7:221493/1‑149 (MQ=255)
                                                                                              aaTTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAAt                                              >  1:82918/1‑149 (MQ=255)
                                                                                               aTTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCAc                                                                                                       >  6:113247/1‑91 (MQ=255)
                                                                                               aTTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCAc                                                                                                       <  5:113247/91‑1 (MQ=255)
                                                                                               aTTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAAtt                                             <  2:82918/149‑1 (MQ=255)
                                                                                                ttAAGAGGATCTGGGTGGAAGCGGGGCCAACGCTCGGTGGCGCATTGCTGCAGTCGGGTGTACTCGATGAGCTGATTGTCTATATCGCACCTACACTATTTGGCAGCGACGCCCGCGGATTTTGCACGCTGCCAGGGCTTGAGAAATTa                                            <  2:6102/149‑1 (MQ=255)
                                                                                                   aCAGCATCTGGGGGGAAGCGGGGCCAACGCGCGCTGGCGCATTGCTGCAGGCGGGTGTAGTCGTGGCGCTGATTGTCTATATCGCACCTAAAGTATTAGGCAGCGACGGCCGCGGATTATGCACGGTGCCAGGGCTTGAGAAATTAGcc                                         <  2:272186/149‑1 (MQ=255)
                                                                                                    cAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCg                                        >  8:50688/1‑149 (MQ=255)
                                                                                                      gCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGAc                                      <  7:44462/149‑1 (MQ=255)
                                                                                                           tGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATAGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGCGCTTGAGAAATTAGCCGACGcccc                                 >  5:69719/1‑149 (MQ=255)
                                                                                                                   aGCGGGGCCAACACTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAAt                                              >  7:949/1‑128 (MQ=255)
                                                                                                                   aGCGGGGCCAACACTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAAt                                              <  8:949/128‑1 (MQ=255)
                                                                                                                    gCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGc                                                                                 <  8:220320/92‑1 (MQ=255)
                                                                                                                    gCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGc                                                                                 >  7:220320/1‑92 (MQ=255)
                                                                                                                      ggggCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAAt                      <  1:273548/149‑1 (MQ=255)
                                                                                                                        ggCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGGCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTc                    <  1:85361/149‑1 (MQ=255)
                                                                                                                         gCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTc                    >  2:341561/1‑148 (MQ=255)
                                                                                                                           cAACGCTCGCTGTCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCaaa                 <  1:341561/149‑1 (MQ=255)
                                                                                                                           cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCaaa                 <  8:221493/149‑1 (MQ=255)
                                                                                                                             aCGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAaga               <  7:44856/149‑1 (MQ=255)
                                                                                                                                   gctGGCGCATTGCTGGAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATa           <  4:284682/147‑1 (MQ=255)
                                                                                                                                   gctGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGc                                                       <  3:326419/103‑1 (MQ=255)
                                                                                                                                   gctGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGc                                                       >  4:326419/1‑103 (MQ=255)
                                                                                                                                     tGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATACGt        >  4:94999/1‑148 (MQ=255)
                                                                                                                                     gggCGCATTGCTGCAGGCGGGTGTAGTCGATGAGCTGATTGGCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATACGTc       <  2:337541/148‑1 (MQ=255)
                                                                                                                                          cATTGCTGCAGGCTGGTTTAGTCGATGAGCTGATTTTATATATCGCACCTAAAGTATTAGGCAGCGACGCCCGCGTATTATGCCCTCTGCCAGGGCTTGAGAAATTTGCCGACGCCCCCCAATTTAAATTCAAAGATATACGTCATGTa  <  8:352662/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GTGAGTGGCCGGAAACGGTGCGTACCTTGCTGATTCCAGAGCATAAAGGTCATCTGGATCTGGTTGTACTGATGATGCAACTGGGTAAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACTATTAGGCAGCGACGCCCGCGGATTATGCACGCTGCCAGGGCTTGAGAAATTAGCCGACGCCCCCCAATTTAAATTCAAAGAGATACGTCATGTA  >  NC_000913/434236‑434522

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