Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,231 A→T 15.8% I416F (ATC→TTC)  proY → proline‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,2310AT15.8% 165.9 / 3.3 69I416F (ATC→TTC) proYproline‑specific permease
Reads supporting (aligned to +/- strand):  ref base A (14/43);  new base T (6/6);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 9.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  NC_000913/421996‑422466
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
tgtgccgtccggggtgttgTTTGCGGGGTTTTTGAACTACATCATGCCGGAAAACGGCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCCCGGTGTGGGTGTGGATTATGAGCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATc                                                                                                                                                                                                                                              <  2:862129/221‑1 (MQ=255)
gatgaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATc                                                                                                                                                                                                                                              <  1:409601/238‑1 (MQ=255)
      tACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCGGGTTTTCCTGCTCTTTATTTTTGGGtgtg                                                                                                                                                                                                                                       >  1:159248/1‑236 (MQ=255)
       tccggggttctgTTTGGGGTGTTACCGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGAt                                                                                                                                                                                                                                      <  1:396602/230‑1 (MQ=255)
       tccGCGCTTTTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTTTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGAt                                                                                                                                                                                                                                      <  2:334528/238‑1 (MQ=255)
       tccGCGCGGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGGGTTCGCTTCGCTGGGAACTTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGAt                                                                                                                                                                                                                                      <  1:871150/238‑1 (MQ=255)
        ccGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTCGTCTTTTTTATCgctttttct                                                                                                                                                                                                                                     >  1:73032/1‑231 (MQ=255)
                       acaggtttCTTAACTTCATCTTGCCGGAAACCGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGa                                                                                                                                                                                                                      <  2:335202/231‑1 (MQ=255)
                                  aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTTCCGCTCTTTTTTATTGTGTTGATTGTGTTTCACCCGGGTACGCGGAttt                                                                                                                                                                                                           >  1:358183/1‑236 (MQ=255)
                                   ggtatgtttataagggacaggccTTCCTGGTTATCCCTTCGCTGGCAACCTTCGCCACGGTGGGGGTGTGGATTATGATCCTGATGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  1:379289/217‑1 (MQ=255)
                                   acccccctctGCCGGGAAACGTTTTCTTTGGTATCTCTTCGCTGGCAACCTTCGCCACGGTGGGGGGGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  2:470361/230‑1 (MQ=255)
                                   accaactcctGCCGGAAAACGTTTTCTTGGTTATCGCTTCGCTGGCAACCTTCGCCACGGGGTGGGTTTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  2:43225/230‑1 (MQ=255)
                                    cTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTCCTGCTCTTTTTTATCCGGTTGTTTGGGTTTCTCCCGGAGATGCGCAtctc                                                                                                                                                                                                          >  1:62171/1‑237 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAGGTTAAGGCGCTGAAATTTAAAGTGCCGGGCGGGGCAGCAAGGACCATAGGCGGGGCGTGTTTTCTGGTCTTTTTTTATCGGGTTATTGGTGTACATCCCGGATACGCGTAtctc                                                                                                                                                                                                          >  2:172896/1‑236 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTGCTCTTTATTATCGGGGTGATTGGTTATTACCCGGATACGCGGAttcctttt                                                                                                                                                                                                      >  1:335202/1‑232 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGGGGGGTAGAACCGACCATCGGCGGGGTGGTTTTCCTGCGCTTTTTTTTTGGGGTGATTTGTTTTCACCCGGGTACGCGGATATCGccg                                                                                                                                                                                                      >  2:919905/1‑237 (MQ=255)
                                            tGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGGGATTTTCCTTCTCTTTTTTTTCGGGGTTATTTGTTTTCACCCCGTTAccgggttttcgctttattgc                                                                                                                                                                                                >  1:904204/1‑220 (MQ=255)
                                                    aaCTTCTTCCTGGTGGTTGTTTTGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTCTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  2:358183/239‑1 (MQ=255)
                                                    aaCTTCTTCCTGGTGATCGCCTCGCTGGCAACCTTCGCCACGGTGTTGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  1:24447/239‑1 (MQ=255)
                                                       gTCTTCCTGGTGATCGCTGCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTAGGATCCTGCTGTCGCAAATTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTg                                                                                                                                                                                     <  1:251731/239‑1 (MQ=255)
                                                        tCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTg                                                                                                                                                                                     <  2:440203/238‑1 (MQ=255)
                                                         cTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGa                                                                                                                                                                                   <  1:741344/239‑1 (MQ=255)
                                                            ccTGGTTATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttg                                                                                                                                                                                <  1:464870/239‑1 (MQ=255)
                                                            ccTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTTTTATCGGGTTGATTTGTTATCACCCGGATACGCGTATTTCCCTGTATGTTGGTTTCTCGGGGGttt                                                                                                                                                                                >  1:607044/1‑238 (MQ=255)
                                                              tggttgtcgctctgtgggCAACCTTCGCCACGGTGTGGGTTTGGATTATGTGCCTGCTTTCGCAAATTGCCTTCCGGCGCCGGTGGCCGGCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttgtt                                                                                                                                                                              <  2:401817/223‑1 (MQ=255)
                                                              tGGTGTTCGCTTCGCTGGCAACCTTCGCCTCGGTGTGGGTTTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttgtt                                                                                                                                                                              <  2:45501/239‑1 (MQ=255)
                                                                  tttcgcttcgttGGCAACCTTCGCCCCGGGGGGGGGGTGGATTATGATCCGGCTGGCGCAAATTGCCTTCCGTCGCCGTTTGCCGGCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                          <  2:468412/237‑1 (MQ=255)
                                                                   aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGGTACGCGTATTTCGCTGTATGTCGGTTTTGGGTGGGTTGTTTTGCt                                                                                                                                                                         >  1:45501/1‑239 (MQ=255)
                                                                          cgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCggg                                                                                                                                                                                                                                           <  1:541080/167‑1 (MQ=255)
                                                                                 aaCCTTCGCCACGGTGTGGGTGTGGATTATGGCCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTgg                                                                                                                                                           <  1:219375/238‑1 (MQ=255)
                                                                                         ccACGGGGTGGGTGTGGATTTTGATCCTGCTGTCGCAAATGGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTaaa                                                                                                                                                  <  1:364660/239‑1 (MQ=255)
                                                                                              gtgtgggtgtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCCTCCCTTTTTTATCGGGGTGGTTGGTTTTTACCCGGGTAGACCGATTTCCCTGTATGTCGGGTTTTGGTGGGTTTTTGTGTGGTTGtttttgtggggggtttaaaagc                                                                                                                                              >  2:438579/1‑219 (MQ=255)
                                                                                               tgtgggtgtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTAt                                                                                                                                                                                                             <  1:872126/175‑1 (MQ=255)
                                                                                               tgtgggtgtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTg                                                                                                                                                            <  1:115368/223‑1 (MQ=255)
                                                                                                  ggttgtggATTTTGATCCGGCGTTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGGGGTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATGGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                          <  2:766446/235‑1 (MQ=255)
                                                                                                  gggtgtggATTATGATCCTGCTGTCGAAAATTGCCTTCCGTCGGCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                          <  1:438579/238‑1 (MQ=255)
                                                                                                    gtttgtattgtttttctggttgcGCAATTTGCCTTCCGGCGCCGTTGGCCCCCAGAAGAAGTTAAGGGCTTTAAATTTAAAGGGCCGGGTGGGGTAGCAACGAGCATCGGCGGGCTGATTTTCCTGCTCTTTATGATCGGGTTGATTGGTTATCTCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACg                                                                                                                                        <  2:436948/216‑1 (MQ=255)
                                                                                                    gggtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACg                                                                                                                                        <  1:26787/236‑1 (MQ=255)
                                                                                                         gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGGGATTTTCCTGCTCTTTATTTTCGGGGTGGTTGGTTTTACCCCCGATATCCGGTCCTCCGCGTATGTCGGTTTGTCGTGGGTTTGTGTGCTGTTTTTTGGTGGGGTGTTTTAACGCCCCCCCGCTCggg                                                                                                                                  >  2:716029/1‑237 (MQ=255)
                                                                                                         gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGTGTTTTCTGCTCTTTATTTTGTGGTTGATTGGTTTTCACCCGGGTACGCGTATCTCGCTGTTTGTCGGGTTCGCGTGGGGTGTTTTGGTGTTGTTTTGGTGGGGGTTTAAACGCCCCCCCGAGCgcc                                                                                                                                  >  2:691401/1‑237 (MQ=255)
                                                                                                         gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGGCTCTCGCTGTATGTTGGGTTCGTGTGGGTTGTTTTGTTGTTGGTTGGCGGGGTGTTTaaaactgtatcttttaaaa                                                                                                                                  >  2:115368/1‑223 (MQ=255)
                                                                                                             aggactctgctgTCTCAAATTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCGCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAgct                                                                                                                              <  1:54631/233‑1 (MQ=255)
                                                                                                                    tggtgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGaa                                                                                                                      <  1:123786/236‑1 (MQ=255)
                                                                                                                    tgctgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGGTGATTGGTTATCACCCGGATACGCGTATCTCGCCTTATGTCCGTTTTGCGTGGGTTGTTGTGCTGTTGGTTGGCTTGGTGTTTTAACGCCGCCACGGTCGTCACCTTGCTGaa                                                                                                                      >  1:808258/1‑239 (MQ=255)
                                                                                                                      ctgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGTGAGTTTCATGCGCTTTATTATTGGGGTGATTGGTTATGACCCGGGAACGCGTATCTCCCTTGATGTCGGGTTTGCGTGGGGTGTTGTGGTGTTGATTGGGTGGGGTTTTAAAAGCCCCCACGagaggcgggcgggggaaaa                                                                                                                     >  2:517412/1‑221 (MQ=255)
                                                                                                                                               gTTTGCCGGCAGAAGAAGTTAAGGCGTTGAAATTTAAAGGGGCGGGTGGGGTAGCAACGACCATCGGGGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCGCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGc                                                                                           <  1:285816/239‑1 (MQ=255)
                                                                                                                                                tttGCCGCCAGAAGAAGTTAGGGGGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGcc                                                                                          <  1:784581/239‑1 (MQ=255)
                                                                                                                                                      gccAGAAGAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCa                                                                                    <  1:691401/239‑1 (MQ=255)
                                                                                                                                                      gccAAAAGAATTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCCTCGGCGGGCTGATTTTCCTGTTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCa                                                                                    <  1:896274/239‑1 (MQ=255)
                                                                                                                                                                    aggtgtttaaaTTTAAAGTGCGGGGTGGGGTGGCAACGACCTTGGGCGGGCTGTTTTTTCTGATCTTTATTTTGGTGTTTATTGTTTATCCCCCGGATTCGCGTCTCGTGCTGTAGGTCTGTTTCGGGTGGATTGTTGTGCTGTTGATGGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGa                                                                      <  2:808258/231‑1 (MQ=255)
                                                                                                                                                                          tGAAATTTAAGGTGCCGGGTTGGGTAGCAAGGACCATCGGGGGTCTGATTTTCCTGCTTTTTATTATCGGGTTTATTGGTTATCACCCGGATACGCTTATCTCGTTGTTTGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGc                                                                <  2:738237/239‑1 (MQ=255)
                                                                                                                                                                                       gCCGGGGTGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTTTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAg                                                     <  1:544255/237‑1 (MQ=255)
                                                                                                                                                                                          ggggtgggTAGCACCGACCAGCGGCGGGCTTATTTTCCTGCTTTTTATTATCGGGTTTATTGGTTATCACCCGGATACGCTTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATc                                                <  1:919905/234‑1 (MQ=255)
                                                                                                                                                                                          gggTGGGGGAGCAACGACCATCGGCGGGCTGATTTTCCGGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGTTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATc                                                <  1:641685/239‑1 (MQ=255)
                                                                                                                                                                                             tGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTAGCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCg                                              <  1:467147/238‑1 (MQ=255)
                                                                                                                                                                                             aggggTAGCAACGACCCTCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCg                                              <  1:8276/237‑1 (MQ=255)
                                                                                                                                                                                                ggTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTAtct                                                                                                                                                                                                           <  1:310192/80‑1 (MQ=255)
                                                                                                                                                                                                      aaCGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTTTCGGGTTGATTGGGTATCACCCCGATACGCGTCTCTCGCTGTTTTTTGGGTTTGCGTGGGTTGTTGTGCGGTTGTTTGGCTGGTTGTTTAAAACCCCCCCCGGTCTGCTGTTGGTCGGAAACCCGAAATCCTTCAGTTTGCTAAACCGGTGTCGGGGGAAAACCTGACCCAGCCATTCCGCAAGATATCccccccgcccc                                     >  2:692217/1‑233 (MQ=255)
                                                                                                                                                                                                         cgccAGCGGCGGGCTGATTTTCCTGCTTTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTTTCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcct                                 <  1:623091/237‑1 (MQ=255)
                                                                                                                                                                                                          aCCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGGTGATTGGTTATCACCCGGCTACGCGTATCTCGCTGTATGTCGGTTTTGCGTGGGTTGTTTTGTTGTTGTTTTGCTGGATGTTTAAACGCCCCCCCGATCGTCAACTTGCTGTAAACCCAGAAACCTTTGTTTTCTTTACCCGGTTTCCGGGGGAATCCCTGCTCGGCCTTTTCCCAAGTTATCCAccccccccccccc                                 >  2:8276/1‑230 (MQ=255)
                                                                                                                                                                                                           tctccggggggttttttttGCTGTTTTTTATTTTGGGGGTGATTGGTTTTCCCCCGGATTCGCTGTTCGTGGTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                               <  2:159248/224‑1 (MQ=255)
                                                                                                                                                                                                           ccATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAt                                                                                                                                                                                  <  1:332152/94‑1 (MQ=255)
                                                                                                                                                                                                                     gttGATTTTCCTGCTTTTTATTTTCGGGTTTATTGGTTATCACCCGGATACGCGTATCGCGCTGTAGTTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                      <  1:172896/236‑1 (MQ=255)
                                                                                                                                                                                                                     gCTGATTTTCCTGCTTTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTTGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                      <  1:61972/238‑1 (MQ=255)
                                                                                                                                                                                                                         aTTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGGTATCCGCCCCCATCCTCCTCCCCCAAAAATCCTTc                 >  2:682852/1‑239 (MQ=255)
                                                                                                                                                                                                                           tttCCTGTTTTTTTTTTTCGGGTTGTTTGGTTTTCCCCCGGTTTCGGTTATTGCGTTGTATGTCGGTTTCGCGGGGATTGTTGGGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAg               <  2:116120/239‑1 (MQ=255)
                                                                                                                                                                                                                               cTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAACCGCCGCCACGATCGTCAGCTGGCTGCAAACCAGTAATCCTTGCGTTCGGTAAGCCGGTCTAGGGGGAAATCCTGCTCCGGCTATTCCGAAAGTTAACGGCCCCCGTCCTCCACCCCCAAAAACCCTCCaaaaaa           >  2:9042/1‑234 (MQ=255)
                                                                                                                                                                                                                                   tctTTATTTGCGGGTTGTTTGGTTATCCGCGCGGTGCGGGTATCTCGTGGTTTGTCTGTTGGGGGTGGATTGTTTGGTGTTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  2:529868/239‑1 (MQ=255)
                                                                                                                                                                                                                                    cttttAGTTTCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  1:827050/237‑1 (MQ=255)
                                                                                                                                                                                                                                        tattatCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTg     >  1:116120/1‑237 (MQ=255)
                                                                                                                                                                                                                                        tattatCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa    >  2:761456/1‑238 (MQ=255)
                                                                                                                                                                                                                                        tattatCGGGTTGATTGGGTATCACCCGGATATCCGTGTCTCGCTGTATGTCGGTTTCGCGTGGCTTGTTGTGCTGTTGATTTGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGGTGAAAACCAGTAACCCTTGCTGTTTGCAAGCCGGTTCCGGTGGAAACCCTGATCCGGCCATTGGGCACGCTTTCCCCCCCCGTCACCCCCCCCCCAAACCCCTTCAgatacatacacg    >  2:555753/1‑229 (MQ=255)
                                                                                                                                                                                                                                         attatCGGGTTGATTGGTTATCACCCGGATCCGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc  <  1:107425/239‑1 (MQ=255)
                                                                                                                                                                                                                                         attatCGGGTTGATTGGTTATCACCCGGATACGCTTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc  <  1:517412/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  NC_000913/421996‑422466

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: