Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,942 (G)5→4 coding (214/783 nt) flgG → flagellar basal‑body rod protein FlgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,9380G.99.2% 1390.3 / 1.9 263coding (210/783 nt)flgGflagellar basal‑body rod protein FlgG
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (114/147);  total (116/147)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 285 total aligned reads displayed.

AGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTT  >  NC_000913/1134838‑1135034
                                                                                                    |                                                                                                
aGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACg                                                                                                  <  1:1659458/101‑1 (MQ=255)
 gcgcgcGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACgg                                                                                                 >  2:1659740/1‑101 (MQ=255)
 gcgcgcGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACgg                                                                                                 <  2:1659470/101‑1 (MQ=255)
     gcgGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACggggta                                                                                             >  1:1659746/1‑99 (MQ=255)
      cgGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACggggtac                                                                                            >  2:1659748/1‑98 (MQ=255)
       ggTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACggggtacg                                                                                           <  1:1659401/101‑5 (MQ=255)
          gtTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGccc                                                                                       <  2:1659661/101‑1 (MQ=255)
          gtTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGccc                                                                                       >  1:1659753/1‑101 (MQ=255)
          gtTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGccc                                                                                       <  1:1659443/101‑1 (MQ=255)
            ttGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCgg                                                                                     >  2:1659756/1‑101 (MQ=255)
            ttGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCgg                                                                                     <  2:1659358/101‑1 (MQ=255)
            ttGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCgg                                                                                     <  2:1659282/101‑1 (MQ=255)
             tGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGt                                                                                    <  1:1659377/101‑1 (MQ=255)
               aaGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCg                                                                                  >  1:1659763/1‑101 (MQ=255)
               aaGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCg                                                                                  >  1:1659760/1‑101 (MQ=255)
               aaGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCg                                                                                  >  2:1659761/1‑101 (MQ=255)
                aGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGc                                                                                 <  2:1659291/101‑1 (MQ=255)
                  aTCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTCCCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCa                                                                               <  1:1659140/101‑1 (MQ=255)
                  aTCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCa                                                                               >  1:1659769/1‑101 (MQ=255)
                  aTCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCa                                                                               >  2:1659766/1‑101 (MQ=255)
                      gCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGa                                                                           >  2:1659773/1‑101 (MQ=255)
                        tttATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAAc                                                                         <  1:1659484/101‑1 (MQ=255)
                        tttATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAAc                                                                         >  1:1659775/1‑101 (MQ=255)
                         ttATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACg                                                                        >  1:1659778/1‑101 (MQ=255)
                         ttATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACg                                                                        <  2:1659319/101‑1 (MQ=255)
                         ttATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACg                                                                        <  1:1659256/101‑1 (MQ=255)
                             cAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTa                                                                    >  2:1659782/1‑101 (MQ=255)
                             cAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTa                                                                    <  1:1659321/101‑1 (MQ=255)
                             cAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTa                                                                    >  1:1659783/1‑101 (MQ=255)
                               aaCCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACTACCTTACCCTCCGGATTACAAATCGGCACgggg                                                                                               >  1:1659789/1‑73 (MQ=255)
                               aaCCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTaca                                                                  >  2:1659788/1‑101 (MQ=255)
                                   aTTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGc                                                              >  2:1659794/1‑101 (MQ=255)
                                   aTTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGc                                                              >  1:1659796/1‑101 (MQ=255)
                                       gCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAgg                                                          <  2:1659559/101‑1 (MQ=255)
                                        ccAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTCCCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACg                                                                        <  1:1659333/86‑1 (MQ=255)
                                         cAGCCGGGGGCACAGTCTTCCGAACAAACCCCCTTCCCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGa                                                        <  1:1659429/101‑1 (MQ=255)
                                         cAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGa                                                        <  2:1659255/101‑1 (MQ=255)
                                         cAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGa                                                        >  2:1659800/1‑101 (MQ=255)
                                         cAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGa                                                        <  1:1659464/101‑1 (MQ=255)
                                          aGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGaa                                                       >  1:1659802/1‑101 (MQ=255)
                                            ccGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAAc                                                     <  1:1659574/101‑1 (MQ=255)
                                             cGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACgggg                                                                                               >  2:1659808/1‑59 (MQ=255)
                                             cGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAAcc                                                    >  2:1659805/1‑101 (MQ=255)
                                                 ggCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTc                                                <  1:1659420/101‑1 (MQ=255)
                                                 ggCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTc                                                <  2:1659197/101‑1 (MQ=255)
                                                   cacaGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGc                                              <  1:1659220/101‑1 (MQ=255)
                                                    acaGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCa                                             >  2:1659813/1‑101 (MQ=255)
                                                    acaGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCa                                             <  1:1659265/101‑1 (MQ=255)
                                                    acaGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCa                                             <  1:1659308/101‑1 (MQ=255)
                                                    acaGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCa                                             >  1:1659812/1‑101 (MQ=255)
                                                    acaGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCa                                             <  1:1659496/101‑1 (MQ=255)
                                                      aGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGa                                           <  2:1659445/101‑1 (MQ=255)
                                                       gTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGa                                                                           >  2:1659818/1‑68 (MQ=255)
                                                        tCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGAcc                                         >  1:1659821/1‑101 (MQ=255)
                                                         cTTCCGAACAAACCCCCTTGCCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACca                                        <  2:1659528/101‑1 (MQ=255)
                                                          ttCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaa                                       >  1:1659824/1‑101 (MQ=255)
                                                          ttCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaa                                       >  2:1659823/1‑101 (MQ=255)
                                                           tCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaac                                      <  1:1659190/101‑1 (MQ=255)
                                                           tCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaac                                      <  1:1659622/101‑1 (MQ=255)
                                                            ccGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaac                                      <  2:1659433/100‑1 (MQ=255)
                                                             cGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACGCAGCCAGGGAAACCTGTCGCAGACcaacaa                                    <  2:1659573/101‑1 (MQ=255)
                                                                aCAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAg                                 >  2:1659829/1‑101 (MQ=255)
                                                                 cAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGc                                >  1:1659830/1‑101 (MQ=255)
                                                                 cAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGc                                >  2:1659834/1‑101 (MQ=255)
                                                                  aaaCCACCTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCa                               >  2:1659836/1‑101 (MQ=255)
                                                                    accaccTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCaaa                             >  1:1659838/1‑101 (MQ=255)
                                                                     ccaccTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAg                            <  2:1659234/101‑1 (MQ=255)
                                                                      caccTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGa                           >  2:1659845/1‑101 (MQ=255)
                                                                      caccTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGa                           >  1:1659843/1‑101 (MQ=255)
                                                                       accTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGAt                          <  2:1659174/101‑1 (MQ=255)
                                                                       accTTACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGAt                          >  1:1659847/1‑101 (MQ=255)
                                                                           tACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTcg                      <  2:1659330/101‑1 (MQ=255)
                                                                           tACCCTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTcg                      <  2:1659603/101‑1 (MQ=255)
                                                                             cccTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTcgcg                    >  1:1659852/1‑101 (MQ=255)
                                                                             cccTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTcg                      <  2:1659346/99‑1 (MQ=255)
                                                                              ccTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCaa                              >  1:1659854/1‑90 (MQ=255)
                                                                              ccTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGAt                          <  2:1659366/94‑1 (MQ=255)
                                                                              ccTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGa                   <  2:1659456/101‑1 (MQ=255)
                                                                              ccTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGa                   <  2:1659402/101‑1 (MQ=255)
                                                                               cTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGAt                  <  2:1659393/101‑1 (MQ=255)
                                                                               cTCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGAt                  >  2:1659856/1‑101 (MQ=255)
                                                                                tCCGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGAtt                 <  1:1659357/101‑1 (MQ=255)
                                                                                 ccGGATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTa                >  2:1659859/1‑101 (MQ=255)
                                                                                   ggATTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTaaa              >  1:1659861/1‑101 (MQ=255)
                                                                                     aTTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCa                               >  2:1659864/1‑82 (MQ=255)
                                                                                     aTTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagg            <  1:1659743/101‑1 (MQ=255)
                                                                                     aTTACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagg            <  1:1659390/101‑1 (MQ=255)
                                                                                      ttACAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAaggg           <  1:1659510/101‑1 (MQ=255)
                                                                                        aCAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggca         <  2:1659098/101‑1 (MQ=255)
                                                                                        aCAAATCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggca         <  1:1659590/101‑1 (MQ=255)
                                                                                          aaaTCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggcagg       >  1:1659869/1‑101 (MQ=255)
                                                                                          aaaTCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggcagg       <  1:1659541/101‑1 (MQ=255)
                                                                                           aaTCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggcaggg      <  2:1659314/101‑1 (MQ=255)
                                                                                           aaTCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggcaggg      >  2:1659870/1‑101 (MQ=255)
                                                                                            aTCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAagggcagggc     <  1:1659506/101‑1 (MQ=255)
                                                                                             tCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCt    <  2:1659498/101‑1 (MQ=255)
                                                                                             tCGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCt    <  2:1659409/101‑1 (MQ=255)
                                                                                              cGGCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCtt   <  1:1659674/101‑1 (MQ=255)
                                                                                               ggCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCttt  <  2:1659476/101‑1 (MQ=255)
                                                                                               ggCACGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCttt  <  1:1659281/101‑1 (MQ=255)
                                                                                                    |                                                                                                
AGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTT  >  NC_000913/1134838‑1135034

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: