Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,017,469 G→A T354T (ACG→ACA)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,4690GA99.6% 867.8 / NA 242T354T (ACG→ACA) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (116/125);  total (116/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCC  >  NC_000913/2017379‑2017558
                                                                                          |                                                                                         
cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACt                                                                                         >  1:125419/1‑93 (MQ=255)
cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACt                                                                                         >  1:4728980/1‑93 (MQ=255)
cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACt                                                                                         <  2:3027204/93‑1 (MQ=255)
cTTTAGTGGTCAGAAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACt                                                                                         <  2:6048094/93‑1 (MQ=255)
  ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTc                                                                                       <  2:6559637/93‑1 (MQ=255)
    aGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTcct                                                                                     <  1:5886554/93‑2 (MQ=255)
     gTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCgg                                                                                    <  2:5757830/93‑1 (MQ=255)
     gTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCgg                                                                                    >  2:4900411/1‑93 (MQ=255)
      tGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGt                                                                                   >  1:3021575/1‑93 (MQ=255)
        gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTc                                                                                 >  1:4837683/1‑93 (MQ=255)
         tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTcc                                                                                >  1:2609990/1‑93 (MQ=255)
          cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTccc                                                                               >  1:4187170/1‑93 (MQ=255)
          cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTccc                                                                               <  1:921451/93‑1 (MQ=255)
            gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGCAGACGACGATACACTTCCGGCTCCCGt                                                                             >  2:6180459/1‑93 (MQ=255)
            gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGt                                                                             >  1:6500874/1‑93 (MQ=255)
             cagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtc                                                                            <  2:727364/93‑1 (MQ=255)
               gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                          <  2:2084629/93‑1 (MQ=255)
               gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                          >  2:946078/1‑93 (MQ=255)
                cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTctct                                                                         <  1:4434427/93‑1 (MQ=255)
                cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTctct                                                                         <  2:167027/93‑1 (MQ=255)
                cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTctct                                                                         <  1:6704518/93‑1 (MQ=255)
                 agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCtt                                                                        <  1:5379970/93‑1 (MQ=255)
                  ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGCCGGCGATACACTTCCGGATCCCGTCTCttt                                                                       >  1:3891036/1‑93 (MQ=255)
                  ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCttt                                                                       >  1:1307985/1‑93 (MQ=255)
                   gCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTa                                                                      >  2:2832291/1‑93 (MQ=255)
                    ccGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTAc                                                                     <  1:4611803/93‑1 (MQ=255)
                     cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACa                                                                    >  2:596918/1‑93 (MQ=255)
                      gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACaa                                                                   >  2:468580/1‑93 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACTATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                  <  1:4945404/93‑1 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                  >  2:6387132/1‑93 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                  <  1:691870/93‑1 (MQ=255)
                        ttCCCAGCAACAGCAAAGCCTACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAgg                                                                 <  2:4437079/93‑1 (MQ=255)
                        ttCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAgg                                                                 >  2:1635797/1‑93 (MQ=255)
                           ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCg                                                              <  1:1786959/93‑1 (MQ=255)
                             agcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtg                                                            <  2:2175808/93‑1 (MQ=255)
                               caacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTa                                                          >  2:5244513/1‑93 (MQ=255)
                                   agcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAg                                                      >  2:872462/1‑93 (MQ=255)
                                   agcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAg                                                      >  2:5540387/1‑93 (MQ=255)
                                   agcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAg                                                      <  1:5060663/93‑1 (MQ=255)
                                    gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAgg                                                     <  2:645531/93‑1 (MQ=255)
                                    gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAAACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAgg                                                     >  1:6625878/1‑93 (MQ=255)
                                      aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCa                                                   <  2:30662/93‑1 (MQ=255)
                                      aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCa                                                   >  2:1977740/1‑93 (MQ=255)
                                       aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTGTTTACAAGGGCGTGTAACAGGCaa                                                  <  1:6465481/93‑1 (MQ=255)
                                       aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCaa                                                  >  2:1359019/1‑93 (MQ=255)
                                       aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCaa                                                  >  2:831301/1‑93 (MQ=255)
                                       aaGCCAACGCACAACAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCaa                                                  <  2:1152499/93‑1 (MQ=255)
                                        aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                 <  1:4865995/93‑1 (MQ=255)
                                         gCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACa                                                <  2:3326161/93‑1 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  1:5962630/93‑1 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  1:3558274/93‑1 (MQ=255)
                                             aCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcgg                                            >  2:5835524/1‑93 (MQ=255)
                                              cGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggc                                           <  2:4232887/93‑1 (MQ=255)
                                              cGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggc                                           >  1:2326325/1‑93 (MQ=255)
                                               gCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggcg                                          >  1:3179745/1‑93 (MQ=255)
                                                cacaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGt                                         <  2:1580606/93‑1 (MQ=255)
                                                 acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                        >  2:1878942/1‑93 (MQ=255)
                                                   aGCAACCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                      <  1:2823720/93‑1 (MQ=255)
                                                    gCAAACCATGAACCTCGGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                     <  2:2336922/93‑1 (MQ=255)
                                                     cAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGata                                    >  2:4751625/1‑93 (MQ=255)
                                                     cAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGata                                    <  2:35642/93‑1 (MQ=255)
                                                      aaaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGatat                                   <  2:4828838/93‑1 (MQ=255)
                                                      aaaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGatat                                   <  1:2664424/93‑1 (MQ=255)
                                                      aaaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGatat                                   >  2:5114764/1‑93 (MQ=255)
                                                        aCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                 <  2:3093536/93‑1 (MQ=255)
                                                         ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                <  2:3429521/93‑1 (MQ=255)
                                                          cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                               >  2:3146371/1‑93 (MQ=255)
                                                          cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                               <  2:3881127/93‑1 (MQ=255)
                                                            tGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                             >  1:6192413/1‑93 (MQ=255)
                                                             gAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGcc                            <  1:1890631/93‑1 (MQ=255)
                                                             gAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGcc                            >  2:2551602/1‑93 (MQ=255)
                                                              aaCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCt                           <  1:3682449/93‑1 (MQ=255)
                                                              aaCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCt                           <  2:2632290/93‑1 (MQ=255)
                                                               aCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa                          <  1:62339/93‑1 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        <  2:3539822/93‑1 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        >  2:5427838/1‑93 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        <  1:4351161/93‑1 (MQ=255)
                                                                     ggCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa                    >  2:2431189/1‑93 (MQ=255)
                                                                      gggggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCag                   <  1:2402163/91‑1 (MQ=255)
                                                                        gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                 <  2:1536889/93‑1 (MQ=255)
                                                                         ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg                >  1:380704/1‑93 (MQ=255)
                                                                         ggggAAGACGACGATACACTTCCGGTTCCCGCCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAGCGTCAGAgg                >  2:3292785/1‑93 (MQ=255)
                                                                          gggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGt               >  2:69481/1‑93 (MQ=255)
                                                                          gggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGt               <  2:3804342/93‑1 (MQ=255)
                                                                            gAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAg             <  1:5680539/93‑1 (MQ=255)
                                                                             aaGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGc            <  2:2866743/93‑1 (MQ=255)
                                                                              aGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCa           >  1:2875329/1‑93 (MQ=255)
                                                                              aGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCa           >  1:1625745/1‑93 (MQ=255)
                                                                               gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc          >  2:1049696/1‑93 (MQ=255)
                                                                               gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc          >  2:5901378/1‑93 (MQ=255)
                                                                                  gacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGt       >  2:6598016/1‑93 (MQ=255)
                                                                                  gacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGt       >  1:1516823/1‑93 (MQ=255)
                                                                                   acgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTa      >  1:5149486/1‑93 (MQ=255)
                                                                                   acgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCCAACGTCAGAGGTAGCACCGTa      <  1:3101866/93‑1 (MQ=255)
                                                                                    cgaTACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa     <  2:4377943/93‑1 (MQ=255)
                                                                                    cgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa     >  2:4294377/1‑93 (MQ=255)
                                                                                      aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAAta   >  2:6434383/1‑92 (MQ=255)
                                                                                      aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc   <  2:6292740/93‑1 (MQ=255)
                                                                                      aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc   <  2:1810581/93‑1 (MQ=255)
                                                                                       tACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc  >  2:4106330/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCC  >  NC_000913/2017379‑2017558

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