Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,018,590 G→A V69I (GTT→ATT)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,018,5900GA97.9% 324.1 / ‑2.2 93V69I (GTT→ATT) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (45/46);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCA  >  NC_000913/2018491‑2018683
                                                                                                   |                                                                                             
gCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTg                                                                                               <  2:1452393/100‑1 (MQ=255)
   ttAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTAttt                                                                                            >  1:2299146/1‑100 (MQ=255)
    tAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTc                                                                                           <  1:1137138/100‑1 (MQ=255)
     aGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTcc                                                                                          >  1:1057750/1‑100 (MQ=255)
     aGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTcc                                                                                          <  2:1921714/100‑1 (MQ=255)
     aGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTcc                                                                                          >  2:2099224/1‑100 (MQ=255)
      gCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCg                                                                                         >  1:444477/1‑100 (MQ=255)
      gCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCg                                                                                         >  2:2229850/1‑100 (MQ=255)
       cAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGa                                                                                        >  2:1002921/1‑100 (MQ=255)
         ggCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGAcg                                                                                      <  2:596457/100‑1 (MQ=255)
            gagaCAACTACCGACAAAGCGGCCGCCCAGTTGCGGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACgcgc                                                                                   <  1:166293/100‑1 (MQ=255)
             agaCAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGca                                                                                  <  1:1121146/100‑1 (MQ=255)
             agaCAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGca                                                                                  <  2:292854/100‑1 (MQ=255)
                cAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGcaaca                                                                               >  2:1357521/1‑100 (MQ=255)
                   cTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGc                                                                            <  2:2865726/100‑1 (MQ=255)
                      ccGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTaa                                                                         >  2:2683907/1‑100 (MQ=255)
                      ccGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTaa                                                                         >  2:303577/1‑100 (MQ=255)
                       cGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAAt                                                                        >  2:930587/1‑100 (MQ=255)
                       cGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAAt                                                                        >  1:1206781/1‑100 (MQ=255)
                        gACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAAtt                                                                       <  2:154876/100‑1 (MQ=255)
                        gACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAAtt                                                                       >  2:2036685/1‑100 (MQ=255)
                           aaaGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTAc                                                                    <  2:119467/100‑1 (MQ=255)
                            aaGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACt                                                                   <  1:205785/100‑1 (MQ=255)
                            aaGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACt                                                                   <  2:1906527/100‑1 (MQ=255)
                            aaGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACt                                                                   <  1:1536823/100‑1 (MQ=255)
                            aaGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACt                                                                   >  2:2547345/1‑100 (MQ=255)
                            aaGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACt                                                                   <  1:1373510/100‑1 (MQ=255)
                             aGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTg                                                                  <  2:569113/100‑1 (MQ=255)
                              gCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGa                                                                 >  2:1601093/1‑100 (MQ=255)
                                ggCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATc                                                               <  1:773209/100‑1 (MQ=255)
                                  ccccccAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATccc                                                             >  1:1272577/1‑100 (MQ=255)
                                     cccAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCtgt                                                          <  1:887455/100‑1 (MQ=255)
                                      ccAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAGGCTAATTTACTGATCCCtgtg                                                         <  2:692772/100‑1 (MQ=255)
                                       tcttccgatctggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTgg                                                        <  1:1978553/90‑1 (MQ=255)
                                       cAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTgg                                                        >  2:556136/1‑100 (MQ=255)
                                       cAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTgg                                                        <  1:542107/100‑1 (MQ=255)
                                       cAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTgg                                                        <  1:543637/100‑1 (MQ=255)
                                         gTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGAt                                                      <  1:2546437/100‑1 (MQ=255)
                                            gCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGa                                                    <  1:2038846/99‑1 (MQ=255)
                                              tggtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAac                                                 <  1:2683883/100‑1 (MQ=255)
                                              tggtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAac                                                 <  2:611427/100‑1 (MQ=255)
                                               ggtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAaca                                                >  2:95957/1‑100 (MQ=255)
                                               ggtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAaca                                                >  1:2069125/1‑100 (MQ=255)
                                                gtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAacac                                               >  1:105529/1‑100 (MQ=255)
                                                gtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAacac                                               <  1:303576/100‑1 (MQ=255)
                                                gtggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAacac                                               >  1:315540/1‑100 (MQ=255)
                                                  ggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCg                                             <  1:2544965/100‑1 (MQ=255)
                                                  ggCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGagatcggaaga                                             >  2:1978570/1‑90 (MQ=255)
                                                   gCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGACACACCGc                                            >  2:2822336/1‑100 (MQ=255)
                                                   gCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGc                                            >  2:2865712/1‑100 (MQ=255)
                                                   gCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGc                                            >  2:1736145/1‑100 (MQ=255)
                                                    cAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGcc                                           <  1:365517/100‑1 (MQ=255)
                                                    cAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGcc                                           <  2:2812599/100‑1 (MQ=255)
                                                    cAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGcc                                           >  2:2112346/1‑100 (MQ=255)
                                                    cAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAACCTAATTTACTGATCCCTGTGGATGAAACACCGcc                                           >  1:1351441/1‑100 (MQ=255)
                                                        aGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtc                                       <  2:2299164/100‑1 (MQ=255)
                                                        aGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtc                                       <  2:1109743/100‑1 (MQ=255)
                                                         gATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtca                                      <  1:118006/100‑1 (MQ=255)
                                                          aTAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACCAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtcat                                     <  1:829402/100‑1 (MQ=255)
                                                           tAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtcatc                                    <  2:551085/100‑1 (MQ=255)
                                                           tAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtcatc                                    >  1:1146938/1‑100 (MQ=255)
                                                             aGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCaa                                  <  2:1061315/100‑1 (MQ=255)
                                                              gCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAac                                 <  2:2462418/100‑1 (MQ=255)
                                                               cccACGACAAAAGGCGAGCCGCTGATCAGCAATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAac                                 >  1:1895876/1‑99 (MQ=255)
                                                                 cACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAacgac                              <  1:134854/100‑1 (MQ=255)
                                                                 cACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAacgac                              <  2:1032941/100‑1 (MQ=255)
                                                                 cACGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAacgac                              >  1:2565318/1‑100 (MQ=255)
                                                                  aCGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAacgacg                             <  2:315542/100‑1 (MQ=255)
                                                                  aCGACAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAacgacg                             <  1:1601081/100‑1 (MQ=255)
                                                                     aCAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAAc                          >  2:1911782/1‑100 (MQ=255)
                                                                     aCAAAAGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAAc                          >  2:2792273/1‑100 (MQ=255)
                                                                       cttccgatctccGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACaa                        <  2:2508952/90‑1 (MQ=255)
                                                                        acaaacGAGCCGCTGAAGAGCGATATTATTTCCGACGCGGAACACGCTAATTTACTGTTCCCTCTGGACGAAACACAGCCTGTCATCAACGacaaacaa                        >  1:1322477/6‑99 (MQ=255)
                                                                          aGGCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATcc                     <  1:2229832/100‑1 (MQ=255)
                                                                           ggCGAGCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCca                    >  2:1584132/1‑100 (MQ=255)
                                                                                gCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCa               <  1:1062684/100‑1 (MQ=255)
                                                                                gCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCa               <  2:444479/100‑1 (MQ=255)
                                                                                gCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCa               >  1:2733874/1‑100 (MQ=255)
                                                                                gCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCa               >  2:1494195/1‑100 (MQ=255)
                                                                                gCCGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCa               >  1:1507336/1‑100 (MQ=255)
                                                                                 ccGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACaaagatcggaag              >  1:2508930/1‑90 (MQ=255)
                                                                                  cGCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAac             <  1:1357514/100‑1 (MQ=255)
                                                                                   gCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAaca            >  2:335320/1‑100 (MQ=255)
                                                                                   gCTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAaca            >  2:1065694/1‑100 (MQ=255)
                                                                                    cTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAacac           >  2:2528620/1‑100 (MQ=255)
                                                                                    cTGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAacac           <  2:832729/100‑1 (MQ=255)
                                                                                     tGATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACAcc          <  2:1146943/100‑1 (MQ=255)
                                                                                      gATCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCg         >  2:2492542/1‑100 (MQ=255)
                                                                                       aTCAGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGt        >  2:1904507/1‑100 (MQ=255)
                                                                                          aGCGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTaa     >  2:531769/1‑100 (MQ=255)
                                                                                            cGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAacc   <  2:1057755/100‑1 (MQ=255)
                                                                                            cGATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAacc   >  1:741330/1‑100 (MQ=255)
                                                                                             gATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAacca  >  1:1610699/1‑100 (MQ=255)
                                                                                             gATATTATTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAacca  >  1:1960638/1‑100 (MQ=255)
                                                                                                   |                                                                                             
GCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCA  >  NC_000913/2018491‑2018683

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