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

GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  NC_000913/421996‑422466
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
tgtgccgtccggggtgttgTTTGCGGGGTTTTTGAACTACATCATGCCGGAAAACGGCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCCCGGTGTGGGTGTGGATTATGAGCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATc                                                                                                                                                                                                                                              <  2:862129/221‑1 (MQ=255)
gatgaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATc                                                                                                                                                                                                                                              <  1:409601/238‑1 (MQ=255)
      tACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCGGGTTTTCCTGCTC‑TTTATTTTTGGGtgtg                                                                                                                                                                                                                                       >  1:159248/1‑236 (MQ=255)
       tccggggttctgTTTGGGGTGTTACCGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGAt                                                                                                                                                                                                                                      <  1:396602/230‑1 (MQ=255)
       tccGCGCTTTTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTTTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGAt                                                                                                                                                                                                                                      <  2:334528/238‑1 (MQ=255)
       tccGCGCGGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGGGTTCGCTTCGCTGGGAACTTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGAt                                                                                                                                                                                                                                      <  1:871150/238‑1 (MQ=255)
        ccGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTCGTC‑TTTTTTATCgctttttct                                                                                                                                                                                                                                     >  1:73032/1‑231 (MQ=255)
                       acaggtttCTTAACTTCATCTTGCCGGAAACCGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGa                                                                                                                                                                                                                      <  2:335202/231‑1 (MQ=255)
                                  aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTTCCGCTC‑TTTTTTATTGTGTTGATTGTG‑TTTCACCCGGGTACGCGGAttt                                                                                                                                                                                                           >  1:358183/1‑236 (MQ=255)
                                   ggtatgtttataagggacaggccTTCCTGGTTATCCCTTCGCTGGCAACCTTCGCCACGGTGGGGGTGTGGATTATGATCCTGATGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  1:379289/217‑1 (MQ=255)
                                   acccccctctGCCGGGAAACGTTTTCTTTGGTATCTCTTCGCTGGCAACCTTCGCCACGGTGGGGGGGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  2:470361/230‑1 (MQ=255)
                                   accaactcctGCCGGAAAACGTTTTCTTGGTTATCGCTTCGCTGGCAACCTTCGCCACGGGGTGGGTTTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  2:43225/230‑1 (MQ=255)
                                    cTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTCCTGCTC‑TTTTTTATCCGGTTGTTTGGG‑TTTCTCCCGGAGATGCGCAtctc                                                                                                                                                                                                          >  1:62171/1‑237 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAGGTTAAGGCGCTGAAATTTAAAGTGCCGGGCGGGGCAGCAAGGACCATAGGCGGGGCGTGTTTTCTGGTCTTTTTTTATCGGGTT‑ATTGGTGTACATCCCGGATACGCGTAtctc                                                                                                                                                                                                          >  2:172896/1‑236 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTGCTC‑TTTATTATCGGGGTGATTGGT‑TATTACCCGGATACGCGGAttcctttt                                                                                                                                                                                                      >  1:335202/1‑232 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGGGGGGTAGAACCGACCATCGGCGGGGTGGTTTTCCTGCGC‑TTTTTTTTTGGGGTGATTTGT‑TTTCACCCGGGTACGCGGATATCGccg                                                                                                                                                                                                      >  2:919905/1‑237 (MQ=255)
                                            tGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGGGATTTTCCTTCTC‑TTTTTTTTCGGGGTTATTTGT‑TTTCACCCCGTTAccgggttttcgctttattgc                                                                                                                                                                                                >  1:904204/1‑220 (MQ=255)
                                                    aaCTTCTTCCTGGTGGTTGTTTTGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTCTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  2:358183/239‑1 (MQ=255)
                                                    aaCTTCTTCCTGGTGATCGCCTCGCTGGCAACCTTCGCCACGGTGTTGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  1:24447/239‑1 (MQ=255)
                                                       gTCTTCCTGGTGATCGCTGCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTAGGATCCTGCTGTCGCAAATTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTg                                                                                                                                                                                     <  1:251731/239‑1 (MQ=255)
                                                        tCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTg                                                                                                                                                                                     <  2:440203/238‑1 (MQ=255)
                                                         cTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGa                                                                                                                                                                                   <  1:741344/239‑1 (MQ=255)
                                                            ccTGGTTATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttg                                                                                                                                                                                <  1:464870/239‑1 (MQ=255)
                                                            ccTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTTTTATCGGGTTGATTTGT‑TATCACCCGGATACGCGTATTTCCCTGTATGTTGGTTTCTCGGGGGttt                                                                                                                                                                                >  1:607044/1‑238 (MQ=255)
                                                              tggttgtcgctctgtgggCAACCTTCGCCACGGTGTGGGTTTGGATTATGTGCCTGCTTTCGCAAATTGCCTTCCGGCGCCGGTGGCCGGCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttgtt                                                                                                                                                                              <  2:401817/223‑1 (MQ=255)
                                                              tGGTGTTCGCTTCGCTGGCAACCTTCGCCTCGGTGTGGGTTTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttgtt                                                                                                                                                                              <  2:45501/239‑1 (MQ=255)
                                                                  tttcgcttcgttGGCAACCTTCGCCCCGGGGGGGGGGTGGATTATGATCCGGCTGGCGCAAATTGCCTTCCGTCGCCGTTTGCCGGCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                          <  2:468412/237‑1 (MQ=255)
                                                                   aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGGTACGCGTATTTCGCTGTATGTCGGTTTTGGGTGGGTTGTTTTGCt                                                                                                                                                                         >  1:45501/1‑239 (MQ=255)
                                                                          cgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCggg                                                                                                                                                                                                                                           <  1:541080/167‑1 (MQ=255)
                                                                                 aaCCTTCGCCACGGTGTGGGTGTGGATTATGGCCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTgg                                                                                                                                                           <  1:219375/238‑1 (MQ=255)
                                                                                         ccACGGGGTGGGTGTGGATTTTGATCCTGCTGTCGCAAATGGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTaaa                                                                                                                                                  <  1:364660/239‑1 (MQ=255)
                                                                                              gtgtgggtgtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCCTCCC‑TTTTTTATCGGGGTGGTTGGT‑TTTTACCCGGGTAGACCGATTTCCCTGTATGTCGGGTTTTGGTGGGTTTTTGTGTGG‑TTGtttttgtggggggtttaaaagc                                                                                                                                              >  2:438579/1‑219 (MQ=255)
                                                                                               tgtgggtgtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTAt                                                                                                                                                                                                             <  1:872126/175‑1 (MQ=255)
                                                                                               tgtgggtgtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTg                                                                                                                                                            <  1:115368/223‑1 (MQ=255)
                                                                                                  ggttgtggATTTTGATCCGGCGTTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGGGGTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATGGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                          <  2:766446/235‑1 (MQ=255)
                                                                                                  gggtgtggATTATGATCCTGCTGTCGAAAATTGCCTTCCGTCGGCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                          <  1:438579/238‑1 (MQ=255)
                                                                                                    gtttgtattgtttttctggttgcGCAATTTGCCTTCCGGCGCCGTTGGCCCCCAGAAGAAGTTAAGGGCTTTAAATTTAAAGGGCCGGGTGGGGTAGCAACGAGCATCGGCGGGCTGATTTTCCTGCTC‑TTTATGATCGGGTTGATTGGT‑TATCTCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACg                                                                                                                                        <  2:436948/216‑1 (MQ=255)
                                                                                                    gggtggATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACg                                                                                                                                        <  1:26787/236‑1 (MQ=255)
                                                                                                         gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGGGATTTTCCTGCTC‑TTTATTTTCGGGGTGGTTGGT‑TTTACCCCCGATATCCGGTCCTCCGCGTATGTCGGTTTGTCGTGGGTTTGTGTGCTG‑TTTTTTGGTGGGGTGTTTTAACGCCCCCCCGCTCggg                                                                                                                                  >  2:716029/1‑237 (MQ=255)
                                                                                                         gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGTGTTTTCTGCTC‑TTTATTTTGTGGTTGATTGGT‑TTTCACCCGGGTACGCGTATCTCGCTGTTTGTCGGGTTCGCGTGGGGTGTTTTGGTG‑TTGTTTTGGTGGGGGTTTAAACGCCCCCCCGAGCgcc                                                                                                                                  >  2:691401/1‑237 (MQ=255)
                                                                                                         gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGGCTCTCGCTGTATGTTGGGTTCGTGTGGGTTGTTTTGTTG‑TTGGTTGGCGGGGTGTTTaaaactgtatcttttaaaa                                                                                                                                  >  2:115368/1‑223 (MQ=255)
                                                                                                             aggactctgctgTCTCAAATTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCGC‑TTTATTATCGGGTTGATTGGT‑TATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAgct                                                                                                                              <  1:54631/233‑1 (MQ=255)
                                                                                                                    tggtgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑aa                                                                                                                      <  1:123786/236‑1 (MQ=255)
                                                                                                                    tgctgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGGTGATTGGT‑TATCACCCGGATACGCGTATCTCGCCTTATGTCCGTTTTGCGTGGGTTGTTGTGCTG‑TTGGTTGGCTTGGTGTTTTAACGCCGCCACGGTCGTCACCTTGCTG‑aa                                                                                                                      >  1:808258/1‑239 (MQ=255)
                                                                                                                      ctgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGTGAGTTTCATGCGC‑TTTATTATTGGGGTGATTGGT‑TATGACCCGGGAACGCGTATCTCCCTTGATGTCGGGTTTGCGTGGGGTGTTGTGGTG‑TTGATTGGGTGGGGTTTTAAAAGCCCCCACGagaggcgggcgggggaaaa                                                                                                                     >  2:517412/1‑221 (MQ=255)
                                                                                                                                               gTTTGCCGGCAGAAGAAGTTAAGGCGTTGAAATTTAAAGGGGCGGGTGGGGTAGCAACGACCATCGGGGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCGCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGc                                                                                           <  1:285816/239‑1 (MQ=255)
                                                                                                                                                tttGCCGCCAGAAGAAGTTAGGGGGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGcc                                                                                          <  1:784581/239‑1 (MQ=255)
                                                                                                                                                      gccAGAAGAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCa                                                                                    <  1:691401/239‑1 (MQ=255)
                                                                                                                                                      gccAAAAGAATTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCCTCGGCGGGCTGATTTTCCTGTTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCa                                                                                    <  1:896274/239‑1 (MQ=255)
                                                                                                                                                                    aggtgtttaaaTTTAAAGTGCGGGGTGGGGTGGCAACGACCTTGGGCGGGCTGTTTTTTCTGATC‑TTTATTTTGGTGTTTATTGTT‑TATCCCCCGGATTCGCGTCTCGTGCTGTAGGTCTGTTTCGGGTGGATTGTTGTGCTG‑TTGATGGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGa                                                                      <  2:808258/231‑1 (MQ=255)
                                                                                                                                                                          tGAAATTTAAGGTGCCGGGTTGGGTAGCAAGGACCATCGGGGGTCTGATTTTCCTGCTT‑TTTATTATCGGGTTTATTGGT‑TATCACCCGGATACGCTTATCTCGTTGTTTGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGc                                                                <  2:738237/239‑1 (MQ=255)
                                                                                                                                                                                       gCCGGGGTGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTT‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAg                                                     <  1:544255/237‑1 (MQ=255)
                                                                                                                                                                                          ggggtgggTAGCACCGACCAGCGGCGGGCTTATTTTCCTGCTT‑TTTATTATCGGGTTTATTGGT‑TATCACCCGGATACGCTTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATc                                                <  1:919905/234‑1 (MQ=255)
                                                                                                                                                                                          gggTGGGGGAGCAACGACCATCGGCGGGCTGATTTTCCGGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGTTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATc                                                <  1:641685/239‑1 (MQ=255)
                                                                                                                                                                                             tGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TAGCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCg                                              <  1:467147/238‑1 (MQ=255)
                                                                                                                                                                                             aggggTAGCAACGACCCTCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCg                                              <  1:8276/237‑1 (MQ=255)
                                                                                                                                                                                                ggTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTAtct                                                                                                                                                                                                           <  1:310192/80‑1 (MQ=255)
                                                                                                                                                                                                      aaCGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTTTCGGGTTGATTGGG‑TATCACCCCGATACGCGTCTCTCGCTGTTTTTTGGGTTTGCGTGGGTTGTTGTGCGG‑TTGTTTGGCTGGTTGTTTAAAACCCCCCCCGGTCTGCTGTTGGTCG‑GAAACCCGAAATCCTTCAGTTTGCTAAACCGGTGTCGGGGGAAAACCTGACCCAGCCATTCCGCAAGATATCccccccgcccc                                     >  2:692217/1‑233 (MQ=255)
                                                                                                                                                                                                         cgccAGCGGCGGGCTGATTTTCCTGCTT‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTTTCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcct                                 <  1:623091/237‑1 (MQ=255)
                                                                                                                                                                                                          aCCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGGTGATTGGT‑TATCACCCGGCTACGCGTATCTCGCTGTATGTCGGTTTTGCGTGGGTTGTTTTGTTG‑TTGTTTTGCTGGATGTTTAAACGCCCCCCCGATCGTCAACTTGCTGTAAACCCAGAAACCTTTGTTTTCTTTACCCGGTTTCCGGGGGAATCCCTGCTCGGCCTTTTCC‑CAAGTTATCCAccccccccccccc                                 >  2:8276/1‑230 (MQ=255)
                                                                                                                                                                                                           tctccggggggttttttttGCTGTTT‑TTTATTTTGGGGGTGATTGGT‑TTTCCCCCGGATTCGCTGTTCGTGGTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                               <  2:159248/224‑1 (MQ=255)
                                                                                                                                                                                                           ccATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAt                                                                                                                                                                                  <  1:332152/94‑1 (MQ=255)
                                                                                                                                                                                                                     gttGATTTTCCTGCTT‑TTTATTTTCGGGTTTATTGGT‑TATCACCCGGATACGCGTATCGCGCTGTAGTTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                      <  1:172896/236‑1 (MQ=255)
                                                                                                                                                                                                                     gCTGATTTTCCTGCTT‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTTGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                      <  1:61972/238‑1 (MQ=255)
                                                                                                                                                                                                                         aTTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGGTATCCGCCCCCATCCTCCTCCCCCAAAAATCCTTc                 >  2:682852/1‑239 (MQ=255)
                                                                                                                                                                                                                           tttCCTGTTT‑TTTTTTTTCGGGTTGTTTGGT‑TTTCCCCCGGTTTCGGTTATTGCGTTGTATGTCGGTTTCGCGGGGATTGTTGGGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAg               <  2:116120/239‑1 (MQ=255)
                                                                                                                                                                                                                               cTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAACCGCCGCCACGATCGTCAGCTGGCTG‑CAAACCAGTAATCCTTGCGTTCGGTAAGCCGGTCTAGGGGGAAATCCTGCTCCGGCTATTCCGAAAGTTAACGGCCCCCGTCCTCCACCCCCAAAAACCCTCCaaaaaa           >  2:9042/1‑234 (MQ=255)
                                                                                                                                                                                                                                   tc‑tTTATTTGCGGGTTGTTTGGT‑TATCCGCGCGGTGCGGGTATCTCGTGGTTTGTCTGTTGGGGGTGGATTGTT‑TGGTGTTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  2:529868/239‑1 (MQ=255)
                                                                                                                                                                                                                                    cttttAGTTTCGGGTTGATTGGT‑TATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  1:827050/237‑1 (MQ=255)
                                                                                                                                                                                                                                        tattatCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTg     >  1:116120/1‑237 (MQ=255)
                                                                                                                                                                                                                                        tattatCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa    >  2:761456/1‑238 (MQ=255)
                                                                                                                                                                                                                                        tattatCGGGTTGATTGGG‑TATCACCCGGATATCCGTGTCTCGCTGTATGTCGGTTTCGCGTGGCTTGTTGTGCTG‑TTGATTTGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGGTG‑AAAACCAGTAACCCTTGCTGTTTGCAAGCCGGTTCCGGTGGAAACCCTGATCCGGCCATTGGGCACGCTTTCCCCCCCCGTCACCCCCCCCCCAAACCCCTTCAgatacatacacg    >  2:555753/1‑229 (MQ=255)
                                                                                                                                                                                                                                         attatCGGGTTGATTGGT‑TATCACCCGGATCCGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc  <  1:107425/239‑1 (MQ=255)
                                                                                                                                                                                                                                         attatCGGGTTGATTGGT‑TATCACCCGGATACGCTTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc  <  1:517412/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTC‑TTTATTATCGGGTTGATTGGT‑TATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTG‑TTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTG‑AAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  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: ‑