Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,216,157 A→C 5.2% V107G (GTT→GGT)  yehW ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,216,1570AC 5.2% 256.8 / 3.8 90V107G (GTT→GGT) yehWputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (48/37);  new base C (3/2);  total (51/39)
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.99e-01

ACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAA  >  NC_000913/2216009‑2216301
                                                                                                                                                    |                                                                                                                                                
aCCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACCCCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGaa                                                                                                                                                  >  6:49133/1‑149 (MQ=255)
  cACGACTCATTACCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGAGCAAGGCGCTAATAGATGGTTGCAGACCCAAGCCGATCACCGGACCg                                                                                                                                                >  5:13029/1‑149 (MQ=255)
   aCGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACgg                                                                                                                                               >  5:205296/1‑149 (MQ=255)
          aTTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGc                                                                                                                                        >  4:273436/1‑149 (MQ=255)
           ttCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCGCCGTTCCGGTGATCggc                                                                                                                                       <  1:424365/149‑3 (MQ=255)
              ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGCACGGCGATCGCCAg                                                                                                                                     >  5:321975/1‑148 (MQ=255)
              ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAg                                                                                                                                     >  4:60755/1‑148 (MQ=255)
              ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAg                                                                                                                                     <  5:49133/148‑1 (MQ=255)
              ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAg                                                                                                                                     <  3:38547/148‑1 (MQ=255)
              ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAg                                                                                                                                     <  6:235393/148‑1 (MQ=255)
              ccATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGAAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCCATCCCCGGAACGGCGATCGcccg                                                                                                                                     >  1:318686/1‑146 (MQ=255)
                  aCCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCGCCGTTCCGGTGATCg                                                                                                                                         >  1:211342/1‑140 (MQ=255)
                  aCCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCGCCGTTCCGGTGATCg                                                                                                                                         <  2:211342/140‑1 (MQ=255)
                    cTTTCGCAACTTCTGTCACGATGGCATCAATCGCTCCCAGACCGGCAAGTGTCGCATGCAGGACGGGCAGCACACCGTAAAGGATCAACGCGAAAATCGCTGGTTGCAGACCAAAGCCGATCACCGGACCGGCGATCGCCAGCAATGc                                                                                                                               >  1:80405/1‑148 (MQ=255)
                       tCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGa                                                                                                                             <  7:304660/147‑1 (MQ=255)
                           aaCTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCg                                                                                                                       >  1:25254/1‑149 (MQ=255)
                                cTGTCACGCTGGCTTCATACGCTCCCAGCCCGGCAAGTGTCTCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGAACAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGaa                                                                                                                    <  5:38536/147‑1 (MQ=255)
                                   tCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGTACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGt                                                                                                                <  2:122155/148‑1 (MQ=255)
                                    cACGCGGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCGGCAGGACGGGCAGCACCCCGTAAAGGATCAAGGCGAGAATCGCTGGTTGCCGACCCAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCt                                                                                                              >  5:311214/1‑149 (MQ=255)
                                     aCGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGCCCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTg                                                                                                             >  1:341553/1‑149 (MQ=255)
                                      cGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGt                                                                                                            <  2:25254/149‑1 (MQ=255)
                                        cTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTcc                                                                                                          <  2:398375/149‑1 (MQ=255)
                                        cTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTcc                                                                                                          <  6:205296/149‑1 (MQ=255)
                                           gCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAc                                                                                                       >  8:420420/1‑149 (MQ=255)
                                             aTCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcgg                                                                                                     >  4:326440/1‑149 (MQ=255)
                                             aTCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcgg                                                                                                     <  7:420420/149‑1 (MQ=255)
                                                aaTCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggcg                                                                                                   <  3:273436/148‑1 (MQ=255)
                                                 aTCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggcgg                                                                                                  >  2:123721/1‑148 (MQ=255)
                                                     cTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCaa                                                                                               >  3:354961/1‑147 (MQ=255)
                                                      tCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGCCCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCACCGGCGGCAATAg                                                                                            >  1:349295/1‑149 (MQ=255)
                                                        cgaGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCACGGCGATAATCGCTGGTTGCAGACCAAACCCGATCACCGGAACGGCGATCGCCAGCCCTGCGACGGGCGGAAAAGGCGGTCCAACGGCGGCAATAGtt                                                                                          <  5:137734/147‑1 (MQ=255)
                                                        ccAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGtt                                                                                          >  3:260508/1‑149 (MQ=255)
                                                        ccAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGCCGGCAATAGtt                                                                                          >  5:15945/1‑149 (MQ=255)
                                                        ccAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGGCTGTCCAACGGCGGCAATAGtt                                                                                          >  6:236727/1‑149 (MQ=255)
                                                          aGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTc                                                                                        <  3:60755/149‑1 (MQ=255)
                                                           gCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTcc                                                                                       <  5:372204/149‑1 (MQ=255)
                                                                  aaGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGt                                                                                                                <  3:76989/117‑1 (MQ=255)
                                                                  aaGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGt                                                                                                                >  4:76989/1‑117 (MQ=255)
                                                                    gtgtCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTggg                                                                              >  8:198378/1‑149 (MQ=255)
                                                                    gtgtCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGACGGCAATAGTTTCCACCAGTggg                                                                              <  6:15945/149‑1 (MQ=255)
                                                                     gttggCCTGCAGGACTGGCAGTACACCGTAAAGGATCAAGGCGATAATAGGTGGTTTCAGTCCAAAGCAGATCAGCTGAACGGGGATCGAGAGCACTGTGACGGGCGGAATAGTCTGTCCATCGGCGGCAATTGTTTCCACCAGTggg                                                                              <  6:321975/144‑1 (MQ=255)
                                                                         gCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggc                                                                                                    >  8:357666/1‑122 (MQ=255)
                                                                         gCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggc                                                                                                    <  7:357666/122‑1 (MQ=255)
                                                                         gCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGaaa                                                                         <  5:236727/149‑1 (MQ=255)
                                                                         gCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGaaa                                                                         >  8:53732/1‑149 (MQ=255)
                                                                          ccTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAAt                                                                        >  2:111437/1‑149 (MQ=255)
                                                                               aGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCcgc                                                                   >  2:403932/1‑149 (MQ=255)
                                                                                ggACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCcgcg                                                                  <  1:106367/149‑1 (MQ=255)
                                                                                 gACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGCCGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCgcgc                                                                 >  5:156153/1‑149 (MQ=255)
                                                                                 gACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTcc                                                                     <  4:56662/145‑1 (MQ=255)
                                                                                 gACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTcc                                                                     >  3:56662/1‑145 (MQ=255)
                                                                                    gcaggacacaCCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGcccc                                                              <  4:340745/143‑1 (MQ=255)
                                                                                       cagcaCACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACg                                                           >  6:357064/1‑149 (MQ=255)
                                                                                         gcaCACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACggg                                                         <  1:123721/149‑1 (MQ=255)
                                                                                             acCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAg                                                     >  6:137948/1‑149 (MQ=255)
                                                                                              ccGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGt                                                    >  5:177150/1‑149 (MQ=255)
                                                                                               cGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCcgc                                                                   >  3:232978/1‑133 (MQ=255)
                                                                                               cGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCcgc                                                                   <  4:232978/133‑1 (MQ=255)
                                                                                                   aaGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAg                                               <  6:311214/149‑1 (MQ=255)
                                                                                                      gATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCaa                                            >  2:200549/1‑149 (MQ=255)
                                                                                                       aTCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAAt                                           >  7:240341/1‑149 (MQ=255)
                                                                                                          aaGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGAt                                                                                                                                           <  6:94480/50‑1 (MQ=255)
                                                                                                          aaGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGAt                                                                                                                                           >  5:94480/1‑50 (MQ=255)
                                                                                                          aaGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTcc                                        <  6:177150/149‑1 (MQ=255)
                                                                                                           aGGCGATATTCGCGGGGTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCATCGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCg                                       <  2:318686/149‑1 (MQ=255)
                                                                                                            ggCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggc                                                                                                    <  1:59491/87‑1 (MQ=255)
                                                                                                            ggCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggc                                                                                                    >  2:59491/1‑87 (MQ=255)
                                                                                                             gCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGc                                     <  7:53732/149‑1 (MQ=255)
                                                                                                                   aTCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggcg                                                                                                   >  4:136457/1‑81 (MQ=255)
                                                                                                                   aTCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAAcggcg                                                                                                   <  3:136457/81‑1 (MQ=255)
                                                                                                                    tCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTg                              >  3:77801/1‑149 (MQ=255)
                                                                                                                    tCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTg                              >  7:96414/1‑149 (MQ=255)
                                                                                                                         ggTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAAt                         >  6:91387/1‑149 (MQ=255)
                                                                                                                          gTTGCAGCCCAACGCCGATCACCGGAACGGCTATCGCCATCATTGCGTCGGGCGGAAAAGTCTGGCCTACGGCGGCAAAAGGTTCCACCAGTGGGCGAAATTCCGCGCCCCAAGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAAt                         <  6:13029/148‑1 (MQ=255)
                                                                                                                             gCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGAt                      >  2:194581/1‑148 (MQ=255)
                                                                                                                              cAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCa                    <  4:354961/149‑1 (MQ=255)
                                                                                                                              cAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCa                    <  5:357064/149‑1 (MQ=255)
                                                                                                                               aGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCAc                   >  2:219972/1‑149 (MQ=255)
                                                                                                                               aGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCAc                   >  1:13934/1‑149 (MQ=255)
                                                                                                                                 aCCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCg                 >  7:23790/1‑149 (MQ=255)
                                                                                                                                   cAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGc                >  8:64182/1‑148 (MQ=255)
                                                                                                                                   cAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGc                >  1:276390/1‑148 (MQ=255)
                                                                                                                                      aGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATa                                                                                             >  6:374023/1‑68 (MQ=255)
                                                                                                                                      aGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATa                                                                                             <  5:374023/68‑1 (MQ=255)
                                                                                                                                      aGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGcaaa             >  3:332081/1‑148 (MQ=255)
                                                                                                                                          gATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCaaacaaac        >  4:308056/1‑149 (MQ=255)
                                                                                                                                            tCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTc      >  5:84364/1‑149 (MQ=255)
                                                                                                                                            tCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTc      >  5:4542/1‑149 (MQ=255)
                                                                                                                                               ccGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGaa   <  5:137948/149‑1 (MQ=255)
                                                                                                                                               ccGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGaa   <  7:198378/149‑1 (MQ=255)
                                                                                                                                                cGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGaaa  >  3:153803/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACCACGACTCATTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAA  >  NC_000913/2216009‑2216301

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