Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,798 A→G 24.2% S118P (TCT→CCT)  nudE ← adenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,525,7980AG24.2% 280.2 / 102.2 177S118P (TCT→CCT) nudEadenosine nucleotide hydrolase; Ap3A/Ap2A/ADP‑ribose/NADH hydrolase
Reads supporting (aligned to +/- strand):  ref base A (74/58);  new base G (25/17);  total (99/78)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

Only 100 of 181 total aligned reads displayed.

GCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACC  >  NC_000913/3525699‑3525897
                                                                                                   |                                                                                                   
gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGa                                                                                                     <  2:702246/100‑1 (MQ=255)
  cagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCAATTTGCTGGAAAAGTAAGACg                                                                                                   >  1:3782809/1‑100 (MQ=255)
   agcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACgg                                                                                                  <  1:2211394/100‑1 (MQ=255)
       cACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGcc                                                                                              <  1:3908044/100‑1 (MQ=255)
        aCCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCa                                                                                             <  1:2670458/100‑1 (MQ=255)
          cTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGACATg                                                                                           >  2:1692909/1‑100 (MQ=255)
           tgtgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTATCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGc                                                                                          <  2:2327412/100‑1 (MQ=255)
             ttgTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTTCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTg                                                                                        <  1:1006491/98‑1 (MQ=255)
              ggTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGa                                                                                       >  2:3295814/1‑100 (MQ=255)
              ggTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGa                                                                                       >  1:2177941/1‑100 (MQ=255)
               gTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAg                                                                                      <  2:2083208/100‑1 (MQ=255)
                  gCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGctt                                                                                   >  2:1175762/1‑100 (MQ=255)
                   cGGCTCTGGCTCATCGCCTTCCAGTGATTTCGGGTAGTGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGcttc                                                                                  <  1:2111541/100‑1 (MQ=255)
                   cGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGcttc                                                                                  >  2:283467/1‑100 (MQ=255)
                    ggCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGTCGGTGCCATGCTGAGcttc                                                                                  <  2:590256/99‑1 (MQ=255)
                     gCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGTAAAAGTAAGGCGGTGCCATGCTAAGctt                                                                                   >  1:956149/1‑97 (MQ=255)
                     gCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAATTAAGACGGTGCCATGCTGAGcttctt                                                                                <  1:2375320/100‑1 (MQ=255)
                      ctctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCttcttc                                                                               <  2:3122789/100‑1 (MQ=255)
                        ctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCaa                                                                             >  1:3064234/1‑100 (MQ=255)
                         tGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCaaa                                                                            >  2:3474275/1‑100 (MQ=255)
                         tGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCaaa                                                                            >  1:1465540/1‑100 (MQ=255)
                          ggCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCaaaa                                                                           <  1:2064291/100‑1 (MQ=255)
                          ggCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCaaaa                                                                           >  1:1482909/1‑100 (MQ=255)
                           gCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCaaaaa                                                                          <  2:823461/100‑1 (MQ=255)
                           gCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCaaaaa                                                                          >  1:2953556/1‑100 (MQ=255)
                            cTCATCTTCTTCCTGTGATTCCGGGTAGTGATCTTGCGCTTCCACGATATTCATTTTGCTGGAAAAGTAAGCCGGTGCCATGCTGAGCTTCTTCaaaaa                                                                          <  1:3031746/99‑1 (MQ=255)
                            cTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCaaaaaa                                                                         >  1:3884099/1‑100 (MQ=255)
                             tCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAg                                                                        <  1:3480177/100‑1 (MQ=255)
                              cATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGt                                                                       >  2:41800/1‑100 (MQ=255)
                                   ccTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGAt                                                                  <  2:2251795/100‑1 (MQ=255)
                                   ccTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGAt                                                                  <  2:1235099/100‑1 (MQ=255)
                                   ccTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGAt                                                                  <  2:2650848/100‑1 (MQ=255)
                                    cTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGAt                                                                  <  2:166451/99‑1 (MQ=255)
                                    cTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATc                                                                 >  1:3840391/1‑100 (MQ=255)
                                      tCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGt                                                               >  2:1902171/1‑100 (MQ=255)
                                      tCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGt                                                               >  2:844303/1‑100 (MQ=255)
                                      tCCAGTGATTCCGGGTAGAGATATTGCGCTACCACGATATTCATATTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCg                                                                >  2:958190/1‑99 (MQ=255)
                                       ccAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGtt                                                              >  1:2541617/1‑100 (MQ=255)
                                       ccAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGtt                                                              >  1:248660/1‑100 (MQ=255)
                                       ccAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGtt                                                              >  1:1082287/1‑100 (MQ=255)
                                        cAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTc                                                             >  2:1911836/1‑100 (MQ=255)
                                        cAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTc                                                             >  1:1575277/1‑100 (MQ=255)
                                         aGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGATGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGttag                                                            >  2:3377432/1‑98 (MQ=255)
                                          gTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGc                                                           >  2:236075/1‑100 (MQ=255)
                                            gATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCtn                                                         <  2:1436248/100‑2 (MQ=255)
                                             aTTCCGGGTAGAGATCTTGCGCTACCACTATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTc                                                         <  1:912302/99‑1 (MQ=255)
                                             aTTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTcc                                                        <  2:712044/100‑1 (MQ=255)
                                              ttCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCa                                                       >  1:379188/1‑100 (MQ=255)
                                              ttCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCa                                                       >  2:1543141/1‑100 (MQ=255)
                                                ccGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCaaa                                                     <  2:3840352/100‑1 (MQ=255)
                                                 cGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAAt                                                    >  2:3411295/1‑100 (MQ=255)
                                                   ggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATcc                                                  <  2:1058478/100‑1 (MQ=255)
                                                   ggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATcc                                                  <  2:2442545/100‑1 (MQ=255)
                                                    gTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATcc                                                  <  2:1031453/99‑1 (MQ=255)
                                                     tAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCaa                                                <  2:2971395/100‑1 (MQ=255)
                                                      agagATCTTGCGCTACCACGATATTCATTTTGCTTGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAAc                                               <  1:1309970/100‑1 (MQ=255)
                                                      agagATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAAc                                               <  1:1997027/100‑1 (MQ=255)
                                                       gagaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAAcc                                              <  2:3639174/100‑1 (MQ=255)
                                                       gagaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGATGAGCTTATTCAAAAAAATCAGATCGTTCGCTCCAAATCCAAcc                                              >  1:1519323/1‑100 (MQ=255)
                                                         gaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACctc                                            >  2:2526945/1‑100 (MQ=255)
                                                         gaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACctc                                            >  1:50543/1‑100 (MQ=255)
                                                           tCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCtctt                                          >  2:2019392/1‑100 (MQ=255)
                                                           tCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCtctt                                          >  2:9051/1‑100 (MQ=255)
                                                            cTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCtcttc                                         >  2:2658832/1‑100 (MQ=255)
                                                              tGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTcttctt                                       >  2:2901587/1‑100 (MQ=255)
                                                                cgcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCtttt                                     <  2:3039742/100‑1 (MQ=255)
                                                                 gcTACCATGATATTCATTTTTCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTa                                    <  1:862421/100‑1 (MQ=255)
                                                                   tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGc                                  >  1:2581012/1‑100 (MQ=255)
                                                                    aCCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCt                                 >  1:641653/1‑100 (MQ=255)
                                                                    aCCACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCt                                 >  2:1807315/1‑100 (MQ=255)
                                                                     ccACGATATTCATTTTGCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTc                                >  2:3056240/1‑100 (MQ=255)
                                                                     ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTc                                >  1:2144166/1‑100 (MQ=255)
                                                                      cACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                               >  1:3119664/1‑100 (MQ=255)
                                                                      cACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                               >  2:267850/1‑100 (MQ=255)
                                                                       aCGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcgc                              <  1:1752112/100‑1 (MQ=255)
                                                                        ctatatTCACTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcgag                             <  2:1744777/98‑3 (MQ=255)
                                                                        cGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcgcg                             >  1:2428146/1‑100 (MQ=255)
                                                                          atatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGt                           <  2:1265119/100‑1 (MQ=255)
                                                                           tatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGtt                          >  2:1965560/1‑100 (MQ=255)
                                                                            atTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTGCGCTCTAAATCCAACCTCTTCTTTTAGCTCGCGGTTa                         >  2:2636057/1‑100 (MQ=255)
                                                                               cATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCg                      <  2:2763321/100‑1 (MQ=255)
                                                                               cATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCg                      <  2:3420776/100‑1 (MQ=255)
                                                                                aTTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCgg                     >  1:1657724/1‑100 (MQ=255)
                                                                                  tttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCt                   >  1:1249598/1‑100 (MQ=255)
                                                                                    tGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTc                 <  1:3095698/100‑1 (MQ=255)
                                                                                     gCTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCg                <  2:481836/100‑1 (MQ=255)
                                                                                     gCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCg                >  1:2270402/1‑100 (MQ=255)
                                                                                      cTGGAAAAGTAAGGCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTAGCGGTTAGCGGCTTCGt               >  1:1816004/1‑100 (MQ=255)
                                                                                      cTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGt               >  1:3447319/1‑100 (MQ=255)
                                                                                        ggAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAANAGTNAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAg             >  1:352736/1‑100 (MQ=255)
                                                                                        ggAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAg             >  2:3562725/1‑100 (MQ=255)
                                                                                         gAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAgg            >  2:3165466/1‑99 (MQ=255)
                                                                                         gAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGa            <  2:1407223/100‑1 (MQ=255)
                                                                                          aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGAc           <  1:3259137/100‑1 (MQ=255)
                                                                                             aGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCt        >  2:1867848/1‑100 (MQ=255)
                                                                                               tAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCttt      <  2:2544321/100‑1 (MQ=255)
                                                                                                taGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTc     <  2:214688/99‑1 (MQ=255)
                                                                                                aaGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTc     >  2:1933974/1‑100 (MQ=255)
                                                                                                  gACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAc   >  1:3600458/1‑100 (MQ=255)
                                                                                                   aCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAcc  <  2:3826369/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACC  >  NC_000913/3525699‑3525897

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: