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.2% 811.7 / ‑6.1 240T354T (ACG→ACA) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (128/108);  total (129/111)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 247 total aligned reads displayed.

CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCC  >  NC_000913/2017379‑2017558
                                                                                          |                                                                                         
cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACt                                                                                         >  1:7649760/1‑93 (MQ=255)
 tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGCAGACGACGATACACtt                                                                                        >  1:5233597/1‑93 (MQ=255)
  ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTc                                                                                       <  2:5792646/93‑1 (MQ=255)
  ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTc                                                                                       >  1:946885/1‑93 (MQ=255)
  ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTc                                                                                       <  1:2867154/93‑1 (MQ=255)
   tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTcc                                                                                      <  2:1603018/93‑1 (MQ=255)
    aGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCg                                                                                     >  2:1647105/1‑93 (MQ=255)
       ggTCAGCAGCAGGCCGCTGCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGtt                                                                                  <  2:6890151/93‑1 (MQ=255)
         tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTcc                                                                                >  1:2304077/1‑93 (MQ=255)
         tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTcc                                                                                >  1:3321691/1‑93 (MQ=255)
            gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGt                                                                             >  2:4612826/1‑93 (MQ=255)
               gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                          >  1:4870640/1‑93 (MQ=255)
                 agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCtt                                                                        >  2:377447/1‑93 (MQ=255)
                   gCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTa                                                                      <  2:5839071/93‑1 (MQ=255)
                     cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACa                                                                    >  2:4255883/1‑93 (MQ=255)
                     cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAAACTCTGGCGGGGGAGGACGAAGATACACTTCCGGTTCCCCTCTCTTTTCa                                                                    >  2:4951695/1‑93 (MQ=255)
                      gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACaa                                                                   <  2:1242675/93‑1 (MQ=255)
                      gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACaa                                                                   <  2:6637348/93‑1 (MQ=255)
                      gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGAGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACaa                                                                   >  2:3573906/1‑93 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                  <  1:7494877/93‑1 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGCTTCTCGTCTCTTTAAAAg                                                                  >  2:7582003/1‑93 (MQ=255)
                        ttCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTACGGTTCCCGTCTCTTTACAAgg                                                                 >  2:7349016/1‑93 (MQ=255)
                         tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                                <  1:727485/93‑1 (MQ=255)
                         tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                                <  2:4720737/93‑1 (MQ=255)
                         tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                                >  1:6938058/1‑93 (MQ=255)
                           ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCg                                                              >  1:3928052/1‑93 (MQ=255)
                           ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCg                                                              <  1:5572783/93‑1 (MQ=255)
                           ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCg                                                              >  2:1819858/1‑93 (MQ=255)
                            cagcagaagcaaAGCCATGGCACAGCAAACCATGAACCTCTGGAGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgt                                                             <  2:86714/93‑1 (MQ=255)
                            cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgt                                                             >  1:3546090/1‑93 (MQ=255)
                            cagcaaaagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgt                                                             <  2:5972508/93‑1 (MQ=255)
                             agcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtg                                                            >  2:253268/1‑93 (MQ=255)
                             agcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtg                                                            >  2:2130816/1‑93 (MQ=255)
                               caacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTa                                                          <  2:3002866/93‑1 (MQ=255)
                               caacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTa                                                          <  2:5225962/93‑1 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGTGCGTGTaa                                                         >  2:7269109/1‑93 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         <  2:5703941/93‑1 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         >  2:7678439/1‑93 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         <  2:5172264/93‑1 (MQ=255)
                                 acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACCCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                        >  2:6107503/1‑93 (MQ=255)
                                 acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                        >  2:618774/1‑93 (MQ=255)
                                 acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                        <  1:3222113/93‑1 (MQ=255)
                                 acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                        <  2:1760246/93‑1 (MQ=255)
                                 acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                        >  2:5350812/1‑93 (MQ=255)
                                  cagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACa                                                       <  2:7201577/93‑1 (MQ=255)
                                   agcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAg                                                      <  2:2300095/93‑1 (MQ=255)
                                   agcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAg                                                      >  1:4654853/1‑93 (MQ=255)
                                      aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCa                                                   <  2:692611/93‑1 (MQ=255)
                                       aaGCCACCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCaa                                                  <  2:7836302/93‑1 (MQ=255)
                                       aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCaa                                                  >  1:4264606/1‑93 (MQ=255)
                                        aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                 <  2:5687387/93‑1 (MQ=255)
                                         gCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACa                                                >  2:2011778/1‑93 (MQ=255)
                                          ccAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTTTTTACAAGGGCGTGTAACAGGCAACAg                                               >  2:222908/1‑93 (MQ=255)
                                           cAACGCACAGCAAACCATGACCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  2:1300617/93‑1 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              >  1:1948873/1‑93 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  2:1438951/93‑1 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              >  1:3739280/1‑93 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  2:6586588/93‑1 (MQ=255)
                                             aCGAACAGCAAACCATGAACCTCTGGCTGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcgg                                            <  2:1612635/93‑1 (MQ=255)
                                              cGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggc                                           <  1:2793593/93‑1 (MQ=255)
                                              cGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggc                                           >  1:7433638/1‑93 (MQ=255)
                                                 acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                        >  1:4301816/1‑93 (MQ=255)
                                                   aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                      <  1:1233985/93‑1 (MQ=255)
                                                    gCAACCCATGAACCTCT‑GCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGata                                    <  2:709914/93‑1 (MQ=255)
                                                    gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                     <  1:5114051/93‑1 (MQ=255)
                                                       aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCGCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                  >  2:4677881/1‑93 (MQ=255)
                                                       aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                  >  1:7357803/1‑93 (MQ=255)
                                                        aCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                 >  2:6990466/1‑93 (MQ=255)
                                                        aCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                 >  1:3484740/1‑93 (MQ=255)
                                                          cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                               <  1:7707043/93‑1 (MQ=255)
                                                          cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                               <  1:841475/93‑1 (MQ=255)
                                                           aTGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCg                              >  1:4339881/1‑93 (MQ=255)
                                                            tGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                             >  1:7812401/1‑93 (MQ=255)
                                                             gAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGcc                            >  2:2617368/1‑93 (MQ=255)
                                                               cccTCTGCCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa                          >  2:4798415/2‑93 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        <  1:889147/93‑1 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        >  1:63214/1‑93 (MQ=255)
                                                                  tctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTacat                       >  1:5123745/1‑90 (MQ=255)
                                                                  tctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACg                       <  1:6968589/93‑1 (MQ=255)
                                                                    cggCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTc                     <  1:6211802/92‑1 (MQ=255)
                                                                     ggCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa                    <  2:5582474/93‑1 (MQ=255)
                                                                      gCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCag                   >  1:4457186/1‑93 (MQ=255)
                                                                       cGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga                  <  2:7460519/93‑1 (MQ=255)
                                                                          gggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATCTTTTCGCCTAACGTCAGAGGt               <  1:1409558/93‑1 (MQ=255)
                                                                          gggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGt               >  2:2875756/1‑93 (MQ=255)
                                                                          gggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGt               >  2:4499915/1‑93 (MQ=255)
                                                                           ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa              >  1:5277271/1‑93 (MQ=255)
                                                                           ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa              <  1:2546084/93‑1 (MQ=255)
                                                                           ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa              >  1:1314975/1‑93 (MQ=255)
                                                                               gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc          >  1:1994570/1‑93 (MQ=255)
                                                                                  gacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGt       >  1:3680682/1‑93 (MQ=255)
                                                                                  gacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGt       <  1:2542010/93‑1 (MQ=255)
                                                                                   acgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTa      >  1:2462540/1‑93 (MQ=255)
                                                                                   acgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTa      >  2:3253134/1‑93 (MQ=255)
                                                                                    cgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa     <  1:5383150/93‑1 (MQ=255)
                                                                                    cgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa     >  1:7091062/1‑93 (MQ=255)
                                                                                    cgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa     >  2:3260488/1‑93 (MQ=255)
                                                                                     gaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTatc    >  2:3622800/1‑91 (MQ=255)
                                                                                      aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc   >  2:4073803/1‑93 (MQ=255)
                                                                                       tACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc  <  2:1021219/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCC  >  NC_000913/2017379‑2017558

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