Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,511,417 A→G 5.5% intergenic (+90/+51) ypeC → / ← mntH DUF2502 family putative periplasmic protein/manganese/divalent cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,511,4170AG 5.5% 272.5 / 4.3 90intergenic (+90/+51)ypeC/mntHDUF2502 family putative periplasmic protein/manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (39/46);  new base G (3/2);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GCTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTG  >  NC_000913/2511269‑2511561
                                                                                                                                                    |                                                                                                                                                
gCTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATaa                                                                                                                                                 >  2:913921/1‑150 (MQ=255)
 cTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAg                                                                                                                                                <  1:872150/150‑1 (MQ=255)
  tGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAgg                                                                                                                                               <  1:468598/150‑1 (MQ=255)
  gggCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAgg                                                                                                                                               <  1:365448/149‑1 (MQ=255)
  gggCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAgg                                                                                                                                               <  2:1552114/149‑1 (MQ=255)
      gcgATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGtt                                                                                                                                           >  2:881235/1‑150 (MQ=255)
      gcgAGCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGtt                                                                                                                                           <  1:1244272/150‑1 (MQ=255)
         aTCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTACCGCGTAGGGCGGATAAGGCGTTTAc                                                                                                                                        >  2:1338432/1‑150 (MQ=255)
         aTCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAc                                                                                                                                        <  1:913921/150‑1 (MQ=255)
         aTCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAc                                                                                                                                        <  1:851769/150‑1 (MQ=255)
          tCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATCGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcg                                                                                                                                       <  1:1383010/150‑1 (MQ=255)
            gTGACTATCCTCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCTTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgcc                                                                                                                                     <  1:1477190/150‑1 (MQ=255)
            gTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTAcgcc                                                                                                                                     <  2:17546/150‑1 (MQ=255)
              gACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTATTATCAGGCCTACAAACCGGCTCTTTAACGCCTAGGGCGGATAAGGCGTTTAcgacgc                                                                                                                                   >  1:1516981/1‑150 (MQ=255)
                  atcatcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATGCGGCGTAAACGCCTTATcc                                                                                                                               <  2:23721/150‑1 (MQ=39)
                  atcatcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATcc                                                                                                                               <  1:1311020/150‑1 (MQ=255)
                    catcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCCTCCGc                                                                                                                             >  2:1525276/1‑150 (MQ=255)
                    catcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGc                                                                                                                             <  2:80715/150‑1 (MQ=255)
                    aatcGCGGAAATTGCCGGGGTCATTTGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGAGGCTAGTTCGTGTTATCAGGCCTACAAACCGGCTCTTTAACGCGTATTGCTGATAATGCGTTTACGCCGCATCCGc                                                                                                                             <  2:1282914/149‑1 (MQ=255)
                     atcGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGCCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGcc                                                                                                                            >  1:296042/1‑150 (MQ=255)
                        gcgGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAAt                                                                                                                         <  2:1238431/150‑1 (MQ=255)
                          ggAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAg                                                                                                                       >  1:673493/1‑150 (MQ=255)
                            aaaaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTg                                                                                                                     <  2:95477/150‑1 (MQ=255)
                            aaaaGGCCGGGGTCATGGGCACCGCCAGTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGTCGTATGCGGCGTAAACGCCTTATCCGccctacgcg                                                                                                                     <  1:1379485/150‑8 (MQ=37)
                             aaaGGCCGGGGTCCTGGGCACCGCCATTACACGCGTCAGCCAATGGAGCAGAATGCCTGATGCGACGCTAGAGCGACTTATCAGGCCTACAAACAGTCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACTCCGCATCCGCCAATAGTGc                                                                                                                    <  1:292421/150‑1 (MQ=37)
                              aaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGcc                                                                                                                   >  1:81566/1‑150 (MQ=255)
                              aaGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGcc                                                                                                                   >  1:1204541/1‑150 (MQ=255)
                                   cGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATg                                                                                                              >  1:69799/1‑150 (MQ=255)
                                     gggTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCg                                                                                                            >  2:1511986/1‑150 (MQ=255)
                                       gTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGAc                                                                                                          <  1:1627716/150‑1 (MQ=255)
                                          aTGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCt                                                                                                       <  2:1214603/150‑1 (MQ=255)
                                          aTGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCt                                                                                                       <  2:1637828/150‑1 (MQ=255)
                                             ggCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCAt                                                                                                    >  1:491558/1‑150 (MQ=255)
                                               cACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTc                                                                                                  <  2:1103032/150‑1 (MQ=255)
                                                  cGCCATAAAACGCGGCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTATACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAAc                                                                                               <  1:307558/150‑1 (MQ=255)
                                                   gCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACt                                                                                              <  1:323345/150‑1 (MQ=255)
                                                   gCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACt                                                                                              >  1:491942/1‑150 (MQ=255)
                                                   gCCATTAAACGCGTCCGCCAAAGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTATGGCGGATAAGGCGTTTAAGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACt                                                                                              <  2:915043/150‑1 (MQ=255)
                                                         aaaCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATc                                                                                        >  1:57647/1‑150 (MQ=255)
                                                          aaCGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATcc                                                                                       <  1:1484594/150‑1 (MQ=255)
                                                                tCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAgcgc                                                                                 <  2:151278/150‑1 (MQ=255)
                                                                  cGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACCATCCCAGCGCCg                                                                               >  1:1487440/1‑150 (MQ=255)
                                                                   gCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCACAGCGCCGt                                                                              >  2:530253/1‑150 (MQ=255)
                                                                        tGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCcca                                                                         <  1:1336472/150‑1 (MQ=255)
                                                                        tGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCcca                                                                         <  2:330965/150‑1 (MQ=255)
                                                                          gAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCccacc                                                                       <  2:369124/150‑1 (MQ=255)
                                                                          gAGCACAATACATGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCccacc                                                                       >  2:636135/1‑150 (MQ=255)
                                                                           aGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCccacca                                                                      >  2:136261/1‑150 (MQ=255)
                                                                             cacaATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCACAGCGCCGTCCCCACcaac                                                                    >  2:979120/1‑150 (MQ=255)
                                                                              acaATGCCTGATGCGACGCTAGTGCGTCTTATCATGCCTACAAACCGTCTCTTTAACGCGTAGGGCGTATAAGGCGTTTACGCCGCATCCGACAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaaca                                                                   <  2:1190676/150‑1 (MQ=255)
                                                                               caATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaacaa                                                                  >  2:1714316/1‑150 (MQ=255)
                                                                               caATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACcaacaa                                                                  >  2:544266/1‑150 (MQ=255)
                                                                                 aTGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                <  2:602926/150‑1 (MQ=255)
                                                                                   gCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAg                                                              <  1:1443346/150‑1 (MQ=255)
                                                                                   gCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGAGGATAAGGCGATTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACACTCCCAGCGCCGTCCCCACCAACAACCAg                                                              >  1:535427/1‑150 (MQ=255)
                                                                                    ccTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGa                                                             <  1:1245313/150‑1 (MQ=255)
                                                                                      tGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGata                                                           >  1:1423476/1‑150 (MQ=255)
                                                                                                cTAGTGCGTCTTATCAGTCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCTTTTACGCCGCATACGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCa                                                 <  1:1402960/150‑1 (MQ=255)
                                                                                                  aGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACg                                               >  2:1264198/1‑150 (MQ=255)
                                                                                                  aGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACg                                               <  2:1204541/150‑1 (MQ=255)
                                                                                                  aGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACg                                               <  2:81566/150‑1 (MQ=255)
                                                                                                       gTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGCTATTCAGCGCCACGACCAg                                          >  1:3319/1‑150 (MQ=255)
                                                                                                       gTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGTGCCACGACCAg                                          <  1:881235/150‑1 (MQ=255)
                                                                                                        tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGc                                         >  1:1622750/1‑150 (MQ=255)
                                                                                                        tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGc                                         >  2:298440/1‑150 (MQ=255)
                                                                                                         cTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCCCCAACAACCAGATATTCAGCGCCACGCCCAGcc                                        >  2:1447416/1‑149 (MQ=255)
                                                                                                         cTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGca                                        <  2:1577863/150‑1 (MQ=255)
                                                                                                         cTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGca                                        >  2:1454156/1‑150 (MQ=255)
                                                                                                           tATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcacc                                      <  2:1355962/150‑1 (MQ=255)
                                                                                                           tATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATAATCAGCGCCACGACCAGcacc                                      <  1:258931/150‑1 (MQ=255)
                                                                                                              cAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGACGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACcaca                                   >  2:1627312/1‑150 (MQ=255)
                                                                                                               aGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACaa                                  <  2:528317/150‑1 (MQ=255)
                                                                                                                 gCCTACAAACCGGCTCTTTAACTCTTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATGCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATc                                <  2:535427/150‑1 (MQ=255)
                                                                                                                 gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATc                                >  1:583625/1‑150 (MQ=255)
                                                                                                                   cTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCACCAACCAGATATTCAGCGCCACGACCAGCACCACAATCAc                              >  2:1495704/1‑150 (MQ=255)
                                                                                                                    tACAAACCGGCTCTTTAACGCGTAGGGCTGATAAGGCGTTTACGCCGCCTCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCCCCAACAACCAGGTATTAAGCGCCCCGCCCCGCACAACAATCAcc                             >  1:53586/1‑150 (MQ=255)
                                                                                                                              ctctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtct                   <  2:583082/150‑1 (MQ=255)
                                                                                                                              ctctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACAATAGTGCCAGATGCGACGCACATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgacg                   >  1:1382324/1‑147 (MQ=255)
                                                                                                                               tctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAc                              <  1:1220472/138‑1 (MQ=255)
                                                                                                                               tctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAc                              >  2:1220472/1‑138 (MQ=255)
                                                                                                                                  ttAACGAGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGAAGCTCATTCAACTACAATCCCAGCGCAGTCCCCACCAACATCCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGttt               <  1:204585/150‑1 (MQ=255)
                                                                                                                                    aaCGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCCCCACCATCACCCAGCCTGTCTGTTTTa             >  1:63375/1‑150 (MQ=255)
                                                                                                                                    aaCGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACAACAATCACGCAGCCTGTCGGTTTTa             >  2:981888/1‑150 (MQ=255)
                                                                                                                                          attggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGCCGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTGTTACGCGtt       >  2:929504/4‑150 (MQ=255)
                                                                                                                                           ttggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACACCCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGttt      >  1:685056/3‑150 (MQ=255)
                                                                                                                                           aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGttt      >  2:165128/1‑150 (MQ=255)
                                                                                                                                           aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGAAGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCACCAACCAGATATTCAGCGCCACGACCCGCACCACCATCACCCAGCCTGTCTGATTTACGCGttt      >  1:203307/1‑150 (MQ=39)
                                                                                                                                             ggCGGATGCGGCGTAAACGCCTTATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGc    >  1:1464596/1‑150 (MQ=255)
                                                                                                                                             ggCGGATAAGGCGTTTACGGCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGc    <  2:1516981/150‑1 (MQ=255)
                                                                                                                                              gCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCCCCAACAACCAGATATTCAGCGCCACGACCAGCACCACCATCCCCCGGCCTGTCTGTTTTACGCGTTTGCt   >  1:1511130/1‑150 (MQ=255)
                                                                                                                                               cgtaTAAGGCGTTTACGCAGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTg  <  2:1382324/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCTGGCGCGATCGTGACGATCATCGCGGAAAAGGCCGGGGTCATGGGCACCGCCATTAAACGCGTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTG  >  NC_000913/2511269‑2511561

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