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

TGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCTC  >  NC_000913/1567587‑1568054
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
tGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGGTTTACTAAATACGCCCCCCCTTTTTCCCACACGCCCTAAACCGCGGCACTCCCTGAGAACGGCCCGTAAATTTTGGGGGCAATGTAATcc                                                                                                                                                                                                                                       >  2:20107/1‑239 (MQ=255)
                cacaGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTcc                                                                                                                                                                                                                        <  2:869298/238‑1 (MQ=255)
                 acaGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCAc                                                                                                                                                                                                                      <  2:805012/239‑1 (MQ=255)
                        cGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGa                                                                                                                                                                                                               <  2:542213/239‑1 (MQ=255)
                               aTCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCTGGGCGGTTTTACGACAACCGCCCCCCACTTTTCTACCACATCCTAAAGAGCGGCCCCCCCCGAGAACGGCGCGTAACTTTTGGGTGGAAAGGTATCCACCAACCCTTGTTTCCCCCCTCGACGGgggt                                                                                                                                                                                                        >  1:421924/1‑239 (MQ=255)
                                tggggcggctTTTGCAGATCCTGTCCGCCGCTGTTAATCATCCAGTCTGGCTCTAACTTTGAAGGTTCCCCCACTGTATAGAAGGCGGTAACGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACTATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCGTGTTCCACCCATCGACGGgtgt                                                                                                                                                                                                        <  1:84544/231‑1 (MQ=255)
                                       gctTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGGCCCCGTCCGGTTTAAGAAATACGCCCCCCACTTTTTCAACACGTCATAACAAGCGGCCCCCCTGGGGCACCGCCACGACAATTTGGGGGAAAAGTAATTAAGCAATCCCTGTTTCACCCCTCGCCgcggggatacggag                                                                                                                                                                                                >  1:878543/1‑226 (MQ=255)
                                          ttttCAGTTCCCGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTGTGAAGGTTCACCCACTTTATAGAAGGCGATTCCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAt                                                                                                                                                                                             <  2:183166/239‑1 (MQ=255)
                                          ttttCAGTACCGGGTCGCCGCCGTTAATCACCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAGGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGAt                                                                                                                                                                                             <  1:1179744/239‑1 (MQ=255)
                                                         cgccgcCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTTACGACAACCGCCCCCCATTTTTCCAACCCATCCTAACGCGCGGCGCTCCGCTGGGAACGCCCGGAAATTTTGGGAGCAAAGGAAATCACCAATCCTTGTTTCTCCCCTCGACGGGGGGCCGCCTAGGATTCGGTATAGCCTgg                                                                                                                                                                               >  1:190039/1‑237 (MQ=255)
                                                               cgTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCt                                                                                                                                                                        <  2:930000/239‑1 (MQ=255)
                                                                 ttAATCATCCAGTCTGGCTCTATCGTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                                                                                                                      <  2:255505/239‑1 (MQ=255)
                                                                       atcCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACg                                                                                                                                                                                                             <  2:862969/194‑1 (MQ=255)
                                                                                            gAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGa                                                                                                                                           <  2:1136760/239‑1 (MQ=255)
                                                                                              aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGAt                                                                                                                                          <  2:878543/238‑1 (MQ=255)
                                                                                              aGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGAt                                                                                                                                          <  2:1047734/238‑1 (MQ=255)
                                                                                                 ttCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGGTTTACGACATCCGCCCACCATTTTTCCCACACGGCATAACGCGCGGCACTCCCGGGGAACGGCCAGCAAATTTTGGGGGGAATGTCATTCAACAATCCTTGTTTCCCCCACCGGCGGGGGGTAGCGTGAGGTTCGTCCGAGGCTGCCGCCCCTCCGCGATCGGGACCGGGCGGAtgggtggattggt                                                                                                                                      >  1:1047734/1‑227 (MQ=255)
                                                                                                          tttATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAg                                                                                                                                                       <  2:867818/213‑1 (MQ=255)
                                                                                                          tttATAGAAGGCGATACCGATATACAGGCGGGGCCTGGTCGGTTTTAACACATCCCCCCCCCCCTTTTTCAAAACCTCCTAAAACGCGGCCCTCCCTGCGGACAGCCAGTAAATTTGGGGGGGAAAGTAATCCAGCAAACCCTTTTTCCCCCACCACCCGGGGTCCGCGCAGGATTTGTCATAGGCTTCCCCCCCCCGCGTATCGGAACCGAGCGGAGCgtgcgtactgttctgcccc                                                                                                                              >  1:371957/1‑222 (MQ=255)
                                                                                                           ttATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACCAGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCGCCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCacac                                                                                                                            <  2:729629/239‑1 (MQ=255)
                                                                                                                 aaGGCGATCCCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGc                                                                                                                      <  1:20107/239‑1 (MQ=255)
                                                                                                                 aaGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCTTAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCTGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGc                                                                                                                      <  2:697489/239‑1 (MQ=255)
                                                                                                                  aGGCGAGACCTATATACGGGCGGGTCCTGGGCGGTTTAACGACAGCCGCCCACCATTGTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGc                                                                                                                      <  1:882054/238‑1 (MQ=255)
                                                                                                                                          tCCTGGTCGGTTTAACGACATCCGCCCACCATTTGGCCAACACGTCATAACGCGCGGCACGCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTcc                                                                                              <  1:1029155/238‑1 (MQ=255)
                                                                                                                                               gTCGGTTTAACGACATCCGCCCACCATTTTGACAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGTCTGCCGCGCCTCGGGTATCGGAACCGAGCGGAACGTGTGATTGGTTACGCAACATGCCCGCCGGGCTAACACCAACTTCGACTCCCGGCtt                                                                                        >  1:964703/1‑239 (MQ=255)
                                                                                                                                                  cGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  2:190039/238‑1 (MQ=255)
                                                                                                                                                  cGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  2:1044698/238‑1 (MQ=255)
                                                                                                                                                  cGGTTCAACGACATCCGCCCACCATTTTTCCAACACGTCATACCGCGCGGCACGCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  1:317329/238‑1 (MQ=255)
                                                                                                                                                  cGGTGTAACGTCATCCGCCTCCCATTTTGTCAACACGTCGTAACGCGCGGCACCCCCGTATAACGGCCTGTAAATGTGGTGAGCAATGGAATACAGCAATCCTTGTTCCGCTCAACGACGGGTGCCAGCGTAGGATTCGTCATAGGCTGCCGCGCCCCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCGTGCCTGGCTAACACCAAATTCGACTCCCGGCTTa                                                                                       <  1:63688/238‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCGCCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  1:166952/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  2:821592/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  2:457833/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGAAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCAGGCTTAAt                                                                                     <  1:1060638/239‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGAGGGTCACGCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAAGCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  2:710392/239‑1 (MQ=255)
                                                                                                                                                   ggTTAAACCACATCCGCCCACCAGTTTGCCAACAGCTCATAACGCGCGGCGCCCCGGAAGAACGGCCAGTAAATTTGGGGAGGAAAGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGGCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTTTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  1:673601/239‑1 (MQ=255)
                                                                                                                                                    gTTTAACGAAATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAAt                                                                                     <  1:528144/238‑1 (MQ=255)
                                                                                                                                                     tttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                                                                                                            >  2:445215/1‑182 (MQ=255)
                                                                                                                                                     tttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGtgtg                                                                                                                                            <  1:445215/182‑1 (MQ=255)
                                                                                                                                                     tttAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGGCTCGGAACCGAGCGGATCGTGTGATAGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGCCTTAAt                                                                                     >  2:673601/1‑237 (MQ=255)
                                                                                                                                                           gACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTaa                                                                             <  2:12587/239‑1 (MQ=255)
                                                                                                                                                             cATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATg                                                                           <  2:1035307/239‑1 (MQ=255)
                                                                                                                                                              aTCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          <  2:762215/239‑1 (MQ=255)
                                                                                                                                                              aTCCGCCCAACATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          <  2:1117521/239‑1 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                                                                                                                      <  2:203875/146‑1 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCg                                                                                                                                                                      >  1:203875/1‑146 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACCCGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  1:964533/1‑239 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  2:329657/1‑239 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  1:29755/1‑239 (MQ=255)
                                                                                                                                                               tCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGgt                                                                         >  2:851778/1‑239 (MQ=255)
                                                                                                                                                                      accaTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAGCGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGAt                                                                  <  1:422552/239‑1 (MQ=255)
                                                                                                                                                                      accaTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGAt                                                                  <  2:421924/239‑1 (MQ=255)
                                                                                                                                                                        caTTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATAc                                                                <  1:329657/239‑1 (MQ=255)
                                                                                                                                                                          ttttGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCt                                                              <  2:597633/239‑1 (MQ=255)
                                                                                                                                                                                  acacGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAAAGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGa                                                                                                                                           >  1:557916/1‑154 (MQ=255)
                                                                                                                                                                                  acacGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAAAGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGa                                                                                                                                           <  2:557916/154‑1 (MQ=255)
                                                                                                                                                                                        cATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGATCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCTGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActg                                                 <  2:1112431/238‑1 (MQ=255)
                                                                                                                                                                                        cATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActg                                                 <  2:1115216/238‑1 (MQ=255)
                                                                                                                                                                                         aTACCGCGCCGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTAGCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActgct                                               <  2:513677/239‑1 (MQ=255)
                                                                                                                                                                                         aTAACGCGCGGCACGCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATTTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGTACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAActgct                                               <  2:964703/239‑1 (MQ=255)
                                                                                                                                                                                               cgcgGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTACCTGCTGAGTAt                                         >  2:156644/1‑239 (MQ=255)
                                                                                                                                                                                                gcgGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGCTACCTTTGACATTAACTTCTGAGTAtt                                        >  1:19967/1‑239 (MQ=255)
                                                                                                                                                                                                          cgggagaACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGcgccg                               <  2:621152/236‑1 (MQ=255)
                                                                                                                                                                                                             gagaACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCa                            >  2:607058/1‑239 (MQ=255)
                                                                                                                                                                                                             gagaACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCa                            >  2:32804/1‑239 (MQ=255)
                                                                                                                                                                                                                  cGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCt                        >  2:253175/1‑238 (MQ=255)
                                                                                                                                                                                                                  cGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCt                        >  1:626240/1‑238 (MQ=255)
                                                                                                                                                                                                                    gCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          >  2:857715/1‑186 (MQ=255)
                                                                                                                                                                                                                    gCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATgg                                                                          <  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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
TGCCGCTAATTTCGGGCACAGCATCGTTAAGATCGAGCTGCTTTTTCAGTTCCGGTACGCCGCCGTTAATCATCCAGTCTGGCTCTATCTTTGAAGGTTCACCCACTTTATAGAAGGCGATACCGATATACAGGCGGGTCCTGGTCGGTTTAACGACATCCGCCCACCATTTTGCCAACACGTCATAACGCGCGGCACTCCGTGAGAACGGCCAGTAAATTTGGGGAGCAATGTAATCCAGCAATCCTTGTTCCACCCATCGACGGGTGTCAGCGTAGGATTCGTCATAGGCTGCCGCGCCTCGGGTATCGGAACCGAGCGGATCGTGTGATCGGTTACGCCACACGCCTGCCGGGCTAACACCAAATTCGACTCCCGGCTTAATGCTTTTAATGGTGTGCGATACCTTTGCAATTAACTGCTGAGTATTGTTGCGCCGCCAGTCTGCTTTTGACGCAAATGCGCCTC  >  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: