Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,018,423 A→G 5.6% D13G (GAC→GGC)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,018,4230AG 5.6% 207.8 / 2.7 71D13G (GAC→GGC) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAG  >  NC_000913/2018274‑2018566
                                                                                                                                                     |                                                                                                                                               
ccTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCACCGCGCCGCCATGAGGAAACCTGAATGATTCCCTTAGCGCCCTTGATTACCGCCGACGTTGa                                                                                                                                                 >  1:513988/1‑150 (MQ=255)
  tGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGaca                                                                                                                                               <  1:1477493/150‑1 (MQ=255)
   ggCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATTGATGAGTTCGCCCAGCGCGCCGCCATTAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGacac                                                                                                                                              <  1:1549408/150‑1 (MQ=255)
   ggCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGacac                                                                                                                                              >  2:138776/1‑150 (MQ=255)
       gACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGAcaccacc                                                                                                                                          <  1:1071049/150‑1 (MQ=255)
       gACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTAGCCCAGCGCGCCTCCATGATTAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGAcaccacc                                                                                                                                          <  2:1839775/150‑1 (MQ=255)
             gCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTg                                                                                                                                    >  1:205843/1‑150 (MQ=255)
              cAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGc                                                                                                                                   <  1:1953784/150‑1 (MQ=255)
               aGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGcc                                                                                                                                  >  2:1458314/1‑150 (MQ=255)
                ggAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCt                                                                                                                                 >  1:92211/1‑150 (MQ=255)
                  aaCGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTgg                                                                                                                               >  1:1958751/1‑150 (MQ=255)
                   aCGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTggc                                                                                                                              <  2:406492/150‑1 (MQ=255)
                    cGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTggcg                                                                                                                             >  1:133901/1‑150 (MQ=255)
                      gCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGGCACCACCACATTGCCTggcggc                                                                                                                           >  1:1597606/1‑150 (MQ=255)
                      gCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGGCACCACCACATTGCCTggcggc                                                                                                                           >  1:1604541/1‑150 (MQ=255)
                       cAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCa                                                                                                                          <  1:1859774/150‑1 (MQ=255)
                          tCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAg                                                                                                                       <  1:47603/150‑1 (MQ=255)
                           ccACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGc                                                                                                                      <  1:780901/150‑1 (MQ=255)
                           ccACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGc                                                                                                                      <  2:92211/150‑1 (MQ=255)
                                 cggcACTGCTTGCAGAAAACCGCCTCGATGAGAAAAAGATGGATGAGTTCGCCCTGCGCGCCGCCATGAGGAATCCTGAATGATTCGGTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCTCATTGCCTGGCGGCAAAGCCTGCGa                                                                                                                <  2:629400/150‑1 (MQ=255)
                                 cggcACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAg                                                                                                                       <  1:1702212/143‑1 (MQ=255)
                                 cggcACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAg                                                                                                                       >  2:1702212/1‑143 (MQ=255)
                                      cTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAtgctg                                                                                                           >  1:554711/1‑150 (MQ=255)
                                        gCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGca                                                                                                         <  2:240266/150‑1 (MQ=255)
                                        gCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGca                                                                                                         <  2:1666500/150‑1 (MQ=255)
                                        gCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGAGCGCCGCCATGAGGAAACATGAATGATTAGCTTAGCGCCCTTGATTACCGCCCACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGAtgctgaa                                                                                                         >  1:358072/1‑148 (MQ=255)
                                             cAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGa                                                                                                    >  1:1119096/1‑150 (MQ=255)
                                              aGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGAt                                                                                                   >  1:524641/1‑150 (MQ=255)
                                                    ccGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTtctc                                                                                             >  2:599652/1‑150 (MQ=255)
                                                          cGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCgttg                                                                                       >  2:947546/1‑150 (MQ=255)
                                                               aGAAAAAGAGGGCTGAGTTCGCCCAGCGCGCCGCCAGGAGGAAACCGGAATGATGCGCTTAGCGCCCTTGATTACCGCCGACTTTGACACCTCCACATTGCCTTGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAg                                                                                  <  2:1973518/150‑1 (MQ=255)
                                                                gAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGc                                                                                 <  2:198563/150‑1 (MQ=255)
                                                                   aaaGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGaa                                                                              <  1:1325738/150‑1 (MQ=255)
                                                                    aaGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAg                                                                             <  2:843045/150‑1 (MQ=255)
                                                                         ggatgAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTa                                                                        >  2:807920/1‑150 (MQ=255)
                                                                           atgAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCACCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGc                                                                      >  1:455939/1‑150 (MQ=255)
                                                                                ttCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCg                                                                 >  2:1006562/1‑150 (MQ=255)
                                                                                  cGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCgag                                                               <  2:205843/150‑1 (MQ=255)
                                                                                   gCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCgaga                                                              >  2:913802/1‑150 (MQ=255)
                                                                                   gCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCgaga                                                              >  1:1641776/1‑150 (MQ=255)
                                                                                    cccAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGAc                                                             <  2:1595169/150‑1 (MQ=255)
                                                                                           cgccgcCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTAcc                                                      >  2:57387/1‑150 (MQ=255)
                                                                                                 cATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACaaa                                                <  1:1458314/150‑1 (MQ=255)
                                                                                                      ggAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGGCACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGc                                           <  2:1604541/150‑1 (MQ=255)
                                                                                                      ggAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGGCACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGc                                           <  2:1597606/150‑1 (MQ=255)
                                                                                                       gAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGcc                                          <  2:1958751/150‑1 (MQ=255)
                                                                                                        aaaGCTGAATGTTTCGCTTAGCGCCCTTGATTTCCGCCGAAGTTGAGACGACCACATTGCGTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCTCGTTGTTTAGAGAAGCATTAGCAGGCGAGACAACTACCGACAATGCGGccc                                         <  2:358072/150‑1 (MQ=255)
                                                                                                         aaCCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGcccc                                        >  2:517046/1‑150 (MQ=255)
                                                                                                           ccTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGcccccc                                      >  2:424626/1‑150 (MQ=255)
                                                                                                             tGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAg                                    >  2:1628658/1‑150 (MQ=255)
                                                                                                                    ttCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCtgg                             <  1:57387/150‑1 (MQ=255)
                                                                                                                     tCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCtggt                            <  1:1057038/150‑1 (MQ=255)
                                                                                                                       gCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCtggtgg                          >  1:949310/1‑150 (MQ=255)
                                                                                                                       gCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCtggtgg                          >  1:373317/1‑150 (MQ=255)
                                                                                                                       gCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCtggtgg                          >  2:923656/1‑150 (MQ=255)
                                                                                                                        cTTAGCACCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGc                         <  2:455939/150‑1 (MQ=255)
                                                                                                                           aGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAAc                      <  2:524641/150‑1 (MQ=255)
                                                                                                                            gcgcCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACa                     >  1:70237/1‑150 (MQ=255)
                                                                                                                             cgcCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAg                    <  2:652012/150‑1 (MQ=255)
                                                                                                                             cgcCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAg                    <  1:1559447/150‑1 (MQ=255)
                                                                                                                              gcCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGAc                                                             >  2:481048/1‑108 (MQ=255)
                                                                                                                              gcCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGAc                                                             <  1:481048/108‑1 (MQ=255)
                                                                                                                                 cTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATaa                >  2:1667685/1‑150 (MQ=255)
                                                                                                                                    gATTACCGCCGACGTTGACACCACCACATTCCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGcc             <  2:527989/150‑1 (MQ=255)
                                                                                                                                        aCCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCa                        >  1:772333/1‑135 (MQ=255)
                                                                                                                                        aCCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCa                        <  2:772333/135‑1 (MQ=255)
                                                                                                                                         ccgccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGa        >  1:572461/1‑150 (MQ=255)
                                                                                                                                         ccgccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGa        >  2:1590351/1‑150 (MQ=255)
                                                                                                                                          cgccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGCTGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGAc       >  2:1941871/1‑150 (MQ=255)
                                                                                                                                          cgccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATa                 >  1:1346086/1‑140 (MQ=255)
                                                                                                                                          cgccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATa                 <  2:1346086/140‑1 (MQ=255)
                                                                                                                                          cgccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGAc       <  2:1119096/150‑1 (MQ=255)
                                                                                                                                           gccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTACACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACa      >  2:398723/1‑150 (MQ=255)
                                                                                                                                            ccgACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACaa     >  1:966984/1‑150 (MQ=255)
                                                                                                                                               aCGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAg  <  2:554711/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAG  >  NC_000913/2018274‑2018566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: