Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 267,378 A→G W210R (TGG→CGG)  yafZ ← CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,3780AG99.3% 501.9 / NA 137W210R (TGG→CGG) yafZCP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (65/71);  minor base C (1/0);  total (66/71)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

CGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGG  >  NC_000913/267230‑267523
                                                                                                                                                    |                                                                                                                                                 
cGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGgc                                                                                                                                                 <  2:1641881/151‑1 (MQ=255)
 ggTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGgcg                                                                                                                                                <  1:1880995/151‑1 (MQ=255)
 ggTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGgcg                                                                                                                                                <  1:1173929/151‑1 (MQ=255)
    tAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCggcggc                                                                                                                                              <  2:1176050/150‑1 (MQ=255)
    tAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGgcggcg                                                                                                                                             >  2:737101/1‑151 (MQ=255)
     aaGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGa                                                                                                                                            <  2:1432085/151‑1 (MQ=255)
      aGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAg                                                                                                                                           <  2:2154522/151‑1 (MQ=255)
       gTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAgg                                                                                                                                          >  1:759655/1‑151 (MQ=255)
           cACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGgag                                                                                                                                      <  1:1816159/151‑1 (MQ=255)
              gTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAgg                                                                                                                                   >  1:1710019/1‑151 (MQ=255)
                  ccGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCt                                                                                                                               >  2:147347/1‑151 (MQ=255)
                    gTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGc                                                                                                                             <  2:470217/151‑1 (MQ=255)
                    gTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGc                                                                                                                             <  2:737539/151‑1 (MQ=255)
                    gTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGc                                                                                                                             <  2:1278463/151‑1 (MQ=255)
                    gTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGc                                                                                                                             >  2:2029342/1‑151 (MQ=255)
                      cGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACAGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCg                                                                                                                            <  1:1196978/150‑1 (MQ=255)
                        aTACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGAtt                                                                                                                         <  1:1934120/151‑1 (MQ=255)
                        aTACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGAtt                                                                                                                         <  2:1007684/151‑1 (MQ=255)
                         tACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTc                                                                                                                        >  2:301661/1‑151 (MQ=255)
                          aCCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTc                                                                                                                        <  2:1811391/150‑1 (MQ=255)
                           ccGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTAAGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTc                                                                                                                        >  2:818364/1‑149 (MQ=255)
                              cgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCa                                                                                                                   <  2:936210/151‑1 (MQ=255)
                                 aCGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCg                                                                                                                >  2:1855048/1‑151 (MQ=255)
                                           tATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTc                                                                                                      >  1:1023385/1‑151 (MQ=255)
                                            aTGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATATGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTc                                                                                                      <  2:626897/150‑1 (MQ=255)
                                            aTGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCt                                                                                                     <  1:1430520/151‑1 (MQ=255)
                                                  ccGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCAcc                                                                                                >  1:2153965/1‑150 (MQ=255)
                                                     tcctccTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCaa                                                                                              >  1:1016148/1‑149 (MQ=255)
                                                       ctcctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCg                                                                                          >  2:1197100/1‑151 (MQ=255)
                                                         cctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCgg                                                                                         <  2:1975801/150‑1 (MQ=255)
                                                         cctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGt                                                                                        <  2:1439323/151‑1 (MQ=255)
                                                            ttGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAt                                                                                      <  1:1147153/150‑1 (MQ=255)
                                                              ggCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtgt                                                                                    <  2:1710655/150‑1 (MQ=255)
                                                               gCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGa                                                                                  <  1:812422/151‑1 (MQ=255)
                                                               gCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGa                                                                                  <  1:1245291/151‑1 (MQ=255)
                                                                cATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGTGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAg                                                                                 <  2:1309640/151‑1 (MQ=255)
                                                                   tACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAggg                                                                               >  2:167560/1‑150 (MQ=255)
                                                                    aCGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcg                                                                             <  1:102053/151‑1 (MQ=255)
                                                                    aCGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcg                                                                             <  2:1984696/151‑1 (MQ=255)
                                                                    aCGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcg                                                                             >  1:1715972/1‑151 (MQ=255)
                                                                    aCGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggc                                                                              <  1:551853/150‑1 (MQ=255)
                                                                    aCGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggc                                                                              <  1:149576/150‑1 (MQ=255)
                                                                     cGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg                                                                            <  1:884662/151‑1 (MQ=255)
                                                                     cGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg                                                                            <  1:1914137/151‑1 (MQ=255)
                                                                     cGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg                                                                            >  1:2014951/1‑151 (MQ=255)
                                                                      ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcggc                                                                           >  1:865542/1‑151 (MQ=255)
                                                                      ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcggc                                                                           >  2:128919/1‑151 (MQ=255)
                                                                      ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcg                                                                             >  1:1796179/1‑149 (MQ=255)
                                                                         cACTGAGCCCGCCCTTAATCAGGTGCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCTCCAAAGCGGTATGTGAGGGCGGCttt                                                                        <  2:143908/151‑1 (MQ=255)
                                                                          aCTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCtttt                                                                       <  1:110240/151‑1 (MQ=255)
                                                                          aCTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCtttt                                                                       >  2:892441/1‑151 (MQ=255)
                                                                          aCTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCtttt                                                                       <  2:299753/151‑1 (MQ=255)
                                                                           cTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTg                                                                      >  1:420668/1‑151 (MQ=255)
                                                                              aGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCa                                                                   <  1:655318/151‑1 (MQ=255)
                                                                               gcccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGcc                                                                    >  1:1246217/1‑149 (MQ=255)
                                                                                 ccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTg                                                                >  2:163635/1‑151 (MQ=255)
                                                                                      cTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCctgc                                                           >  2:1587807/1‑151 (MQ=255)
                                                                                         aaTCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCtgt                                                        >  2:1315397/1‑151 (MQ=255)
                                                                                         aaTCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCtgt                                                        <  2:578668/151‑1 (MQ=255)
                                                                                            cAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGcc                                                     >  1:1604307/1‑151 (MQ=255)
                                                                                               gTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGTTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCggg                                                  <  2:495458/151‑1 (MQ=255)
                                                                                               gTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCggg                                                  >  1:246131/1‑151 (MQ=255)
                                                                                                 tctcCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCggggg                                                >  2:1407091/1‑151 (MQ=255)
                                                                                                   tcCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTg                                              >  1:250696/1‑151 (MQ=255)
                                                                                                   tcCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTg                                              >  1:2004924/1‑151 (MQ=255)
                                                                                                    ccTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGt                                               <  1:1881563/149‑1 (MQ=255)
                                                                                                         aTACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAAc                                         >  1:523289/1‑150 (MQ=255)
                                                                                                           aCGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGc                                      <  2:901645/151‑1 (MQ=255)
                                                                                                           aCGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGc                                      <  1:147347/151‑1 (MQ=255)
                                                                                                           aCGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCTGCTGTGCCGGGGGTGGCAACAGc                                      >  1:162740/1‑150 (MQ=255)
                                                                                                             gCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCaa                                    <  1:1738323/151‑1 (MQ=255)
                                                                                                             gCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCa                                     >  2:1885466/1‑150 (MQ=255)
                                                                                                                   aCGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTg                              >  2:101056/1‑151 (MQ=255)
                                                                                                                    cGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGc                             <  2:346125/151‑1 (MQ=255)
                                                                                                                    cGTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGc                             >  2:1882041/1‑151 (MQ=255)
                                                                                                                     gTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGc                             >  2:2167993/1‑150 (MQ=255)
                                                                                                                     gTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGc                             >  2:2093804/1‑150 (MQ=255)
                                                                                                                     gTGGTCCACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCa                            <  1:887586/151‑1 (MQ=255)
                                                                                                                          ccACAGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGc                        <  2:1147456/150‑1 (MQ=255)
                                                                                                                           cacaGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCAt                      <  2:2072798/151‑1 (MQ=255)
                                                                                                                             caGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATcc                    <  1:126231/151‑1 (MQ=255)
                                                                                                                             caGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATcc                    <  1:1166424/151‑1 (MQ=255)
                                                                                                                             caGGTCATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATcc                    <  1:1675957/151‑1 (MQ=255)
                                                                                                                                  cATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGttt               <  1:562774/151‑1 (MQ=255)
                                                                                                                                  cATTGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGtt                <  1:956541/150‑1 (MQ=255)
                                                                                                                                     tGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTtctc            >  1:2130590/1‑151 (MQ=255)
                                                                                                                                     tGCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTtctc            >  2:1517912/1‑151 (MQ=255)
                                                                                                                                      gCTCTCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcc           >  2:598980/1‑151 (MQ=255)
                                                                                                                                       ctctCATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTctcct          >  1:1332733/1‑151 (MQ=255)
                                                                                                                                           cATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCtcctcc        <  1:945396/149‑1 (MQ=255)
                                                                                                                                           cATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCtcctcc        <  1:197006/149‑1 (MQ=255)
                                                                                                                                           cATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCa       <  2:1099139/150‑1 (MQ=255)
                                                                                                                                           cATCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCa       <  2:236620/150‑1 (MQ=255)
                                                                                                                                            aTCCTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCAcc     <  1:1108273/151‑1 (MQ=255)
                                                                                                                                              ccTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCAccc    <  1:909243/150‑1 (MQ=255)
                                                                                                                                              ccTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCAccc    <  2:92446/150‑1 (MQ=255)
                                                                                                                                               cTGCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCgg  >  1:49734/1‑151 (MQ=255)
                                                                                                                                                tGCCCGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGCGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCgg  >  1:708486/1‑150 (MQ=255)
                                                                                                                                                 gCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCgg  <  1:993322/149‑1 (MQ=255)
                                                                                                                                                 gCCGGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCgg  <  2:1084338/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCTTTTGCCAGTGCCTGCTGTGCCGGGGGTGGCAACAGCAACGACTGCATGGCATCCCGTTTCTCCTCCACCCGG  >  NC_000913/267230‑267523

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