Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,898,954 G→A L176L (CTG→CTA yoaD → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,9540GA98.9% 296.5 / ‑5.3 90L176L (CTG→CTAyoaDputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (44/45);  total (44/46)
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.91e-01

AGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCTTCCGCA  >  NC_000913/1898858‑1899049
                                                                                                |                                                                                               
aGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACtt                                                                                              <  1:2738494/100‑1 (MQ=255)
   cAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGGTTAACATCGACTTACTGACGGCTATGCTACTTGAg                                                                                           <  1:1986955/100‑1 (MQ=255)
   cAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAg                                                                                           <  1:369119/100‑1 (MQ=255)
   cAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAg                                                                                           >  1:928712/1‑100 (MQ=255)
   cAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAg                                                                                           <  2:2859336/100‑1 (MQ=255)
    aGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGc                                                                                          >  2:603493/1‑100 (MQ=255)
    aGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGc                                                                                          <  2:304595/100‑1 (MQ=255)
     gTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGcc                                                                                         >  1:478237/1‑100 (MQ=255)
     gTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGcc                                                                                         >  2:1635587/1‑100 (MQ=255)
      tCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCca                                                                                        >  1:2455450/1‑100 (MQ=255)
           ttttGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACt                                                                                   <  1:1952173/100‑1 (MQ=255)
            tttGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTg                                                                                  >  1:1061946/1‑100 (MQ=255)
              tGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGcc                                                                                <  2:2347968/100‑1 (MQ=255)
               gATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCg                                                                               <  2:747225/100‑1 (MQ=255)
                aTTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGGCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGc                                                                              >  1:2418336/1‑100 (MQ=255)
                 ttCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCa                                                                             >  1:1482930/1‑100 (MQ=255)
                  tCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGTCACAACTGCCGCaa                                                                            >  2:849258/1‑100 (MQ=255)
                   cAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCaaa                                                                           >  2:685781/1‑100 (MQ=255)
                   cAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATAATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCaaa                                                                           >  2:2889262/1‑100 (MQ=255)
                    aaTGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCa                                                                             <  2:2624516/97‑1 (MQ=255)
                    aaTGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAAt                                                                          <  1:607129/100‑1 (MQ=255)
                      tGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCTCAAATCa                                                                        <  1:1592910/100‑1 (MQ=255)
                      tGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCa                                                                        <  1:2396123/100‑1 (MQ=255)
                      tGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCa                                                                        >  1:2593888/1‑100 (MQ=255)
                       ggACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAg                                                                       >  2:1393102/1‑100 (MQ=255)
                       ggACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAg                                                                       >  2:814487/1‑100 (MQ=255)
                       ggACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAg                                                                       <  1:1123915/100‑1 (MQ=255)
                        gACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGt                                                                      >  2:1040029/1‑100 (MQ=255)
                         aCGCCTGGAGCGGGCAGTAGCAAGGCTAGGGTCATGGAGATGATGAACATCGACGTACTGACGGCAATGCTACTTGAGCCACAACCGCCGCAAATCAg                                                                       <  2:2959399/98‑1 (MQ=255)
                         aCGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTa                                                                     <  1:691516/100‑1 (MQ=255)
                          cGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAg                                                                    >  1:1191060/1‑100 (MQ=255)
                            ccTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAgcg                                                                  >  1:336685/1‑100 (MQ=255)
                            ccTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATCATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAgcg                                                                  >  2:1867588/1‑100 (MQ=255)
                             cTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAgcgc                                                                 <  2:794849/100‑1 (MQ=255)
                              tGCAGCGGGCAGTAGCAATGTTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAgcgc                                                                 <  1:2489761/99‑1 (MQ=255)
                              tGCAGCGGGCAGTAGCAATGCTGGGGTCGTGGATATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGcc                                                                >  1:2192518/1‑100 (MQ=255)
                               gCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACCTCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGccn                                                               <  1:1174036/100‑2 (MQ=255)
                                  gCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGcc                                                            <  2:521446/100‑1 (MQ=255)
                                  gCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGc                                                             >  2:1084388/1‑99 (MQ=255)
                                      gCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGAc                                                        >  2:2856534/1‑100 (MQ=255)
                                        aGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACgg                                                      >  1:1128749/1‑100 (MQ=255)
                                        aGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACgg                                                      <  1:512792/100‑1 (MQ=255)
                                         ggaGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGt                                                     <  2:2669386/98‑1 (MQ=255)
                                          tAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTg                                                    >  2:811629/1‑100 (MQ=255)
                                           aGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTgg                                                   >  1:1512875/1‑100 (MQ=255)
                                           aGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTgg                                                   <  2:478232/100‑1 (MQ=255)
                                           aGCAAGACTGGGGTCATGGAGATGATTAACATCGACTTCCTGTCGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTgg                                                   <  2:746448/100‑1 (MQ=255)
                                            gCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGa                                                  >  2:346793/1‑100 (MQ=255)
                                            gCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGa                                                  >  2:752651/1‑100 (MQ=255)
                                            gCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGa                                                  <  1:222549/100‑1 (MQ=255)
                                             cAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGAc                                                 <  1:637478/100‑1 (MQ=255)
                                             cAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGAc                                                 <  2:1916441/100‑1 (MQ=255)
                                                tGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACaaa                                              >  1:1674959/1‑100 (MQ=255)
                                                      ggTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCAt                                        >  1:2286249/1‑100 (MQ=255)
                                                       gTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCGGCCTGACGGTGGACAAACGGCAtt                                       >  2:2124056/1‑100 (MQ=255)
                                                       gTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCAtt                                       >  1:870745/1‑100 (MQ=255)
                                                         cATGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTg                                     >  1:1433612/1‑100 (MQ=255)
                                                          aTGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGc                                    >  1:1505125/1‑100 (MQ=255)
                                                          aTGGAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGc                                    <  2:1061936/100‑1 (MQ=255)
                                                            ggAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGctc                                  >  2:1832620/1‑100 (MQ=255)
                                                            ggAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGctc                                  <  2:336680/100‑1 (MQ=255)
                                                            ggAGATGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGctc                                  <  2:2102465/100‑1 (MQ=255)
                                                              agaTGATTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTa                                <  1:752658/100‑1 (MQ=255)
                                                               gatgatTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTAt                               <  1:2807431/100‑1 (MQ=255)
                                                                atgatTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATg                              <  1:1048898/100‑1 (MQ=255)
                                                                 tgatTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATgg                             <  2:1307586/100‑1 (MQ=255)
                                                                 tgatTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATgg                             >  2:2610390/1‑100 (MQ=255)
                                                                 tgatTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATgg                             <  1:2096002/100‑1 (MQ=255)
                                                                 tgatTAACATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATgg                             >  1:1926304/1‑100 (MQ=255)
                                                                      aaCATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATg                        >  2:2568538/1‑100 (MQ=255)
                                                                      aaCATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATg                        >  1:84475/1‑100 (MQ=255)
                                                                       aCATCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATgg                       >  1:504091/1‑100 (MQ=255)
                                                                        cATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATggg                      <  2:92332/100‑1 (MQ=255)
                                                                         aTCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGc                     <  2:2745808/100‑1 (MQ=255)
                                                                         aTCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCATGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGc                     >  2:1644919/1‑100 (MQ=255)
                                                                          tCGACTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCt                    <  2:928704/100‑1 (MQ=255)
                                                                             aCTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGt                 <  1:285676/100‑1 (MQ=255)
                                                                              cTTACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTa                <  1:597514/100‑1 (MQ=255)
                                                                                tACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGa              <  2:504897/100‑1 (MQ=255)
                                                                                tACTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGa              >  1:71751/1‑100 (MQ=255)
                                                                                 aCTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGAt             >  2:1185130/1‑100 (MQ=255)
                                                                                  cTGACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGAtt            <  1:2369360/100‑1 (MQ=255)
                                                                                    gACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTcc          >  1:1552054/1‑100 (MQ=255)
                                                                                    gACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTcc          <  2:1505113/100‑1 (MQ=255)
                                                                                    gACGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTcc          >  1:1325461/1‑100 (MQ=255)
                                                                                      cGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGCTTACCt        >  1:2894680/1‑100 (MQ=255)
                                                                                      cGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCt        <  2:2016896/100‑1 (MQ=255)
                                                                                      cGGCAATGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCAAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCt        <  2:137160/100‑1 (MQ=255)
                                                                                        gcgaTGCTACTTGAGCCACAACTGCAGCAAATCAGTTGCGCCAGCCTGGCGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCtgc      <  2:1915455/97‑3 (MQ=255)
                                                                                          aaTGCTACTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCTTCCg    <  2:1368680/100‑1 (MQ=255)
                                                                                            tgctacTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCTTCCGCa  <  1:1645413/100‑1 (MQ=255)
                                                                                            tgctacTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCTTCCGCa  <  2:1482918/100‑1 (MQ=255)
                                                                                                |                                                                                               
AGGCAGTCCGGTTTTGATTCAATGGACGCCTGCAGCGGGCAGTAGCAATGCTGGGGTCATGGAGATGATTAACATCGACTTACTGACGGCAATGCTGCTTGAGCCACAACTGCCGCAAATCAGTAGCGCCAGCCTGACGGTGGACAAACGGCATTTGCTCTATGGTAATGGGCTGGTAGATTCCCTTCCGCA  >  NC_000913/1898858‑1899049

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