Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,021,787 (C)5→6 coding (284/366 nt) fliO → flagellar biosynthesis protein FliO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,021,7821.C99.6% 1174.3 / ‑3.3 278coding (279/366 nt)fliOflagellar biosynthesis protein FliO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (144/133);  total (145/133)
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.93e-01

Only 100 of 311 total aligned reads displayed.

CTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGCGTCGT  >  NC_000913/2021682‑2021877
                                                                                                     |                                                                                                
cTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACtt                                                                                                   <  1:3687121/101‑1 (MQ=255)
cTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACtt                                                                                                   >  1:3687569/1‑101 (MQ=255)
 ttCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTc                                                                                                 <  2:3687042/101‑1 (MQ=255)
 ttCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTc                                                                                                 >  2:3687571/1‑101 (MQ=255)
  tCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTcc                                                                                               >  1:3687576/1‑101 (MQ=255)
   cACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGGGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTccc                                                                                              >  2:3687577/1‑101 (MQ=255)
     cTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTccccc                                                                                            <  1:3687067/101‑1 (MQ=255)
      tGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTcccccc                                                                                           <  2:3687496/101‑2 (MQ=255)
       gggCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTcccccct                                                                                          <  2:3687376/101‑3 (MQ=255)
         gcgcgcGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTccccccttc                                                                                        <  2:3687156/101‑5 (MQ=255)
         gcgcgcGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTccccccttc                                                                                        <  1:3687097/101‑5 (MQ=255)
         gcgcgcGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTccccccttc                                                                                        >  1:3687582/1‑97 (MQ=255)
          cgcgcgTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTccccccttct                                                                                       <  2:3687238/101‑6 (MQ=255)
               gTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCAcc                                                                                   >  1:3687587/1‑101 (MQ=255)
                tGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCa                                                                                  >  2:3687592/1‑101 (MQ=255)
                tGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCa                                                                                  >  1:3687590/1‑101 (MQ=255)
                  aaaGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAAc                                                                                >  1:3687595/1‑101 (MQ=255)
                    aGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACgg                                                                              <  1:3687363/101‑1 (MQ=255)
                    aGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACgg                                                                              <  1:3687233/101‑1 (MQ=255)
                     gggTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGa                                                                             <  2:3687218/101‑1 (MQ=255)
                      ggTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGaa                                                                            >  1:3687599/1‑101 (MQ=255)
                         tgtgGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAagag                                                                         <  1:3687190/101‑1 (MQ=255)
                           tggtggTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAgaga                                                                        >  2:3687602/1‑100 (MQ=255)
                              tggtCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATAcc                                                                    <  1:3686897/101‑1 (MQ=255)
                               ggtCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCg                                                                   >  1:3687607/1‑101 (MQ=255)
                                gtCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAagag                                                                         >  2:3687609/1‑94 (MQ=255)
                                  cGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAg                                                                >  1:3687613/1‑101 (MQ=255)
                                  cGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAg                                                                <  2:3687138/101‑1 (MQ=255)
                                   gATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGa                                                               <  2:3686995/101‑1 (MQ=255)
                                   gATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGa                                                               >  2:3687614/1‑101 (MQ=255)
                                      gtgGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCg                                                            >  1:3687618/1‑101 (MQ=255)
                                        ggAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGAc                                                              >  1:3687621/1‑97 (MQ=255)
                                         gAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGAtt                                                         <  1:3687258/101‑1 (MQ=255)
                                         gAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGAtt                                                         <  2:3687265/101‑1 (MQ=255)
                                         gAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGAt                                                          <  1:3687152/100‑1 (MQ=255)
                                          aaGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGAttt                                                        >  2:3687623/1‑101 (MQ=255)
                                            gATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTc                                                      >  2:3687627/1‑101 (MQ=255)
                                               gCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGt                                                   >  1:3687631/1‑101 (MQ=255)
                                               gCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGt                                                   <  2:3687087/101‑1 (MQ=255)
                                                cACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTc                                                  >  2:3687634/1‑101 (MQ=255)
                                                 aCGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCg                                                 >  1:3687636/1‑101 (MQ=255)
                                                 aCGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCg                                                 <  1:3687128/101‑1 (MQ=255)
                                                   ggCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGt                                               <  1:3687012/101‑1 (MQ=255)
                                                    gCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTc                                              >  1:3687640/1‑101 (MQ=255)
                                                        gTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGa                                          <  2:3687200/101‑1 (MQ=255)
                                                          gCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGaaa                                        <  2:3687167/101‑1 (MQ=255)
                                                           cTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGaaaa                                       >  1:3687644/1‑101 (MQ=255)
                                                             cGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAAt                                     >  2:3687646/1‑101 (MQ=255)
                                                              ggCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAAtt                                    <  1:3686972/101‑1 (MQ=255)
                                                              ggCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAAtt                                    <  2:3687279/101‑1 (MQ=255)
                                                              ggCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAAtt                                    >  1:3687648/1‑101 (MQ=255)
                                                               gCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAAttt                                   <  1:3687183/101‑1 (MQ=255)
                                                               gCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAAttt                                   >  1:3687652/1‑101 (MQ=255)
                                                                cGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTg                                  >  2:3687655/1‑101 (MQ=255)
                                                                 gTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGc                                 >  1:3687658/1‑101 (MQ=255)
                                                                   tACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCtt                               <  2:3687172/101‑1 (MQ=255)
                                                                   tACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCtt                               >  1:3687662/1‑101 (MQ=255)
                                                                    aCCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTa                              <  1:3687336/101‑1 (MQ=255)
                                                                    aCCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTa                              <  2:3687329/101‑1 (MQ=255)
                                                                     ccGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTaa                             <  2:3687340/101‑1 (MQ=255)
                                                                      cGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAg                            <  2:3687310/101‑1 (MQ=255)
                                                                      cGCAGGTCAAATCAATCTGCTGCATAAACTTCACCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTa                              >  1:3687666/1‑99 (MQ=255)
                                                                       gCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGc                           >  2:3687668/1‑101 (MQ=255)
                                                                        cAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCg                          >  1:3687672/1‑101 (MQ=255)
                                                                        cAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCg                          <  1:3687281/101‑1 (MQ=255)
                                                                        cAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCg                          >  2:3687673/1‑101 (MQ=255)
                                                                        cAGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCg                          <  1:3687394/101‑1 (MQ=255)
                                                                         aGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGt                         >  1:3687680/1‑101 (MQ=255)
                                                                         aGGTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGt                         >  2:3687679/1‑101 (MQ=255)
                                                                          ggTCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTa                        <  2:3687192/101‑1 (MQ=255)
                                                                            tCAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGc                      <  1:3687421/101‑1 (MQ=255)
                                                                             cAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGAt                                                                       <  2:3687076/51‑1 (MQ=255)
                                                                             cAAATCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCg                     >  1:3687685/1‑101 (MQ=255)
                                                                              aaaTCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCgg                    <  1:3687302/101‑1 (MQ=255)
                                                                              aaaTCAATCTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCgg                    >  2:3687688/1‑101 (MQ=255)
                                                                               aatcaatcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCg                          >  1:3687690/1‑94 (MQ=255)
                                                                                 tcaatcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGgag                 <  1:3687374/101‑1 (MQ=255)
                                                                                  caatcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGgaga                >  1:3687697/1‑101 (MQ=255)
                                                                                  caatcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGgaga                >  2:3687692/1‑101 (MQ=255)
                                                                                   aatcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGaa                                         >  2:3687699/1‑75 (MQ=255)
                                                                                    atcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATc              <  2:3687110/101‑1 (MQ=255)
                                                                                     tcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATcc             >  2:3687703/1‑101 (MQ=255)
                                                                                     tcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATcc             >  1:3687704/1‑101 (MQ=255)
                                                                                     tcTGCTGCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCAATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATcc             <  1:3687355/101‑1 (MQ=255)
                                                                                      ctgctgCATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCt            >  1:3687708/1‑101 (MQ=255)
                                                                                       tgctgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTg           <  2:3687250/101‑1 (MQ=255)
                                                                                       tgctgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTg           <  1:3687217/101‑1 (MQ=255)
                                                                                       tgctgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTg           >  1:3687711/1‑101 (MQ=255)
                                                                                         ctgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGAt         >  2:3687714/1‑101 (MQ=255)
                                                                                         ctgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGAt         <  1:3687225/101‑1 (MQ=255)
                                                                                          tgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATg        <  1:3687056/101‑1 (MQ=255)
                                                                                          tgcATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATg        >  1:3687716/1‑101 (MQ=255)
                                                                                            cATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcg      <  2:3687290/101‑1 (MQ=255)
                                                                                            cATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcg      >  1:3687720/1‑101 (MQ=255)
                                                                                            cATAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcg      <  2:3687418/101‑1 (MQ=255)
                                                                                             aTAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcgt     <  1:3687209/101‑1 (MQ=255)
                                                                                              tAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcgtc    <  1:3687397/101‑1 (MQ=255)
                                                                                              tAAACTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcgtc    <  2:3687032/101‑1 (MQ=255)
                                                                                               aaaCTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcgtcg   >  2:3687725/1‑101 (MQ=255)
                                                                                                aaCTTCCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGcgtcgt  >  1:3687728/1‑101 (MQ=255)
                                                                                                     |                                                                                                
CTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGCTGGTGCTCGGCGTTACCGCAGGTCAAATCAATCTGCTGCATAAACTTCCCCCTTCTGCACCAACGGAAGAGATACCGCAGACCGATTTTCAGTCGGTCATGAAAAATTTGCTTAAGCGTAGCGGGAGATCCTGATGCGTCGT  >  NC_000913/2021682‑2021877

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: