Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,019,158 G→A 43.6% S258N (AGT→AAT)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,019,1580GA43.6% 60.9 / 160.0 134S258N (AGT→AAT) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base A (31/27);  total (73/61)
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.01e-01

Only 100 of 138 total aligned reads displayed.

CCGTCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCA  >  NC_000913/2019059‑2019257
                                                                                                   |                                                                                                   
ccGTCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAg                                                                                                     <  1:3674980/100‑1 (MQ=255)
  gtccctgcACCTGTTTTGAGTGTACCGCGGGGTTCTCACGAATGGCAACAAGCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTg                                                                                                   <  2:1419045/94‑1 (MQ=255)
   tCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATAGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGc                                                                                                  <  2:761936/100‑1 (MQ=255)
    cgccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCa                                                                                                 >  1:4034994/1‑100 (MQ=255)
    cgccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCa                                                                                                 >  1:1300993/1‑100 (MQ=255)
    cgccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCa                                                                                                 <  1:1328740/100‑1 (MQ=255)
    cgccgcACCTGATTTGAGTGCAACGCTGGGTTCACACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGAAGAAACAAAATGc                                                                                                  <  1:887418/99‑1 (MQ=255)
     gccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCag                                                                                                <  1:2208344/100‑1 (MQ=255)
      ccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCaga                                                                                               <  2:221360/100‑1 (MQ=255)
       cgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCagag                                                                                              <  1:1998601/100‑1 (MQ=255)
        gcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGt                                                                                             <  2:2513817/100‑1 (MQ=255)
          aCCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTg                                                                                           <  1:1155589/100‑1 (MQ=255)
            cTGTTTTGAGTGCAACGCTGGGTTCTTACGTATCGCTTCTTTCATTACGCCAGCATGTTTTGCGGATTACCCGCTAAGGGGCCCAAATTGCagttttggg                                                                                         >  1:2177294/1‑93 (MQ=255)
             tGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGt                                                                                        >  2:284120/1‑100 (MQ=255)
             tGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGt                                                                                        >  1:3081317/1‑100 (MQ=255)
             tGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGt                                                                                        >  2:287631/1‑100 (MQ=255)
              gTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTc                                                                                       <  2:3049059/100‑1 (MQ=255)
                tttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTg                                                                                     >  2:2562662/1‑100 (MQ=255)
                nttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTg                                                                                     >  2:2403146/2‑100 (MQ=255)
                 ttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGc                                                                                    >  2:627974/1‑100 (MQ=255)
                  tGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGc                                                                                    >  1:2974021/1‑99 (MQ=255)
                     gTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGGAACAAAGTGCAGAGTTGCGTCTGCAc                                                                                  >  1:4265224/1‑98 (MQ=255)
                        cACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCa                                                                             <  2:227140/100‑1 (MQ=255)
                         aCCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAg                                                                            >  2:2579256/1‑100 (MQ=255)
                          ccGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAgg                                                                           >  2:4072068/1‑100 (MQ=255)
                          ccGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTACACCCGCCACGGGCACCAAAATGCAGAGTTGCGTCTGCACCCGCAgg                                                                           >  1:2889919/1‑100 (MQ=255)
                              tGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGAttt                                                                       <  2:3724641/100‑1 (MQ=255)
                              tGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGAttt                                                                       <  1:3554065/100‑1 (MQ=255)
                               gggTCCTCACGAATGGCAACAATCAGTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGTACCCGCAGGATTTa                                                                      <  1:1003263/100‑1 (MQ=255)
                                 gTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAgg                                                                    <  2:2978681/100‑1 (MQ=255)
                                 gTTCTCACGAACGTCAACAACCATTAAGCCAGCAAATTTCGCTGTCCTCCCGCCAGGGGCAACAAAGTGCAGAGGGGCGTCTGCACCCGCAGGATTGAgg                                                                    <  1:835250/100‑1 (MQ=255)
                                  ttCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTAACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGCACTTAggg                                                                   >  1:781431/1‑99 (MQ=255)
                                   tctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTg                                                                  <  2:640384/100‑1 (MQ=255)
                                    ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGa                                                                 >  2:4003281/1‑100 (MQ=255)
                                    ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGa                                                                 >  2:1135706/1‑100 (MQ=255)
                                     tcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGaa                                                                >  2:2405947/1‑100 (MQ=255)
                                      cACGAATGGCAACAATCATTAAGGCAGCATATTTCGCTGTTCACCCGCCAGGGGCAAAAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAg                                                               >  1:2108831/1‑100 (MQ=255)
                                      cACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACGAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAg                                                               >  2:102648/1‑100 (MQ=255)
                                      cACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAg                                                               >  2:1841883/1‑100 (MQ=255)
                                       aCGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAAGGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGt                                                              <  1:2624980/100‑1 (MQ=255)
                                        cGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTg                                                             <  1:2902503/100‑1 (MQ=255)
                                        cGAATGGCAACAATCATTAAGCCAGAATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTg                                                             >  1:641229/1‑100 (MQ=255)
                                         gAATGGCAACAATCATTATGCCGGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGa                                                                 <  1:3497478/95‑1 (MQ=255)
                                         gAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGc                                                            <  2:1459132/100‑1 (MQ=255)
                                            tGGCAACCATCATTAAGCAAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCATAGTTGCGTCTGCACCCGCAGGAGTGAGGTTAAGTGCaaa                                                         >  1:1848024/1‑100 (MQ=255)
                                            tGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCaaa                                                         >  2:1111205/1‑100 (MQ=255)
                                               caacaaTCATTAATCCAGCATATTTCGCTGTTCACCCGCCATGCGCAACAAAGTGCAGAGTTACGTCTGCACCCCCAGGGTTTAGGTGAAGTGCAAAtct                                                      >  1:2859616/1‑100 (MQ=255)
                                               caacaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAtct                                                      >  1:1852470/1‑100 (MQ=255)
                                               caacaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAtc                                                       >  2:183233/1‑99 (MQ=255)
                                                 acaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTcc                                                    >  1:687263/1‑100 (MQ=255)
                                                  caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc                                                   <  2:2665541/100‑1 (MQ=255)
                                                  caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc                                                   <  2:2255406/100‑1 (MQ=255)
                                                     tgaTTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCTGGATTTAGGTGAAGTGCAAATCTCCCTc                                                 <  2:1099014/97‑1 (MQ=255)
                                                      cATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCaa                                               >  1:3148267/1‑100 (MQ=255)
                                                      cATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCaa                                               <  1:3543801/100‑1 (MQ=255)
                                                        ttAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAg                                             <  2:3471390/100‑1 (MQ=255)
                                                          aaGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTg                                           >  1:2495478/1‑100 (MQ=255)
                                                          aaGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTg                                           <  1:3039753/100‑1 (MQ=255)
                                                            gCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGga                                         >  2:2769018/1‑100 (MQ=255)
                                                              cAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgatg                                       <  1:3230315/100‑1 (MQ=255)
                                                              cAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgatg                                       <  2:3764443/100‑1 (MQ=255)
                                                               aGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgatga                                      <  1:2769007/100‑1 (MQ=255)
                                                               aGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgatga                                      <  1:845382/100‑1 (MQ=255)
                                                                   tatTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAAc                                  >  2:52293/1‑100 (MQ=255)
                                                                    atTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAAcc                                 >  2:1332687/1‑100 (MQ=255)
                                                                     tttCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATagatc                                >  1:280972/1‑96 (MQ=255)
                                                                       tCGCTGTTCAACCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGTATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAgg                              >  1:295113/1‑100 (MQ=255)
                                                                        cGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgc                             >  1:1130544/1‑100 (MQ=255)
                                                                        cGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgc                             <  1:3337357/100‑1 (MQ=255)
                                                                          cTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgcgc                           >  1:1222311/1‑100 (MQ=255)
                                                                            gTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCaa                         >  2:1685018/1‑100 (MQ=255)
                                                                            gTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAACTCTCCCTCAAAGTGGATGATAACCAGgcgc                           >  1:4050268/1‑98 (MQ=255)
                                                                            gTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGACTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACAAGgcgc                           >  1:1000030/1‑98 (MQ=255)
                                                                             ttCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCaaa                        >  2:495619/1‑100 (MQ=255)
                                                                              tCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAAt                       >  2:731665/1‑100 (MQ=255)
                                                                               cACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATc                      >  1:1746373/1‑100 (MQ=255)
                                                                                aCCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATcc                     >  1:3288344/1‑100 (MQ=255)
                                                                                 cccGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCa                    <  2:1721678/100‑1 (MQ=255)
                                                                                    gCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGa                  <  1:3826724/99‑1 (MQ=255)
                                                                                     ccAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATg                <  1:3125477/100‑1 (MQ=255)
                                                                                      cAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATgg               <  2:3148279/100‑1 (MQ=255)
                                                                                      cAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATgg               >  1:1454445/1‑100 (MQ=255)
                                                                                        ggggCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGtt             >  2:3871556/1‑100 (MQ=255)
                                                                                        ggggCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGtt             >  1:1134548/1‑100 (MQ=255)
                                                                                           gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTc           >  1:2887457/1‑99 (MQ=255)
                                                                                           gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa          <  2:3336859/100‑1 (MQ=255)
                                                                                           gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa          >  1:2110400/1‑100 (MQ=255)
                                                                                           gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa          >  1:4247785/1‑100 (MQ=255)
                                                                                           gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa          <  2:159842/100‑1 (MQ=255)
                                                                                           gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa          >  2:4152580/1‑100 (MQ=255)
                                                                                             aacaaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCAcc        <  2:1063759/100‑1 (MQ=255)
                                                                                             aacaaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCAcc        <  2:3198957/100‑1 (MQ=255)
                                                                                             aacaaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGAAGAAGTACAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCAcc        <  1:1686439/100‑1 (MQ=255)
                                                                                              acaaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGGAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCg       >  1:2334542/1‑100 (MQ=255)
                                                                                              acaaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCg       >  2:585341/1‑100 (MQ=255)
                                                                                              acaaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATATCCCTAAAACTGGCTGATAACCAGGAGCAAATCCAGATGGTTTCa          >  1:179096/1‑97 (MQ=255)
                                                                                                aaaaTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCa     <  1:1439767/100‑1 (MQ=255)
                                                                                                  aGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATc   <  1:3737002/100‑1 (MQ=255)
                                                                                                   gTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAGCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCa  <  1:4021866/100‑1 (MQ=255)
                                                                                                   gTGCAGAGTTGCGTCTGCACCAGCAGGCATTCGGTGCACTGCAAATCTCACTCCCAGCCGATGATAACCCGCCGCAAATCCAGTTGGTTTCACCGCATCa  <  1:1419941/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCGTCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCA  >  NC_000913/2019059‑2019257

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