Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,567,815 C→T A337T (GCT→ACT)  yddW ← liprotein, glycosyl hydrolase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,567,8150CT93.0% 353.2 / ‑5.6 102A337T (GCT→ACT) yddWliprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (40/53);  minor base G (7/0);  total (48/54)
Fisher's exact test for biased strand distribution p-value = 3.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCC  >  NC_000913/1567671‑1567962
                                                                                                                                                       |                                                                                                                                                   
ctATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATg                                                                                                                                                 <  2:1010506/149‑1 (MQ=255)
  aTCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGt                                                                                                                                                <  2:437223/148‑1 (MQ=255)
  aTCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCCCCTCGGGTTTATTACAACCCCCCACCCCTTTTCCAACAACTCCTAAAGCGCGGCACTCCGTGGGGAAGGCCCCGAAAATTGTGGTGGAATGt                                                                                                                                                >  4:117980/1‑148 (MQ=255)
      ttGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTGGCGGTCGGGTTTAAGACATCCGCCCCCCACTTTTCCAAAACCGTCTAAAACGCGGCACCCCCCGAGAACAGGCCGAAAAATTTGGGGGGAAAttaat                                                                                                                                              >  4:402625/1‑143 (MQ=255)
      ttGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATcc                                                                                                                                           <  2:1417957/149‑1 (MQ=255)
          aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCCGGGTGGGTTTAAGACAACCCCCCCCCCATTTTCCCAAAACCCATAAAACGCGGGACCCCCCGGGAACGGCCCAGAAAATTTGGGGGGGAAGTAATCCACCaa                                                                                                                                      >  2:1492910/1‑150 (MQ=255)
            gTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATcc                                                                                                                                   <  1:998269/151‑1 (MQ=255)
             ttCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCTGGGCCGGTTTAAAACACTCCCCCCCCCATTTTGCAAAACGGCAAAAAAAGCGGGACCCCCGGAGAAAAGGCCGCAAATTTTGGGGGGAATGTaaaacaccaaacc                                                                                                                                   >  1:686513/1‑139 (MQ=255)
               cACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTg                                                                                                                                <  2:267366/151‑1 (MQ=255)
               cACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTg                                                                                                                                <  1:475644/151‑1 (MQ=255)
               cACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTg                                                                                                                                <  2:1460825/151‑1 (MQ=255)
                  ccACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCTGGTCGGGGTTAAGACATCCCCCCCCCCTTTTTCCAACACCTCCTAAAACGCGGGCCTCCCTGGGAACGGGCCGTAAATTTTGGGGGGAATTTAAAcccacaaacc                                                                                                                                   >  2:1523873/1‑137 (MQ=255)
                     cTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTTGTCGGGTTTACGACAACCGCCCCCCAATTTTCCCACACATCATAACGCGCGGCACCCCCTGAGAAAGGCCAGGAAAATTGGGGGGGAAAGTAAACCAGCAAACCTTTTTTCAc                                                                                                                          >  2:628855/1‑151 (MQ=255)
                       ttATAGAAGGCGATACCGATATACAGGCGGGGCCTGGTCGGGTTAACAAAATCCCCCCACCCTTTTTGCCACACCCCATAAAACGCGGCGCCCCGCGGGGAAGGCCAGTAAAATTTGGGGGGGATGTAATCCAGCAATCCTTGTTTCTccc                                                                                                                        >  1:648041/1‑151 (MQ=255)
                       ttATAGAAGGCGATACCGATATACAGGCGGGGCCCCGTCGGGTTTACGACCTCCGCCCACCAATTTTCCAAAACGCCAAAAACCGCGGCGCCCCGGGAGAAAGGCCCGGAAAATTTGGGGGGGAAGTAAACCAACCAACCTTTTTCCAccc                                                                                                                        >  2:609026/1‑151 (MQ=255)
                        tataGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCAccc                                                                                                                        <  2:368981/150‑1 (MQ=255)
                         ataGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCAt                                                                                                                      <  1:1098243/151‑1 (MQ=255)
                           aGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATc                                                                                                                     <  2:1655041/150‑1 (MQ=255)
                           aGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATc                                                                                                                     <  3:362559/150‑1 (MQ=255)
                             aaGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGAc                                                                                                                  <  1:1483579/151‑1 (MQ=255)
                               ggCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACgg                                                                                                                <  1:1604313/151‑1 (MQ=255)
                                gCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACggg                                                                                                               <  2:178038/151‑1 (MQ=255)
                                    tgccGATGTACAGGCAGGTCCTGGTCAGTTTAACGACATCCGCGCACCATTTTACCAGCACGTCATAGCGCGCGGCTCTCCGTGAAAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTACCACCCATCGACGGgtgt                                                                                                            <  4:178062/148‑1 (MQ=255)
                                     aCCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCa                                                                                                          <  3:435135/151‑1 (MQ=255)
                                      ccGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAg                                                                                                         <  1:568855/151‑1 (MQ=255)
                                       cGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCTCCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAAAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGc                                                                                                        <  2:536284/151‑1 (MQ=255)
                                       cGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGc                                                                                                        <  4:408371/151‑1 (MQ=255)
                                         atatACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGt                                                                                                      <  2:1282442/151‑1 (MQ=255)
                                             aCAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGa                                                                                                  <  4:291762/151‑1 (MQ=255)
                                             aCAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGa                                                                                                  <  1:393062/151‑1 (MQ=255)
                                               aGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGGGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAAACCTTGTTCCACCCATCGACGGGTGTCAGCg                                                                                                       >  2:199232/1‑144 (MQ=255)
                                                ggCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTc                                                                                               <  1:706913/151‑1 (MQ=255)
                                                 gCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCg                                                                                              <  1:827628/151‑1 (MQ=255)
                                                 gCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCg                                                                                              <  1:295190/151‑1 (MQ=255)
                                                  cGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGt                                                                                             >  2:1205249/1‑151 (MQ=255)
                                                    gggTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTc                                                                                            >  2:1014007/1‑151 (MQ=255)
                                                     ggTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTc                                                                                            >  3:433831/1‑150 (MQ=255)
                                                      gTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAAAGTAATCCAGCAATCCTTGTTGCACCCGTCGACGGGTGTCAGCGTAGGATCCGTc                                                                                            <  1:637100/149‑1 (MQ=255)
                                                         cTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAgg                                                                                       >  4:283884/1‑151 (MQ=255)
                                                         cTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAgg                                                                                       >  1:438144/1‑151 (MQ=255)
                                                          tGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGc                                                                                      >  2:1703205/1‑151 (MQ=255)
                                                            ggTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGc                                                                                      <  2:488960/150‑1 (MQ=255)
                                                            ggTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCt                                                                                     <  1:880285/151‑1 (MQ=255)
                                                             gTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTGATAGGCTg                                                                                    <  1:39990/151‑1 (MQ=255)
                                                             gTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTg                                                                                    >  2:705148/1‑151 (MQ=255)
                                                                ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCcg                                                                                 >  2:1466188/1‑151 (MQ=255)
                                                                ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCcg                                                                                 >  2:1525560/1‑151 (MQ=255)
                                                                 gTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCcgc                                                                                <  1:1162107/151‑1 (MQ=255)
                                                                   ttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCgcgc                                                                              >  4:107212/1‑151 (MQ=255)
                                                                      aCGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCCGCGTAGGATTCGTCATAGGCTGCCGCGCCTc                                                                           <  1:1190095/151‑1 (MQ=255)
                                                                      aCGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTc                                                                           <  2:1111948/151‑1 (MQ=255)
                                                                       cGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                          >  2:1218746/1‑151 (MQ=255)
                                                                         aCATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCggg                                                                        >  1:1117540/1‑151 (MQ=255)
                                                                         aCATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCggg                                                                        >  1:1340141/1‑151 (MQ=255)
                                                                          cATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGCGTCAGCGTAGGATTCGTCATAGCCTGCCGCGCCTCggg                                                                        >  4:399557/1‑150 (MQ=255)
                                                                                 gCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCAATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCgg                                                                  >  1:829068/1‑151 (MQ=255)
                                                                                  cccACCATTTTGCCAACACGTCATAGCGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTATGATGCGTCATAGGCTGCCGTGCATCGGGt                                                                       <  1:1461082/145‑1 (MQ=255)
                                                                                   ccaccaTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTACTCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCg                                                                   >  1:1436147/1‑148 (MQ=255)
                                                                                   ccaccaTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGa                                                                 >  4:187399/1‑150 (MQ=255)
                                                                                     accaTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAAcc                                                              <  2:1631886/151‑1 (MQ=255)
                                                                                     accaTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAAcc                                                              <  4:180565/151‑1 (MQ=255)
                                                                                        aTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACCGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAGTCGTCATAGGCTGCCGCGCCTCTGGTATGGGAACCg                                                             >  4:124402/1‑149 (MQ=255)
                                                                                               cAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCg                                                    <  1:559226/151‑1 (MQ=255)
                                                                                                aaCACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAGTGTAAACCAGCAATCCTTGTTCCACCCAACGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCgt                                                   <  2:1552037/151‑1 (MQ=255)
                                                                                                 acacGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCgtg                                                  >  2:430451/1‑151 (MQ=255)
                                                                                                     cGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                >  3:127337/1‑150 (MQ=255)
                                                                                                     cGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                >  3:270987/1‑150 (MQ=255)
                                                                                                      gTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGAt                                              <  3:138854/151‑1 (MQ=255)
                                                                                                        cATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCg                                            >  2:1275148/1‑151 (MQ=255)
                                                                                                           aaCGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGtt                                         <  2:57533/151‑1 (MQ=255)
                                                                                                             cgcgcgGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTAc                                       <  2:585274/151‑1 (MQ=255)
                                                                                                              gcgcgGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACg                                      >  1:892832/1‑151 (MQ=255)
                                                                                                                  ggCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCcac                                  <  3:932/151‑1 (MQ=255)
                                                                                                                   gCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCcac                                  >  2:396588/1‑150 (MQ=255)
                                                                                                                      cTCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCCGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCCCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACg                               >  1:568331/1‑150 (MQ=255)
                                                                                                                       tCCGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGcc                             >  4:459681/1‑151 (MQ=255)
                                                                                                                        ccGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCt                            >  1:1290431/1‑151 (MQ=255)
                                                                                                                        ccGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCt                            <  1:93240/151‑1 (MQ=255)
                                                                                                                         cGTGAGAACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGCGTGATCGGTTACGCCACACGCCTg                           >  2:977525/1‑151 (MQ=255)
                                                                                                                            gagaACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCg                        >  2:1341972/1‑151 (MQ=255)
                                                                                                                             agaACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCgg                       >  2:837078/1‑151 (MQ=255)
                                                                                                                              gaACGGCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCggg                      <  4:415589/151‑1 (MQ=255)
                                                                                                                                   ggCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCGCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTaa                  <  1:567757/151‑1 (MQ=255)
                                                                                                                                    gCCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAac                 <  1:1534129/151‑1 (MQ=255)
                                                                                                                                       cAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAacac               >  2:499763/1‑151 (MQ=255)
                                                                                                                                       cAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAacac               <  2:481943/151‑1 (MQ=255)
                                                                                                                                        aGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACAcc              >  2:480039/1‑151 (MQ=255)
                                                                                                                                        aGTAAATGTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACAcc              <  1:24459/151‑1 (MQ=255)
                                                                                                                                          tAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCa             >  1:80191/1‑150 (MQ=255)
                                                                                                                                           aaaTTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCaaa           <  2:1597311/151‑1 (MQ=255)
                                                                                                                                           aaaTTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCaa            <  4:373017/150‑1 (MQ=255)
                                                                                                                                           aaaTTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCaa            >  3:446626/1‑150 (MQ=255)
                                                                                                                                            aaTTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAAt          >  1:1510498/1‑151 (MQ=255)
                                                                                                                                             aTTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAAtt         <  1:992058/151‑1 (MQ=255)
                                                                                                                                              tttGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTc        >  2:225376/1‑151 (MQ=255)
                                                                                                                                               ttGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTc        <  2:430280/150‑1 (MQ=255)
                                                                                                                                                tGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCg       >  3:41093/1‑150 (MQ=255)
                                                                                                                                                 ggggAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGa      <  2:41816/150‑1 (MQ=255)
                                                                                                                                                 ggggAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGAc     >  2:1010476/1‑151 (MQ=255)
                                                                                                                                                    gAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTcc  <  3:526696/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCC  >  NC_000913/1567671‑1567962

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