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.

CTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCC  >  NC_000913/1567671‑1567962
                                                                                                                                                       |                                                                                                                                                   
ctATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATg                                                                                                                                                 <  2:1010506/149‑1 (MQ=255)
  aTCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGt                                                                                                                                                <  2:437223/148‑1 (MQ=255)
  aTCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCC‑CCTCGGGTTTATTAC‑A‑ACCCCCCACCCCTTTTCCAACA‑ACTCCTAAAGCGCGGCACTCCGTGGGGAA‑GG‑CCCCGAAAATTGTGGTGGAATGt                                                                                                                                                >  4:117980/1‑148 (MQ=255)
      ttGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTGGC‑GGTCGGGTTTAAGAC‑A‑TCCGCCCCCCACTTTTCCAAAACCGTC‑TAAAACGCGGCACCCCCCGAGAACAGG‑CCGAAAAATTTGGGGGGAAAttaat                                                                                                                                              >  4:402625/1‑143 (MQ=255)
      ttGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATcc                                                                                                                                           <  2:1417957/149‑1 (MQ=255)
          aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCCG‑GGTGGGTTTAA‑GACAA‑CCCCCCCCCCATTTTCCCAAAA‑C‑CCATAAAACGCGGGACCCCCCGGGAAC‑GGCCCAGAAAATTTGGGGGGGAA‑GTAATCCACCaa                                                                                                                                      >  2:1492910/1‑150 (MQ=255)
            gTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATcc                                                                                                                                   <  1:998269/151‑1 (MQ=255)
             ttCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGGCCTGGGCCGGTTTAA‑AAC‑ACTCCCCCCCCCATTTTG‑CAAAA‑CGGCAAAAAAAGCGGGACCCCCGGAGAAA‑AG‑GCCGCAAATTTTGGGGGGAATGTaaaacaccaaacc                                                                                                                                   >  1:686513/1‑139 (MQ=255)
               cACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTg                                                                                                                                <  2:267366/151‑1 (MQ=255)
               cACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTg                                                                                                                                <  1:475644/151‑1 (MQ=255)
               cACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTg                                                                                                                                <  2:1460825/151‑1 (MQ=255)
                  ccACTTTATAGAAGGCGATACCGATATACAGGC‑GGGGCCT‑GGTCGGGGTTAAGAC‑A‑TCCCCCCCCCCTTTTTCCAACA‑CCTCCTAAAACGCGGGCCTCCCTGGGAAC‑GG‑GCCGTAAATTTTGGGGGGAATTTAAAcccacaaacc                                                                                                                                   >  2:1523873/1‑137 (MQ=255)
                     cTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑TGTCGGGTTTACGAC‑A‑ACCGCCCCCCAATTTTCCCACA‑CATCATAACGCGCGGCACCCCCTGAGAAA‑GG‑CCAGGAAAATTGGGGGGGAAAGTAAACCAGCAAACCTTTTTTCAc                                                                                                                          >  2:628855/1‑151 (MQ=255)
                       ttATAGAAGGCGATACCGATATACAGGC‑GGGGCCT‑GGTCGGGTTAACAAA‑A‑TCCCCCCACCCTTTTTGCCACA‑CCCCATAAAACGCGGCGCCCCGCGGGGAA‑GG‑CCAGTAAAATTTGGGGGGGATGTAATCCAGCAATCCTTGTTTCTccc                                                                                                                        >  1:648041/1‑151 (MQ=255)
                       ttATAGAAGGCGATACCGATATACAGGC‑GGGGCCC‑CGTCGGGTTTACGAC‑C‑TCCGCCCACCAATTTTCCAAAA‑CGCCAAAAACCGCGGCGCCCCGGGAGAAA‑GGCCCGGAAAATTTGGGGGGGAA‑GTAAACCAACCAACCTTTTTCCAccc                                                                                                                        >  2:609026/1‑151 (MQ=255)
                        tataGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCAccc                                                                                                                        <  2:368981/150‑1 (MQ=255)
                         ataGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCAt                                                                                                                      <  1:1098243/151‑1 (MQ=255)
                           aGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATc                                                                                                                     <  2:1655041/150‑1 (MQ=255)
                           aGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATc                                                                                                                     <  3:362559/150‑1 (MQ=255)
                             aaGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGAc                                                                                                                  <  1:1483579/151‑1 (MQ=255)
                               ggCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACgg                                                                                                                <  1:1604313/151‑1 (MQ=255)
                                gCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACggg                                                                                                               <  2:178038/151‑1 (MQ=255)
                                    tgccGATGTACAGGC‑AGGTCCT‑GGTCAGTTTAACGAC‑A‑TCCGCGCACCATTTTACCAGCA‑CGTCATAGCGCGCGGCTCTCCGTGAAAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTACCACCCATCGACGGgtgt                                                                                                            <  4:178062/148‑1 (MQ=255)
                                     aCCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCa                                                                                                          <  3:435135/151‑1 (MQ=255)
                                      ccGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAg                                                                                                         <  1:568855/151‑1 (MQ=255)
                                       cGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCTCCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAAAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGc                                                                                                        <  2:536284/151‑1 (MQ=255)
                                       cGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGc                                                                                                        <  4:408371/151‑1 (MQ=255)
                                         atatACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGt                                                                                                      <  2:1282442/151‑1 (MQ=255)
                                             aCAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGa                                                                                                  <  4:291762/151‑1 (MQ=255)
                                             aCAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGa                                                                                                  <  1:393062/151‑1 (MQ=255)
                                               aGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGGGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAAACCTTGTTCCACCCATCGACGGGTGTCAGCg                                                                                                       >  2:199232/1‑144 (MQ=255)
                                                ggC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTc                                                                                               <  1:706913/151‑1 (MQ=255)
                                                 gC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCg                                                                                              <  1:827628/151‑1 (MQ=255)
                                                 gC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCg                                                                                              <  1:295190/151‑1 (MQ=255)
                                                  c‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGt                                                                                             >  2:1205249/1‑151 (MQ=255)
                                                    gggTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTc                                                                                            >  2:1014007/1‑151 (MQ=255)
                                                     ggTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTc                                                                                            >  3:433831/1‑150 (MQ=255)
                                                      gTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAAAGTAATCCAGCAATCCTTGTTGCACCCGTCGACGGGTGTCAGCGTAGGATCCGTc                                                                                            <  1:637100/149‑1 (MQ=255)
                                                         cT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAgg                                                                                       >  4:283884/1‑151 (MQ=255)
                                                         cT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAgg                                                                                       >  1:438144/1‑151 (MQ=255)
                                                          t‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGc                                                                                      >  2:1703205/1‑151 (MQ=255)
                                                            ggTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGc                                                                                      <  2:488960/150‑1 (MQ=255)
                                                            ggTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCt                                                                                     <  1:880285/151‑1 (MQ=255)
                                                             gTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTGATAGGCTg                                                                                    <  1:39990/151‑1 (MQ=255)
                                                             gTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTg                                                                                    >  2:705148/1‑151 (MQ=255)
                                                                ggTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCcg                                                                                 >  2:1466188/1‑151 (MQ=255)
                                                                ggTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCcg                                                                                 >  2:1525560/1‑151 (MQ=255)
                                                                 gTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCcgc                                                                                <  1:1162107/151‑1 (MQ=255)
                                                                   ttAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCgcgc                                                                              >  4:107212/1‑151 (MQ=255)
                                                                      aCGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCCGCGTAGGATTCGTCATAGGCTGCCGCGCCTc                                                                           <  1:1190095/151‑1 (MQ=255)
                                                                      aCGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTc                                                                           <  2:1111948/151‑1 (MQ=255)
                                                                       cGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                          >  2:1218746/1‑151 (MQ=255)
                                                                         aC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCggg                                                                        >  1:1117540/1‑151 (MQ=255)
                                                                         aC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCggg                                                                        >  1:1340141/1‑151 (MQ=255)
                                                                          c‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGCGTCAGCGTAGGATTCGTCATAGCCTGCCGCGCCTCggg                                                                        >  4:399557/1‑150 (MQ=255)
                                                                                 gCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCAATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCgg                                                                  >  1:829068/1‑151 (MQ=255)
                                                                                  cccACCATTTTGCCAACA‑CGTCATAGCGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTATGATGCGTCATAGGCTGCCGTGCATCGGGt                                                                       <  1:1461082/145‑1 (MQ=255)
                                                                                   ccaccaTTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTACTCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCg                                                                   >  1:1436147/1‑148 (MQ=255)
                                                                                   ccaccaTTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGa                                                                 >  4:187399/1‑150 (MQ=255)
                                                                                     accaTTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAAcc                                                              <  2:1631886/151‑1 (MQ=255)
                                                                                     accaTTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAAcc                                                              <  4:180565/151‑1 (MQ=255)
                                                                                        aTTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑CG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAGTCGTCATAGGCTGCCGCGCCTCTGGTATGGGAACCg                                                             >  4:124402/1‑149 (MQ=255)
                                                                                               cAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCg                                                    <  1:559226/151‑1 (MQ=255)
                                                                                                aaCA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGCAGTGTAAACCAGCAATCCTTGTTCCACCCAACGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCgt                                                   <  2:1552037/151‑1 (MQ=255)
                                                                                                 aca‑cGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCgtg                                                  >  2:430451/1‑151 (MQ=255)
                                                                                                     cGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                >  3:127337/1‑150 (MQ=255)
                                                                                                     cGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                >  3:270987/1‑150 (MQ=255)
                                                                                                      gTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGAt                                              <  3:138854/151‑1 (MQ=255)
                                                                                                        cATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCg                                            >  2:1275148/1‑151 (MQ=255)
                                                                                                           aaCGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGtt                                         <  2:57533/151‑1 (MQ=255)
                                                                                                             cgcgcgGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTAc                                       <  2:585274/151‑1 (MQ=255)
                                                                                                              gcgcgGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACg                                      >  1:892832/1‑151 (MQ=255)
                                                                                                                  ggCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCcac                                  <  3:932/151‑1 (MQ=255)
                                                                                                                   gCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCcac                                  >  2:396588/1‑150 (MQ=255)
                                                                                                                      cTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCCGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCCCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACg                               >  1:568331/1‑150 (MQ=255)
                                                                                                                       tCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGcc                             >  4:459681/1‑151 (MQ=255)
                                                                                                                        ccGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCt                            >  1:1290431/1‑151 (MQ=255)
                                                                                                                        ccGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCt                            <  1:93240/151‑1 (MQ=255)
                                                                                                                         cGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGCGTGATCGGTTACGCCACACGCCTg                           >  2:977525/1‑151 (MQ=255)
                                                                                                                            gagaAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCg                        >  2:1341972/1‑151 (MQ=255)
                                                                                                                             agaAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCgg                       >  2:837078/1‑151 (MQ=255)
                                                                                                                              gaAC‑GG‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCggg                      <  4:415589/151‑1 (MQ=255)
                                                                                                                                   gg‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCGCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTaa                  <  1:567757/151‑1 (MQ=255)
                                                                                                                                    g‑CCAGTAAATTTGGGGAGTAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAac                 <  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)
                                                                                                                                                       |                                                                                                                                                   
CTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGC‑GGGTCCT‑GGTCGGTTTAACGAC‑A‑TCCGCCCACCATTTTGCCAACA‑CGTCATAACGCGCGGCACTCCGTGAGAAC‑GG‑CCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCC  >  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: ‑