Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,904 A→G 5.0% 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,9040AG 5.0% 329.7 / 4.7 99I96I (ATT→ATCflhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (43/51);  new base G (2/3);  total (45/54)
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.17e-01

GCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCA  >  NC_000913/1965756‑1966051
                                                                                                                                                    |                                                                                                                                                   
gCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTa                                                                                                                                                     >  2:100270/1‑149 (MQ=255)
  aaCTCCGCGCCAGTCTGAGCCGAGAACATCCGTTGAATGCCCGGTAGCGGGTTGAGTGTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACAGGAGAAATGAGCGCCACCAGCACCACGCCGCTAAt                                                                                                                                                   >  2:186067/1‑149 (MQ=255)
   aCTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATc                                                                                                                                                  <  1:737444/149‑1 (MQ=255)
     tCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAg                                                                                                                                                >  4:185679/1‑149 (MQ=255)
      ccGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGc                                                                                                                                               <  1:629904/149‑1 (MQ=255)
       cgcgCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCg                                                                                                                                              >  3:670323/1‑149 (MQ=255)
       cgcgCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCg                                                                                                                                              >  1:729409/1‑149 (MQ=255)
       cgcgCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCg                                                                                                                                              <  1:602646/149‑1 (MQ=255)
        gcgcCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCgg                                                                                                                                             >  4:192795/1‑149 (MQ=255)
        gcgcCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCgg                                                                                                                                             >  4:242503/1‑149 (MQ=255)
        gcgcCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCgg                                                                                                                                             >  1:796804/1‑149 (MQ=255)
           ccAGTCTGAGCCGAGAACGTCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgt                                                                                                                                          <  4:838658/149‑3 (MQ=255)
            cAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCggtgtg                                                                                                                                         <  1:771429/149‑5 (MQ=255)
                 tGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAgcgc                                                                                                                                    <  2:411918/149‑1 (MQ=255)
                  gAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGcc                                                                                                                                   >  1:490286/1‑149 (MQ=255)
                  gAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGcc                                                                                                                                   >  4:308034/1‑149 (MQ=255)
                        agaACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATg                                                                                                                             >  1:286981/1‑149 (MQ=255)
                           aCATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAATGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCt                                                                                                                          <  4:747920/149‑1 (MQ=255)
                            cATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCtt                                                                                                                         <  2:729409/149‑1 (MQ=255)
                             aTCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtc                                                                                                                        >  1:885098/1‑149 (MQ=255)
                             aTCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtc                                                                                                                        >  4:25362/1‑149 (MQ=255)
                              tCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtct                                                                                                                       >  1:59588/1‑149 (MQ=255)
                              tCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtct                                                                                                                       <  3:137374/149‑1 (MQ=255)
                                cGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTctct                                                                                                                     >  2:444474/1‑149 (MQ=255)
                                  tttAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGa                                                                                                                   >  2:446648/1‑149 (MQ=255)
                                  tttAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGa                                                                                                                   >  2:228351/1‑149 (MQ=255)
                                  tttAATGCCCGGTAGCGGGTTGAGATTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGa                                                                                                                   >  3:814138/1‑149 (MQ=255)
                                   ttAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGAt                                                                                                                  <  2:59588/149‑1 (MQ=255)
                                    tAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATc                                                                                                                 >  2:433256/1‑149 (MQ=255)
                                    tAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATc                                                                                                                 <  4:118268/149‑1 (MQ=255)
                                     aaTGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCGTGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATca                                                                                                                <  4:233552/149‑1 (MQ=255)
                                        gCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagca                                                                                                             <  3:137586/149‑1 (MQ=255)
                                        gCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagca                                                                                                             >  1:538003/1‑149 (MQ=255)
                                           cGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaa                                                                                                          <  2:762721/149‑1 (MQ=255)
                                                gCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATc                                                                                                     >  3:834248/1‑149 (MQ=255)
                                                gCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATc                                                                                                     >  1:788488/1‑149 (MQ=255)
                                                    gTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGcc                                                                                                 >  4:389899/1‑149 (MQ=255)
                                                        aGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAg                                                                                             <  4:6093/149‑1 (MQ=255)
                                                        aGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAg                                                                                             <  3:25362/149‑1 (MQ=255)
                                                          tttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagca                                                                                                             <  4:387293/131‑1 (MQ=255)
                                                          tttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagca                                                                                                             >  3:387293/1‑131 (MQ=255)
                                                          tttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGa                                                                                           >  3:291724/1‑149 (MQ=255)
                                                            tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtg                                                                                                                                         >  4:389443/1‑98 (MQ=255)
                                                            tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtg                                                                                                                                         <  3:389443/101‑4 (MQ=255)
                                                            tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATc                                                                                         <  4:834248/149‑1 (MQ=255)
                                                            tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATc                                                                                         <  1:433256/149‑1 (MQ=255)
                                                             ggAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                                                                        >  1:91853/1‑149 (MQ=255)
                                                             ggAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                                                                        <  3:325089/149‑1 (MQ=255)
                                                             ggAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATAGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                                                                        <  3:424134/149‑1 (MQ=255)
                                                                    ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTACTCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                 >  4:725089/1‑149 (MQ=255)
                                                                    ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCATCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                 <  1:578941/149‑1 (MQ=255)
                                                                    ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                 <  3:869756/149‑1 (MQ=255)
                                                                    ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                 <  3:442601/149‑1 (MQ=255)
                                                                    ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                 <  2:788786/149‑1 (MQ=255)
                                                                    ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                 <  1:498442/149‑1 (MQ=255)
                                                                      cGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGt                                                                               <  1:712896/149‑1 (MQ=255)
                                                                         cTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCAt                                                                            >  1:249727/1‑149 (MQ=255)
                                                                          tGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCAtt                                                                           <  1:722357/149‑1 (MQ=255)
                                                                           gCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTg                                                                          >  1:524767/1‑149 (MQ=255)
                                                                            cAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGa                                                                         <  4:393784/149‑1 (MQ=255)
                                                                               ggATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGa                                                                         <  3:387127/146‑1 (MQ=255)
                                                                               ggATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGAGTCGGGTCATTGa                                                                         >  4:387127/1‑146 (MQ=255)
                                                                                 aTTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGataata                                                                    <  1:333578/149‑1 (MQ=255)
                                                                                   ttGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACt                                                                  >  2:489289/1‑149 (MQ=255)
                                                                                       cGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGa                                                              <  2:91853/149‑1 (MQ=255)
                                                                                       cGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGa                                                              <  2:885098/149‑1 (MQ=255)
                                                                                       cGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGa                                                              <  1:489289/149‑1 (MQ=255)
                                                                                         cTAAATACCAGCCCTCCCAACATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATc                                                            <  3:812526/149‑1 (MQ=255)
                                                                                             aTACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCaaaa                                                        <  4:882263/149‑1 (MQ=255)
                                                                                              tACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAAt                                                       <  3:573178/149‑1 (MQ=255)
                                                                                               aCCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATg                                                      >  4:829660/1‑149 (MQ=255)
                                                                                                  aGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAg                                                   >  1:752763/1‑149 (MQ=255)
                                                                                                  aGCCCTCCCAGAATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGCTAATACTGTGATCAAAATGCAg                                                   >  3:498826/1‑149 (MQ=255)
                                                                                                       tCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAg                                              >  3:565044/1‑149 (MQ=255)
                                                                                                         ccAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCg                                            >  3:630067/1‑149 (MQ=255)
                                                                                                         ccAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCg                                            >  2:512098/1‑149 (MQ=255)
                                                                                                           aGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGga                                          <  3:156474/149‑1 (MQ=255)
                                                                                                             cATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGgaga                                        <  3:322229/149‑1 (MQ=255)
                                                                                                               tGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGc                                      <  1:739120/149‑1 (MQ=255)
                                                                                                                   cGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGc                                  <  1:889370/149‑1 (MQ=255)
                                                                                                                   cGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCGGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGc                                  <  4:498826/149‑1 (MQ=255)
                                                                                                                     gagaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGccc                                >  4:210599/1‑149 (MQ=255)
                                                                                                                     gagaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGccc                                <  1:32141/149‑1 (MQ=255)
                                                                                                                     gagaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGccc                                >  1:137819/1‑149 (MQ=255)
                                                                                                                       gaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCTCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGa                              <  2:550567/149‑1 (MQ=255)
                                                                                                                           tGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAAt                          >  4:649646/1‑149 (MQ=255)
                                                                                                                           tGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAAt                          >  2:154204/1‑149 (MQ=255)
                                                                                                                           tGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAAt                          >  1:420212/1‑149 (MQ=255)
                                                                                                                            gAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATc                         <  2:500446/149‑1 (MQ=255)
                                                                                                                            gAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATc                         <  3:210599/149‑1 (MQ=255)
                                                                                                                               cgcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGAc                      <  3:99761/149‑1 (MQ=255)
                                                                                                                               cgcCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGAAAATCGAc                      <  2:490286/149‑1 (MQ=255)
                                                                                                                                 ccaccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACgg                    <  4:670323/149‑1 (MQ=255)
                                                                                                                                 ccaccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACgg                    <  1:228351/149‑1 (MQ=255)
                                                                                                                                   accaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGc                  >  3:327971/1‑149 (MQ=255)
                                                                                                                                        caccacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGcatgct             <  1:100270/149‑6 (MQ=255)
                                                                                                                                         ctgctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGa            >  3:80937/6‑149 (MQ=255)
                                                                                                                                         accacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGa                       >  2:289543/1‑138 (MQ=255)
                                                                                                                                         accacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGa                       <  1:289543/138‑1 (MQ=255)
                                                                                                                                                   tAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCa  <  3:389899/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCA  >  NC_000913/1965756‑1966051

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: