Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,567,818 T→A 8.2% I336F (ATT→TTT)  yddW ← liprotein, glycosyl hydrolase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,567,8180TA 8.2% 252.6 / 2.3 74I336F (ATT→TTT) yddWliprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base T (20/47);  new base A (5/2);  total (25/49)
Fisher's exact test for biased strand distribution p-value = 3.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCTC  >  NC_000913/1567587‑1568054
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
tGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGGTTTACTAAATACGCCCCCCCTTTTTCCCACACGCCCTAAACCGCGGCACTCCC‑TGAGAACGGCCC‑GTAAATTTTGGGGGCAATGTAATcc                                                                                                                                                                                                                                       >  2:20107/1‑239 (MQ=255)
                cacaGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTcc                                                                                                                                                                                                                        <  2:869298/238‑1 (MQ=255)
                 acaGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCAc                                                                                                                                                                                                                      <  2:805012/239‑1 (MQ=255)
                        cGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGa                                                                                                                                                                                                               <  2:542213/239‑1 (MQ=255)
                               aTCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCTGGG‑CGGTTTTACGACAACCGCCCCCCACTTTTCTACCACATCCTAAAGAGCGGCCCCCCC‑CGAGAACGGCGC‑GTAACTTTTGGGTGGAAAGGTATCCACCAACCCTTGTTTCCCCCCTCGACGGgggt                                                                                                                                                                                                        >  1:421924/1‑239 (MQ=255)
                                tggggcggctTTTGCAGATCCTGTCCGCCGCTGTTAATCATCCAGTCTGGCTCTAACTTTGAAGGTTCCCCCACTGTATAGAAGGCGGTAACGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACTATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCGTGTTCCACCCATCGACGGgtgt                                                                                                                                                                                                        <  1:84544/231‑1 (MQ=255)
                                       gctTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCCCGTCCGGTTTAA‑GAAATACGCCCCCCACTTTTTCAACACGTCATAACAAGCGGCCCCCCT‑GGGGCACCGCCACGACAATTTGGGG‑GAAAAGTAATTAAGCAATCCCTGTTTCACCCCTCGCCgcggggatacggag                                                                                                                                                                                                >  1:878543/1‑226 (MQ=255)
                                          ttttCAGTTCCCGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTGTGAAGGTTCACCCACTTTATAGAAGGCGATTCCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAt                                                                                                                                                                                             <  2:183166/239‑1 (MQ=255)
                                          ttttCAGTACCGGGTCGCCGCCGTTAATCACCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAGGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAt                                                                                                                                                                                             <  1:1179744/239‑1 (MQ=255)
                                                         cgccgcCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTTACGACAACCGCCCCCCATTTTTCCAACCCATCCTAACGCGCGGCGCTCCGCTGGGAAC‑GCCC‑GGAAATTTTGGGAGCAAAGGAAATCACCAATCCTTGTTTCTCCCCTCGACGGGGGGCCGCCTAGGATTCGGTATAGCCTgg                                                                                                                                                                               >  1:190039/1‑237 (MQ=255)
                                                               cgTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCt                                                                                                                                                                        <  2:930000/239‑1 (MQ=255)
                                                                 ttAATCATCCAGTCTGGCTCTATCGTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                                                                                                                      <  2:255505/239‑1 (MQ=255)
                                                                       atcCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACg                                                                                                                                                                                                             <  2:862969/194‑1 (MQ=255)
                                                                                            gAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGa                                                                                                                                           <  2:1136760/239‑1 (MQ=255)
                                                                                              aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGAt                                                                                                                                          <  2:878543/238‑1 (MQ=255)
                                                                                              aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGAt                                                                                                                                          <  2:1047734/238‑1 (MQ=255)
                                                                                                 ttCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGGTTTACGACATCCGCCCACCATTTTTCCCACACGGCATAACGCGCGGCACTCCC‑GGGGAACGGCCA‑GCAAATTTTGGGGGGAATGTCATTCAACAATCCTTGTTTCCCCCACCGGCGGGGGGTAGCGTGAGGTTCGTCCGAGGCTGCCGCCCCTCCGCGATCGGGACCGGGCGGAtgggtggattggt                                                                                                                                      >  1:1047734/1‑227 (MQ=255)
                                                                                                          tttATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAg                                                                                                                                                       <  2:867818/213‑1 (MQ=255)
                                                                                                          tttATAGAAGGCGATACCGATATACAGGCGGGGCCTGGT‑CGGTTTTAACACATCCCCCCCCCCCTTTTTCAAAACCTCCTAAAACGCGGCCCTCCC‑TGCGGACAGCCA‑GTAAATTTGGGGGGGAAAGTAATCCAGCAAACCCTTTTTCCCCCACCACCCGGGGTCCGCGCAGGATTTGTCATAGGCTTCCCCCCCCCGCGTATCGGAACCGAGCGGAGCgtgcgtactgttctgcccc                                                                                                                              >  1:371957/1‑222 (MQ=255)
                                                                                                           ttATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACCAGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCGCCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCacac                                                                                                                            <  2:729629/239‑1 (MQ=255)
                                                                                                                 aaGGCGATCCCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGc                                                                                                                      <  1:20107/239‑1 (MQ=255)
                                                                                                                 aaGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCTTAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCTGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGc                                                                                                                      <  2:697489/239‑1 (MQ=255)
                                                                                                                  aGGCGAGACCTATATACGGGCGGGTCCTGGG‑CGGTTTAACGACAGCCGCCCACCATTGTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGc                                                                                                                      <  1:882054/238‑1 (MQ=255)
                                                                                                                                          tCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTGGCCAACACGTCATAACGCGCGGCACGCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTcc                                                                                              <  1:1029155/238‑1 (MQ=255)
                                                                                                                                               gT‑CGGTTTAACGACATCCGCCCACCATTTTGACAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGTCTGCCGCGCCTCGGGTATCGGAACCGAGCGGAACGTGTGATTGGTTACGCAACATGCCCGCCGGGCTAACACCAACTTCGACTCCCGGCtt                                                                                        >  1:964703/1‑239 (MQ=255)
                                                                                                                                                  cGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  2:190039/238‑1 (MQ=255)
                                                                                                                                                  cGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  2:1044698/238‑1 (MQ=255)
                                                                                                                                                  cGGTTCAACGACATCCGCCCACCATTTTTCCAACACGTCATACCGCGCGGCACGCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  1:317329/238‑1 (MQ=255)
                                                                                                                                                  cGGTGTAACGTCATCCGCCTCCCATTTTGTCAACACGTCGTAACGCGCGGCACCCCC‑GTATAACGGCCT‑GTAAATGTGGTGAGCAATGGAATACAGCAATCCTTGTTCCGCTCAACGACGGGTGCCAGCGTAGGATTCGTCATAGGCTGCCGCGCCCCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCGTGCCTGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  1:63688/238‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCGCCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  1:166952/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  2:821592/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  2:457833/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGAAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCAGGCTTAAt                                                                                     <  1:1060638/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGAGGGTCACGCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAAGCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  2:710392/239‑1 (MQ=255)
                                                                                                                                                   ggTTAAACCACATCCGCCCACCAGTTTGCCAACAGCTCATAACGCGCGGCGCCCCG‑GAAGAACGGCCA‑GTAAATTTGGGGAGGAAAGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGGCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTTTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  1:673601/239‑1 (MQ=255)
                                                                                                                                                    gTTTAACGAAATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  1:528144/238‑1 (MQ=255)
                                                                                                                                                     tttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                                                                                                            >  2:445215/1‑182 (MQ=255)
                                                                                                                                                     tttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                                                                                                            <  1:445215/182‑1 (MQ=255)
                                                                                                                                                     tttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGGCTCGGAACCGAGCGGATCGTGTGATAGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGCCTTAAt                                                                                     >  2:673601/1‑237 (MQ=255)
                                                                                                                                                           gACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTaa                                                                             <  2:12587/239‑1 (MQ=255)
                                                                                                                                                             cATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATg                                                                           <  2:1035307/239‑1 (MQ=255)
                                                                                                                                                              aTCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          <  2:762215/239‑1 (MQ=255)
                                                                                                                                                              aTCCGCCCAACATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          <  2:1117521/239‑1 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                                                                                                                      <  2:203875/146‑1 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                                                                                                                      >  1:203875/1‑146 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACCCGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  1:964533/1‑239 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  2:329657/1‑239 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  1:29755/1‑239 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  2:851778/1‑239 (MQ=255)
                                                                                                                                                                      accaTTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAGCGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGAt                                                                  <  1:422552/239‑1 (MQ=255)
                                                                                                                                                                      accaTTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGAt                                                                  <  2:421924/239‑1 (MQ=255)
                                                                                                                                                                        caTTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATAc                                                                <  1:329657/239‑1 (MQ=255)
                                                                                                                                                                          ttttGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCt                                                              <  2:597633/239‑1 (MQ=255)
                                                                                                                                                                                  acacGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAAAGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGa                                                                                                                                           >  1:557916/1‑154 (MQ=255)
                                                                                                                                                                                  acacGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAAAGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGa                                                                                                                                           <  2:557916/154‑1 (MQ=255)
                                                                                                                                                                                        cATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGATCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCTGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActg                                                 <  2:1112431/238‑1 (MQ=255)
                                                                                                                                                                                        cATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActg                                                 <  2:1115216/238‑1 (MQ=255)
                                                                                                                                                                                         aTACCGCGCCGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTAGCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActgct                                               <  2:513677/239‑1 (MQ=255)
                                                                                                                                                                                         aTAACGCGCGGCACGCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATTTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGTACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActgct                                               <  2:964703/239‑1 (MQ=255)
                                                                                                                                                                                               cgcgGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTACCTGCTGAGTAt                                         >  2:156644/1‑239 (MQ=255)
                                                                                                                                                                                                gcgGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGCTACCTTTGACATTAACTTCTGAGTAtt                                        >  1:19967/1‑239 (MQ=255)
                                                                                                                                                                                                          cgggagaACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGcgccg                               <  2:621152/236‑1 (MQ=255)
                                                                                                                                                                                                             gagaACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCa                            >  2:607058/1‑239 (MQ=255)
                                                                                                                                                                                                             gagaACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCa                            >  2:32804/1‑239 (MQ=255)
                                                                                                                                                                                                                  cGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCt                        >  2:253175/1‑238 (MQ=255)
                                                                                                                                                                                                                  cGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCt                        >  1:626240/1‑238 (MQ=255)
                                                                                                                                                                                                                    gCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          >  2:857715/1‑186 (MQ=255)
                                                                                                                                                                                                                    gCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          <  1:857715/186‑1 (MQ=255)
                                                                                                                                                                                                                             aTTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGc             <  2:591536/239‑1 (MQ=255)
                                                                                                                                                                                                                                    gAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATgc       >  2:355732/1‑238 (MQ=255)
                                                                                                                                                                                                                                      gCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATgcgc     >  2:822042/1‑238 (MQ=255)
                                                                                                                                                                                                                                        aaTGTAATCCATCCATCCTTGTCCCGCCTAGCAGCGGGTGTAAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCTc  <  1:754124/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCt   >  1:1082097/1‑238 (MQ=255)
                                                                                                                                                                                                                                        aaTGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCTc  <  2:29755/239‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
TGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGT‑CGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCG‑TGAGAACGGCCA‑GTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCTC  >  NC_000913/1567587‑1568054

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