Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,877,945 (C)6→7 intergenic (+327/+624) iap → / ← cas2 alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,877,9391.C98.9% 1431.8 / 6.6 351intergenic (+321/+630)iap/cas2alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (180/167);  total (182/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 406 total aligned reads displayed.

ACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑CCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTC  >  NC_000913/2877839‑2878034
                                                                                                      |                                                                                                 
aCGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATgcg                                                                                                    >  2:5516763/1‑101 (MQ=255)
aCGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATgcg                                                                                                    >  1:5516765/1‑101 (MQ=255)
   cgGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑ccc                                                                                               >  2:5516775/1‑101 (MQ=255)
    ggggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑cccc                                                                                              <  1:5516041/101‑1 (MQ=255)
    ggggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑cccc                                                                                              >  1:5516779/1‑101 (MQ=255)
     gggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATgcg                                                                                                    >  2:5516785/1‑96 (MQ=255)
       gAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGTG‑‑ccccccc                                                                                          <  2:5516235/101‑2 (MQ=255)
        aaCTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑cccccccg                                                                                         <  1:5516070/101‑3 (MQ=255)
         aCTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑cccccccgg                                                                                        >  1:5516791/1‑98 (MQ=255)
          cTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑cccccccggt                                                                                       <  2:5515899/101‑5 (MQ=255)
            cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTGATGCGCCCCCCCGGttt                                                                                      >  1:5516797/1‑101 (MQ=255)
            cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGttt                                                                                      >  2:5516798/1‑101 (MQ=255)
              tAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTAt                                                                                    <  2:5516076/101‑1 (MQ=255)
               aGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATc                                                                                   >  1:5516803/1‑101 (MQ=255)
                gTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATcc                                                                                  >  1:5516809/1‑101 (MQ=255)
                 tCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATccc                                                                                 >  2:5516813/1‑101 (MQ=255)
                 tCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATccc                                                                                 <  2:5516035/101‑1 (MQ=255)
                 tCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATccc                                                                                 <  2:5516157/101‑1 (MQ=255)
                  ccATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATcccc                                                                                >  1:5516818/1‑101 (MQ=255)
                      catTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTg                                                                            >  1:5516824/1‑101 (MQ=255)
                       gttCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTgg                                                                           >  1:5516831/2‑101 (MQ=255)
                        ttCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGgc                                                                          <  1:5516229/101‑1 (MQ=255)
                        ttCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGgc                                                                          <  1:5516092/101‑1 (MQ=255)
                          ccACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTgg                                                                           <  2:5516003/98‑1 (MQ=255)
                            aCCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCgg                                                                      >  2:5516835/1‑101 (MQ=255)
                            aCCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCgg                                                                      <  1:5516263/101‑1 (MQ=255)
                             ccTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑CCCCCCGGTTTATCCCCGCTGGCGCgggg                                                                    <  2:5516189/101‑1 (MQ=255)
                              cTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCgggg                                                                    >  1:5516838/1‑101 (MQ=255)
                                 tGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCTC‑‑CCCCCCGGTTTATCCCCGCTGGCGCGGGGAAc                                                                 <  2:5515951/101‑1 (MQ=255)
                                 tGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAAc                                                                 >  1:5516845/1‑101 (MQ=255)
                                  gTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAAct                                                                >  2:5516849/1‑101 (MQ=255)
                                      gAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑CCCCTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTgg                                                            >  1:5516855/1‑101 (MQ=255)
                                       aaCTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTggg                                                           <  2:5516112/101‑1 (MQ=255)
                                       aaCTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTggg                                                           <  1:5515996/101‑1 (MQ=255)
                                       aaCTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTggg                                                           >  2:5516858/1‑101 (MQ=255)
                                         cgcccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAActc                                                               <  2:5516126/93‑1 (MQ=255)
                                         cTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGggcg                                                         >  1:5516862/1‑101 (MQ=255)
                                          tCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATcccc                                                                                >  2:5516866/1‑77 (MQ=255)
                                            ccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTgg                                                                           <  1:5516335/80‑1 (MQ=255)
                                             cGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTg                                                             <  1:5516437/93‑1 (MQ=255)
                                             cGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGctt                                                     <  2:5516324/101‑1 (MQ=255)
                                              ggTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttg                                                    >  1:5516871/1‑101 (MQ=255)
                                               gTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgc                                                   >  1:5516878/1‑101 (MQ=255)
                                                 ttATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcct                                                 >  2:5516882/1‑101 (MQ=255)
                                                 ttATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcct                                                 >  1:5516883/1‑101 (MQ=255)
                                                  tATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCgggg                                                                    >  2:5516887/1‑81 (MQ=255)
                                                     ccccGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGggcgg                                                        >  2:5516893/1‑90 (MQ=255)
                                                     ccccGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGca                                             <  1:5516411/101‑1 (MQ=255)
                                                        cGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagcc                                          <  1:5516145/101‑1 (MQ=255)
                                                        cGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagcc                                          >  2:5516898/1‑101 (MQ=255)
                                                         gCTGGCGCGGGGACCTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGACCTCTGGGCGGCTTGCCTTGcagcca                                         <  1:5516241/101‑1 (MQ=39)
                                                         gCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagcca                                         >  1:5516904/1‑101 (MQ=255)
                                                         gCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagcca                                         >  2:5516903/1‑101 (MQ=255)
                                                          cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccag                                        <  2:5516608/101‑1 (MQ=255)
                                                          cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccag                                        <  2:5516312/101‑1 (MQ=255)
                                                          cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccag                                        <  2:5516265/101‑1 (MQ=255)
                                                          cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccag                                        <  2:5516386/101‑1 (MQ=255)
                                                          cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccag                                        <  1:5516207/101‑1 (MQ=255)
                                                          cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGcagccag                                        <  1:5516110/101‑1 (MQ=255)
                                                           tGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcct                                                 >  2:5516912/1‑91 (MQ=255)
                                                           tGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑cccccccg                                                                                         >  1:5516913/1‑48 (MQ=255)
                                                              cgcgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTcc                                    <  2:5516211/101‑1 (MQ=255)
                                                               gcgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATcccc                                                                                >  1:5516924/1‑56 (MQ=255)
                                                               gcgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCca                                   <  1:5516283/101‑1 (MQ=255)
                                                               gcgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCca                                   <  2:5516178/101‑1 (MQ=255)
                                                               gcgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCca                                   >  1:5516920/1‑101 (MQ=255)
                                                                cgGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCcag                                  >  2:5516925/1‑101 (MQ=255)
                                                                 ggggAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCg                                                                               <  1:5516191/55‑1 (MQ=255)
                                                                   ggAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGaa                                                                  >  1:5516932/1‑66 (MQ=255)
                                                                    gAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCagcagc                              <  2:5516249/101‑1 (MQ=255)
                                                                    gAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCagcagc                              <  1:5516309/101‑1 (MQ=255)
                                                                    gAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCagcagc                              <  2:5516434/101‑1 (MQ=255)
                                                                     aaCTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATcccc                                                                                >  2:5516934/1‑50 (MQ=255)
                                                                       cTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGt                           <  2:5516357/101‑1 (MQ=255)
                                                                        tCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGtt                          <  1:5516348/101‑1 (MQ=255)
                                                                         cccGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGttt                         >  2:5516939/1‑101 (MQ=255)
                                                                          ccGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTa                        >  1:5516941/1‑101 (MQ=255)
                                                                           cGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGgc                                                                          >  2:5516949/1‑50 (MQ=38)
                                                                           cGGGGGATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTAt                       >  2:5516945/1‑101 (MQ=255)
                                                                            gggggATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCt                                                                             >  2:5516955/1‑46 (MQ=38)
                                                                            gggagATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCca                                   <  1:5515692/88‑1 (MQ=38)
                                                                              gggATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGaa                                                                  >  1:5516958/1‑55 (MQ=25)
                                                                               ggATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATcccc                   <  1:5516377/101‑1 (MQ=255)
                                                                                gATAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCg                  <  1:5516548/101‑1 (MQ=255)
                                                                                 aTAATGTTTACGGTCATGCGCCCCCCCGGTTTATcccc                                                                                >  1:5516967/1‑38 (MQ=38)
                                                                                 aTAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGc                 <  2:5516407/101‑1 (MQ=255)
                                                                                 aTAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGc                 <  1:5516166/101‑1 (MQ=255)
                                                                                 aTAATGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGc                 >  2:5516963/1‑101 (MQ=255)
                                                                                  tAATGTTTACGGTCATGCGCCCCCCCGGTTTATcccc                                                                                >  2:5516970/1‑37 (MQ=38)
                                                                                   aaTGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTg               <  1:5516063/101‑1 (MQ=255)
                                                                                    aTGTTTACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGcttgcc                                                  >  1:5516975/1‑65 (MQ=255)
                                                                                       tttACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGgcgc           >  1:5516980/1‑101 (MQ=255)
                                                                                        ttACGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGcgcg          >  2:5516983/1‑101 (MQ=255)
                                                                                          aCGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCggg        <  2:5516287/101‑1 (MQ=255)
                                                                                          aCGGTCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCggg        >  2:5516988/1‑101 (MQ=255)
                                                                                             gTCATGCGCCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGa      <  2:5516535/101‑1 (MQ=255)
                                                                                              tCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAAc    >  1:5516992/1‑101 (MQ=255)
                                                                                              tCATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAAc    <  1:5516497/101‑1 (MQ=255)
                                                                                               cATGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACt   <  1:5515832/101‑1 (MQ=255)
                                                                                                aTGCGCCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTc  >  1:5516996/1‑101 (MQ=255)
                                                                                                      |                                                                                                 
ACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCG‑‑CCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTC  >  NC_000913/2877839‑2878034

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