Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,904 A→G 15.3% I96I (ATT→ATC flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,9040AG15.3% 374.4 / 59.3 165I96I (ATT→ATCflhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (58/81);  new base G (14/11);  total (72/93)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01

Only 100 of 170 total aligned reads displayed.

GGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGC  >  NC_000913/1965806‑1966003
                                                                                                  |                                                                                                    
gggTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTaa                                                                                                     >  1:3570216/1‑100 (MQ=255)
   ttGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCa                                                                                                  >  2:2704399/1‑100 (MQ=255)
      aGTTTGGAAAACTTCGCCTACAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCg                                                                                               <  1:2168614/100‑1 (MQ=255)
         ttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGgca                                                                                            <  1:3278669/100‑1 (MQ=255)
         ttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGgca                                                                                            >  2:836710/1‑100 (MQ=255)
         ttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGgca                                                                                            <  2:2490186/100‑1 (MQ=255)
         ttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGgca                                                                                            <  2:3567220/100‑1 (MQ=255)
          tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgt                                                                                           <  2:1977421/100‑3 (MQ=255)
          tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGgcag                                                                                           >  1:563173/1‑100 (MQ=255)
            gAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtgg                                                                                         <  1:229874/100‑5 (MQ=255)
            gAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtgg                                                                                         <  1:1452519/100‑5 (MQ=255)
             aaatCTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCTCGCCGCTGATTAGCGgcgtggt                                                                                        <  2:441547/100‑6 (MQ=255)
             aaaaCTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtggt                                                                                        <  1:3202070/100‑6 (MQ=255)
               agctgcGGCTGCAAGGAGTTGCCGCTAAATAGTAGCCCTCCCAGCATGACCGGAGAAAGGAGGGCCACCAGCACCACGCCGCTGATTAGCGgcgtgg                                                                                         <  1:4233088/92‑5 (MQ=255)
                aCTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAgcgc                                                                                     <  1:3348863/100‑1 (MQ=255)
                aCTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAgcgc                                                                                     >  2:462436/1‑100 (MQ=255)
                 cTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGcc                                                                                    >  1:1244317/1‑100 (MQ=255)
                    cGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGc                                                                                 >  2:1695467/1‑100 (MQ=255)
                     ggCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCa                                                                                <  2:2234707/100‑1 (MQ=255)
                       cTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATg                                                                              <  2:3163417/100‑1 (MQ=255)
                       cTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATg                                                                              <  2:4280527/100‑1 (MQ=255)
                         gCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGc                                                                            >  2:940950/1‑100 (MQ=255)
                          cAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCt                                                                           >  1:173800/1‑100 (MQ=255)
                          cAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCt                                                                           <  2:668865/100‑1 (MQ=255)
                          cAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCt                                                                           >  1:978853/1‑100 (MQ=255)
                           aaGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGcct                                                                          >  1:1268195/1‑98 (MQ=255)
                              gATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtctc                                                                       <  2:2292355/100‑1 (MQ=255)
                              gATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtctc                                                                       <  1:4272647/100‑1 (MQ=255)
                                tttGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                     <  1:3971012/100‑1 (MQ=255)
                                tttGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                     >  2:2332496/1‑100 (MQ=255)
                                tttGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                     <  2:2029690/100‑1 (MQ=255)
                                 ttGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGa                                                                    <  2:1399828/100‑1 (MQ=255)
                                 ttGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGa                                                                    <  1:3838824/100‑1 (MQ=255)
                                  tGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGAt                                                                   >  1:878249/1‑100 (MQ=255)
                                   gCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                     <  1:4255238/97‑1 (MQ=255)
                                   gCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGAtn                                                                  <  1:2601081/100‑2 (MQ=255)
                                   gCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATc                                                                  <  2:2832143/100‑1 (MQ=255)
                                    ccGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATca                                                                 >  2:2118685/1‑100 (MQ=255)
                                    ccGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATca                                                                 <  1:463936/100‑1 (MQ=255)
                                     cGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATca                                                                 <  1:1193207/99‑1 (MQ=255)
                                      gCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagc                                                               <  2:184585/100‑1 (MQ=255)
                                        tAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagcag                                                             >  1:1119647/1‑100 (MQ=255)
                                        tAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagcag                                                             <  1:2246847/100‑1 (MQ=255)
                                           aTACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaat                                                          <  2:3793943/100‑1 (MQ=255)
                                           aTACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaat                                                          >  1:2641280/1‑100 (MQ=255)
                                           aTACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaat                                                          <  1:796262/100‑1 (MQ=255)
                                            tACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaata                                                         >  2:1146830/1‑100 (MQ=255)
                                              ccAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaataa                                                        >  2:1666190/1‑99 (MQ=255)
                                               cAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATc                                                      >  1:4092568/1‑100 (MQ=255)
                                                aGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCt                                                     <  1:3159115/100‑1 (MQ=255)
                                                     tCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCg                                                >  1:1046712/1‑100 (MQ=255)
                                                     tCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCg                                                >  2:3785368/1‑100 (MQ=255)
                                                     tCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCg                                                <  1:1802769/100‑1 (MQ=255)
                                                      cccAGCATGACCAGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGc                                                   <  1:2064143/96‑1 (MQ=255)
                                                        cAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAgg                                             >  2:3230883/1‑100 (MQ=255)
                                                            aTGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCGGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                         <  2:1918233/100‑1 (MQ=255)
                                                               aCCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGAt                                      <  2:1259685/100‑1 (MQ=255)
                                                               aCCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGAt                                      >  2:2712799/1‑100 (MQ=255)
                                                               aCCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGAt                                      <  2:2533103/100‑1 (MQ=255)
                                                                 cGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTc                                    <  1:1007225/100‑1 (MQ=255)
                                                                  ggAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTc                                    <  1:3419003/99‑1 (MQ=255)
                                                                  ggAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCg                                   >  1:1594706/1‑100 (MQ=255)
                                                                     gaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCggg                                 >  1:387188/1‑99 (MQ=255)
                                                                     gaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCggg                                 <  1:3257781/99‑1 (MQ=255)
                                                                       aaTGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCa                              >  2:4226961/1‑100 (MQ=255)
                                                                       aaTGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCa                              >  1:1887593/1‑100 (MQ=255)
                                                                       aaTGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTACCCGAGGATCAGATTCGGGTCa                              >  1:2472120/1‑100 (MQ=255)
                                                                        aTGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCAt                             <  2:2787563/100‑1 (MQ=255)
                                                                         tGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCAt                             >  1:1646629/1‑99 (MQ=255)
                                                                           aGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGa                          <  2:3880009/100‑1 (MQ=255)
                                                                            gcgcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGat                         <  2:3557524/100‑1 (MQ=255)
                                                                              gcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGataa                       <  1:1695461/100‑1 (MQ=255)
                                                                              gcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGataa                       >  2:75546/1‑100 (MQ=255)
                                                                                  ccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACt                   >  2:1576395/1‑100 (MQ=255)
                                                                                    gcgctgctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACtgt                 >  1:124119/9‑100 (MQ=255)
                                                                                     gCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACtgtg                >  2:454222/1‑100 (MQ=255)
                                                                                       ctgctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGAt              >  2:3397376/6‑100 (MQ=255)
                                                                                       accacGCCGCTAATCAGCGGCAGCATCGCCAGCATGGCTTCTCTGAGCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGAt              <  1:2540603/100‑1 (MQ=255)
                                                                                        ccacGCCGCTCATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAAGATCAGATTCGGGTCATTGATAATACTGTGATc             <  1:2893705/100‑1 (MQ=255)
                                                                                        ccacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGa               >  1:1519023/1‑98 (MQ=255)
                                                                                        ccacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATc             >  1:862364/1‑100 (MQ=255)
                                                                                         gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCCTTGATAGTACTGTGATCa           >  2:3055264/4‑100 (MQ=255)
                                                                                         gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATc             >  1:2378352/4‑99 (MQ=255)
                                                                                         gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa           >  2:2916979/4‑100 (MQ=255)
                                                                                         gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa           >  1:2463089/4‑100 (MQ=255)
                                                                                         gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa           >  2:219418/4‑100 (MQ=255)
                                                                                         gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa           >  1:1433115/4‑100 (MQ=255)
                                                                                         gctgccgcTAATTAGCGGCAGCAGCACCAGCATGGCTTCTCTGATCAACAGAATAATCTGCCAGAGGATCAGATTCAGGTCATTGATAATACTGTGAt              >  1:1279692/4‑98 (MQ=255)
                                                                                         cacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa           <  1:4050713/100‑1 (MQ=255)
                                                                                            gccgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCaaa         >  1:3639119/1‑99 (MQ=255)
                                                                                             ccgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAAt       <  2:4135441/100‑1 (MQ=255)
                                                                                              cgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATg      <  1:2043241/100‑1 (MQ=255)
                                                                                              cgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATg      <  2:3205611/100‑1 (MQ=255)
                                                                                               gcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc     >  2:3916657/1‑100 (MQ=255)
                                                                                               gcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc     <  2:3202867/100‑1 (MQ=255)
                                                                                               gcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc     <  2:1419678/100‑1 (MQ=255)
                                                                                               gcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc     <  1:700133/100‑1 (MQ=255)
                                                                                               gcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc     <  2:1119651/100‑1 (MQ=255)
                                                                                               gcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGACCGAGGATCAGATGCGGGTCATTGATAATACTGTGATCAAAATGc     >  1:3590658/1‑100 (MQ=255)
                                                                                                  aaTCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAAATGCAg   >  2:2084724/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
GGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGC  >  NC_000913/1965806‑1966003

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: