Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 602,635 A→G 15.6% E226G (GAG→GGG)  pheP → phenylalanine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913602,6350AG15.6% 178.6 / 2.5 75E226G (GAG→GGG) phePphenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base A (13/48);  new base G (11/2);  total (25/50)
Fisher's exact test for biased strand distribution p-value = 3.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAACCTGGTGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCG  >  NC_000913/602404‑602870
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
ttAACCTTGGTAACGTGCCGTTATATGGCGAAACCGGGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTAGGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCCCGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGAtt                                                                                                                                                                                                                                          <  1:995890/239‑1 (MQ=255)
         tGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTAc                                                                                                                                                                                                                                  <  2:256736/238‑1 (MQ=255)
             cGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGGAATTATTTTCTCCCTCGGCGGTCTGGGGCGGGTTTGGGTTAGTGGCg                                                                                                                                                                                                                            >  2:281745/1‑239 (MQ=255)
                 cgcTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCCGTATCGACAACCTCTGGCGGTACGGTGGGTTCTTCGCCACCGGCTGGAAGTGGCCGGTTTTTGTCCGTGCGGGGATTTTTTTCTCCCTTTGGGGTCTGGGGGCGGTTGGGGTTTATGCCGccgg                                                                                                                                                                                                                        >  1:158669/1‑236 (MQ=255)
                       tttGGCGAACCCGATTTTTGTTTTGGTTTTATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTTGTGCTGGTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAgcgcg                                                                                                                                                                                                                  <  1:502215/237‑1 (MQ=255)
                            cGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGGAATTTTGTTTTCCTTCGGCGGGCTGGGGCTTGTTGGGATTACTGCCCCCGAAAcgcgcggt                                                                                                                                                                                                             >  2:294975/1‑236 (MQ=255)
                             gAAACCGAGTTCTGGTTTGCGTTTATTAAAGGGCGGGCAATCATCGGTATGATCGGTTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTAACGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATc                                                                                                                                                                                                            <  2:416330/238‑1 (MQ=255)
                              aaactggtttttGGTTTGGGTTGATTAAAGTTTCGGCAATCACTGGTTTGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATcc                                                                                                                                                                                                           <  1:68573/227‑1 (MQ=255)
                              aaaCCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTCTCGCACACCTCTGGCGCTACGGTGGTTTCTTTGCCCCCGGCTGGAATTGGCTGATTTTTTCGCTGGCCGGCAATTTGTTTTCCCTCCGTGGGTCTGGGGCTGTTTGGGGTTATGCCGccggaaagcgtcgatcg                                                                                                                                                                                                           >  1:120911/1‑223 (MQ=255)
                              aaaCCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTGGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGCCGGGAAATATGTTTTCCTTTGGCGGGCTGGAGGTGGTTTGGGTTACTTCCGGCTGAGCGCGCGATCCg                                                                                                                                                                                                          >  2:24623/1‑239 (MQ=255)
                              aaaCCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGGATGGGCTGGTTTTGTCGCTGGCGGGAATTATTTTTTCCTTCTGGGGGGTGGGGCGGTTTGGGGTTaacagcgctgaaatggcgggaccg                                                                                                                                                                                                            >  1:787940/1‑216 (MQ=255)
                              aaaCCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTTTGTTTTCCTTCGGCGGTCTTGGGCTGATTTGGGTTTCCGCCGCTGAAGCGCGCGATccc                                                                                                                                                                                                          >  2:105074/1‑238 (MQ=255)
                              aaaCCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCCCTGGCGCTACGGGGGTTTCTTCGCCACCCGCTGGAATGGGCTGGTTTTTTCGCTGGCGGGTAATTTGTTTTCCTTTGGCGCTCTGGGGGGGATTGGGGGTACTTCCGCGGAAGCGCGCGATCCg                                                                                                                                                                                                          >  1:356891/1‑239 (MQ=255)
                                     gTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGGAATTTTGTTTTCCTTTGGCCGGCTGGGGCTGCTTTGGGTTCACGCCCCGTAAACGCGCGATCCGGgcaaaa                                                                                                                                                                                                   >  2:844673/1‑239 (MQ=255)
                                       tgtgtgttGCGTTTATTAAAGTGGTGGCAATCCGCGGTATGATCGGTTTTGGCCTGTGGCTACTGTTTTCTGGGCCCGGCGGCGAGAAAGCCAGTGTCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCTGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAg                                                                                                                                                                                                  <  1:988049/232‑1 (MQ=255)
                                       tCTGGTTTGCCTTGATTAAGGTGCGGGCAATTATCGGTATGATCGGTTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAg                                                                                                                                                                                                  <  2:1014170/238‑1 (MQ=255)
                                               gggTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCaaa                                                                                                                                                                                       <  2:24034/237‑1 (MQ=255)
                                               gCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCaaa                                                                                                                                                                                       <  2:620074/239‑1 (MQ=255)
                                               gCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCaaa                                                                                                                                                                                       <  2:397494/239‑1 (MQ=255)
                                               gCGTGGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCGTGTGGCTGGTGTTTGCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCaaa                                                                                                                                                                                       <  2:581435/239‑1 (MQ=255)
                                                cGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTTTCCTTCCGCGGTTTGGGGCTGGTTGGGATTACTGCCGCCGAAACGCCCGATCCGGGAAAAAGCATTCCcaaa                                                                                                                                                                                      >  2:618514/1‑239 (MQ=255)
                                                     tttaaatgtttGGCATTCGTCGGTTGGATTGGTTTTGGCCTGGGGGTGTCGTTTGCTGGTCCCGGCGGCGAGAAAGCCAGTATCGACAACCTTTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGt                                                                                                                                                                                 <  1:146332/229‑1 (MQ=255)
                                                     tttAAAGTTCGGGCAATCATCGGGATGAGCGGCTTTGGCCTGTGGCTGCTGTTTTGTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGt                                                                                                                                                                                 <  1:787979/238‑1 (MQ=255)
                                                      ttAAATTGCTGGCAATTATCGGTTTTATCGGTTTTGGCCTGGGGGTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTAGCGACAACCCCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTa                                                                                                                                                                                <  1:521363/239‑1 (MQ=255)
                                                          aGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATc                                                                                                                                                                            <  2:158669/239‑1 (MQ=255)
                                                             tcgggCAATCCTCGGTATGATCGGCTTTGGCCTGTGGCTGGTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAgg                                                                                                                                                                         <  1:198279/236‑1 (MQ=255)
                                                             gCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAgg                                                                                                                                                                         <  2:815359/239‑1 (MQ=255)
                                                               tGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGGCTGGAGCTGATTGGGGTTACTGCCGCTGAAAGGGCCGATCCGGAGAAAAGCATTCCCAAAACGGGAAATCacggg                                                                                                                                                                       >  2:198279/1‑235 (MQ=255)
                                                               tGGCAATCAGCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAggtg                                                                                                                                                                       <  2:134741/239‑1 (MQ=255)
                                                               gggCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAggtg                                                                                                                                                                       <  2:317048/238‑1 (MQ=255)
                                                                ggCAATCATCGGTATGATCGGCTTGGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAggtg                                                                                                                                                                       <  2:194494/238‑1 (MQ=255)
                                                                           gTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGGTTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGTGCTGATTGGGGTTACTGCCGCTGAAGCGCGCGGTCCGGAAAAAAGCATTCCCAAAGCGGTAAAtacggggggggattgcac                                                                                                                                                            >  2:396433/1‑221 (MQ=255)
                                                                                         tgcttgtgGCTGCTTTTTTTGGTTCTCGGGGGCGAGAAAGCCAGTGTCGCGACCCTCTGGCGCTGCGGTGGTTTCTACGCCACCGGCTGGAATGGGCTGATTGTGTCGCTGGCGGTAAGTATGTTCTCCTTCGGCGGGCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGGGGTGTATCGCATCCTGCTGTTTTACAt                                                                                                                                             <  1:515789/235‑1 (MQ=255)
                                                                                         gtgtttttgttgtttttttCTTGGCCCGGGGGCGAGAAAGCCAGTATCGAACACCTTCGGCGCTCAGGTGGTTTCTTCGCCCCCGGCTGGAATGGGCTTATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACAt                                                                                                                                             <  1:586123/224‑1 (MQ=255)
                                                                                         ggCCTGTGGCTGCTGTTTTCGGGTCCGGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTTCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACAt                                                                                                                                             <  2:483272/239‑1 (MQ=255)
                                                                                              gtgGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTTTCGCTGGGGGTAATTTTTTTTTCCTTCGGCGGGCTGGGGCTGGTTGGGGTTACTGCCCCTGAAGCCCCCGGTCCGGCAAAAAACATTCCAAAAAGCGGAAATCAGGGGGGGGATCggcacccgccggtttaccaaggtt                                                                                                                                        >  1:102075/1‑216 (MQ=255)
                                                                                               tgGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGGAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTTGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAAACAGGGGGGGTATCGCATCCTGCTGTTTTACATCGGtt                                                                                                                                        >  2:793891/1‑238 (MQ=255)
                                                                                                   ttcttttttttCGTCTCCGCGGCGAGAAACCCAGTATTGGTAACCTCTGCGGCGCCGGGGGTTTTCTCGCCCAGGGCTTGAAGGGGGGGAATTTGGCGCTGGGGGTAATTATGGTCTCCTTCGGGGGTCTGGATATGATTGGGATTACTGGCGATGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtg                                                                                                                                   <  1:615222/234‑1 (MQ=255)
                                                                                                     ctgTTTTCTGGTCGCGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtggt                                                                                                                                 <  2:195419/239‑1 (MQ=255)
                                                                                                         tttCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGGGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCCAAAGCGGTAAATCAGGGGGGGTATCGCATCCTGCTGGTTTACATCGGTTCACtggtgggt                                                                                                                             >  2:531637/1‑237 (MQ=255)
                                                                                                                 cccgggggCGGGAAAGCCAGGTTCGCAAACCTTGGGCGCTACGGTGGTTTCTTCGCCACGGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTgg                                                                                                                      <  1:618514/236‑1 (MQ=255)
                                                                                                                 cACGGGGGCGAGAAAGCCAGTATCGCCGTCCGCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTGTGTCGCTGGCGGTAATTACGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTgg                                                                                                                      <  1:24623/238‑1 (MQ=255)
                                                                                                                 cACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTgg                                                                                                                      <  2:163456/238‑1 (MQ=255)
                                                                                                                  aCGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGGAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGGAAAAAGCATTCCCAAAACGGGAAAACCAGGGGGGGATTCGAACCCTCCGGTTTACATCGGGTCACTGGGGGGTTTTCttgg                                                                                                                     >  2:19224/1‑235 (MQ=255)
                                                                                                                    ggcggcGAGAAAGCCAGTTTCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGct                                                                                                                  <  2:255591/239‑1 (MQ=255)
                                                                                                                        gcgAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTa                                                                                                               <  2:970591/238‑1 (MQ=255)
                                                                                                                             aaaGCCATTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCgt                                                                                                          <  2:879170/238‑1 (MQ=255)
                                                                                                                                       tcaaCAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGGCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTg                                                                                               <  2:884946/236‑1 (MQ=255)
                                                                                                                                            aaCCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATc                                                                                          <  2:356891/239‑1 (MQ=255)
                                                                                                                                                  gggggcTTCGGTGTTTTTTCTGCCCACGGCGGGAGTGGGCTGGTTTTGTCGGTGGCGGTAATTAGGTTCTCCTTCGGCGGTCTGGAGCGGATTGGGATTACTGCCGATGAAGCGCGAGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCAATGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAg                                                                                    <  1:268987/235‑1 (MQ=255)
                                                                                                                                                  gggcggttgggttgttttgttgCCACCGGCTGGAATGGGGCGATTTTGGCGTGGGCGGTAAGTATGTTCTCCTTCGGCGGTCTGGAGCTGATGGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTTTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAg                                                                                    <  1:556218/218‑1 (MQ=255)
                                                                                                                                                                     tcGCCCCCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGAtt                                                                 <  2:787940/239‑1 (MQ=255)
                                                                                                                                                                                  gAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTgg                                                                                                                      <  2:206120/173‑1 (MQ=255)
                                                                                                                                                                                  gAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAAtctc                                                     <  2:551598/238‑1 (MQ=255)
                                                                                                                                                                                          cggtgttttgtgtgtgCGGTAATTACGTTCTCCGTCGGCGGTTCGGGTCTGATTGGGAATACCGCCGCTGAAGCGCGCGATCCGGCAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGGTTTACATCGGTTCACTGGTGGGTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCa                                              <  1:472994/223‑1 (MQ=255)
                                                                                                                                                                                           tGATTTTGTCGCTGGCGGTAATTAGGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAAc                                            <  2:102075/239‑1 (MQ=255)
                                                                                                                                                                                            gaatttGGCGGTGGCGGTAATTATGTTCTCCTTCCGCGGTCTGGAGCTGAATGGGATTACTGCCGGTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAAc                                            <  2:447459/235‑1 (MQ=255)
                                                                                                                                                                                                  ttgtcccgtcggttatttATGTTTTCCTTCGGCGGTCTGGAGCTGTTTGGGATTCCTGTCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCGGGTGGTGTATCGCATCCTGCTGTTTTACATGGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGt                                       <  1:820823/223‑1 (MQ=255)
                                                                                                                                                                                                  gTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGGTCCGGAAAAAAGCATTCCAAAAGCGGGAAATCAGGTGGTGCATCGGATCCTGCCGTTTTACATCGGTTCACTTGTGGGTTTACTGGTGCTCTAGCCGTGGGTGGAAGTGAAATCCCACCGCACCCCGTTTGTGAGGGTTTTTCATaaattcgcagccaaagggtgaa                                     >  1:706233/1‑219 (MQ=255)
                                                                                                                                                                                                            cgggtaTTTTGTTCCCCTTCGGCGGGTTGGGGCTGGTTGGGATTACTGCCGGTGAAGGCGCGGAGGCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTTTAACCCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCTCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTgcg                             <  1:710678/234‑1 (MQ=255)
                                                                                                                                                                                                                 atttttttCTCCTTCGGGGGTTTGGAGTTTATTGGGATTTCTGCCGATGAAGCGCGGGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGaa                        <  1:19224/233‑1 (MQ=255)
                                                                                                                                                                                                                     tGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGCTCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGGGGTGTATCGCATCCTGCCGTTTTACATAGGTTCACTGGTGGTATGAATGGCGATCTATCCTCGGGTGGGAGGGAAATCCACTAGTAGACCGTTTGGGATGATTTTCCTTAATCTCGACAGCACAGGGGAAGCTACTgggcgagacttt                    >  1:803503/1‑231 (MQ=255)
                                                                                                                                                                                                                         tccccgtccggggTCTGGAGCTTATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTc                 <  2:297861/228‑1 (MQ=255)
                                                                                                                                                                                                                         ctcCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGGGGTGTATCGCATCCTGCTGTTTTACATCGGTTAACTGTTGGTTTTACTGGCGCTCTTTCCGGGGGGGGAAGGGAAATCCAACAGTAGCCCGTTTGTGATGCTTTTCCATAATCTCGCAACCAACGTGGTAGCTTTTCGGGTGAACtttcgca                >  1:297861/1‑234 (MQ=255)
                                                                                                                                                                                                                          tcCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCAt               <  2:365540/239‑1 (MQ=255)
                                                                                                                                                                                                                          tcCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTAACTGGTGGCTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGGTTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCAt               >  2:860093/1‑239 (MQ=255)
                                                                                                                                                                                                                             cTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGACAGAAGCATTCCAAAAGAGGTAAATCAGGTGGAGTATCGCATCCTGCTGTTTTACATCGGTTCCCTGGCGGTTTTGCTGGCGCTCTATCCGTGGGTCGAAGTGAAATCCAACAGTCTCACGTTTGTGAGGCTTTTCCCCAATCTCGCCCCCACCGAGGTAGCTTCAACGGTGCACTTctccatct             >  1:134567/1‑232 (MQ=255)
                                                                                                                                                                                                                             cTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGACGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGCCAGCAACGGGGTAGCTTCTGCGCTGGACTTCGTCAt               >  2:973947/1‑237 (MQ=255)
                                                                                                                                                                                                                             cTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGGTTGTGAGGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCAt               >  2:1020996/1‑237 (MQ=255)
                                                                                                                                                                                                                                  gcggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACt                                                                                                                                    <  1:923486/115‑1 (MQ=255)
                                                                                                                                                                                                                                  gcggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACt                                                                                                                                    >  2:923486/1‑115 (MQ=255)
                                                                                                                                                                                                                                   cggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTa        <  2:1024500/238‑1 (MQ=255)
                                                                                                                                                                                                                                     gTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAgg                                                                                                                                                                         <  1:339455/75‑1 (MQ=255)
                                                                                                                                                                                                                                     gTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAgg                                                                                                                                                                         >  2:339455/1‑75 (MQ=255)
                                                                                                                                                                                                                                        tGGAGCTGATTGGGATTACTGCCGGTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCg  <  2:469752/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
TTAACCTGGTGAACGTGCGCTTATATGGCGAAACCGAGTTCTGGTTTGCGTTGATTAAAGTGCTGGCAATCATCGGTATGATCGGCTTTGGCCTGTGGCTGCTGTTTTCTGGTCACGGCGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCG  >  NC_000913/602404‑602870

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