Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,231 A→T 18.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,2310AT18.8% 118.1 / 3.2 58I416F (ATC→TTC) proYproline‑specific permease
Reads supporting (aligned to +/- strand):  ref base A (9/36);  new base T (8/4);  total (17/41)
Fisher's exact test for biased strand distribution p-value = 3.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  NC_000913/421996‑422466
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
gttgaCCACCGCGTGTCTGTTCGGGGTGTATCTTAACCACATCATGCCGGAAAACTTCTTCCGGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATc                                                                                                                                                                                                                                             <  2:457881/236‑1 (MQ=255)
 gtggcttccggggtgttttttGCGGTGTTCTTGAACTCCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCg                                                                                                                                                                                                                                            <  1:388554/220‑1 (MQ=255)
   gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTCAGGCGCTGAAATTTAAAGTGCCGGGTGGGTGAGCAACGACCATCGGCGGGCCGCTTTTCTCGCTCTTTATTTTCcggg                                                                                                                                                                                                                                         >  2:477776/1‑238 (MQ=255)
        ccGCGCTGCTGTTTGAGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACATTCGACACGGCGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTACGTCGCCGGTTGCCGCCTGACGAGGTTAAGGAGCTGAAATTTAAATTGACGGGGGGGGGCGCACAGCCACTCGGCCGGCGGTGCTTTCTGCTCTTTTTTATTggagggaga                                                                                                                                                                                                                                    >  1:762078/1‑232 (MQ=255)
         cgcgCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGAAGCAACGACCATCGCCGGGCGGATTTTTCTGCCCCTTATTTCTGGgttggttg                                                                                                                                                                                                                                   >  2:850064/1‑239 (MQ=255)
             ctgctgTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCGGCCAGAAGAAGTTAAGGGGCTGAAATTTAAAGTGCCGGGGGGGGAAGCACCGACCATCGGCGGGGGGTTTTTCTTCCTTTGTTTTTTAGGGTGGATGGGttg                                                                                                                                                                                                                               >  1:693174/1‑238 (MQ=255)
                 tttttGGGGTGTATCTGAACCTAATCATTGCGGAAAACGGCTTCCTGGTGATCGCTTCGGTGGCAACCTTCGCCACGGGGTGGGTGTGGATTATGATTCTGGTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc                                                                                                                                                                                                                           <  1:212803/237‑1 (MQ=255)
                 tgTTGGCGGTTTTTCTGATCTTGAGCATCCCGGAAATCGTCTTCCTGGTGATGGCTGCGCTGGCAACCTGCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc                                                                                                                                                                                                                           <  1:502702/239‑1 (MQ=255)
                  gTTTGTGGTGTGTTTGAACCACACTAGGCCGGAAAACGTCTTCCTGGGGATGGCTTCGCGGGCAACCTTTGCCAGGGTTTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCATTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTAGTATCGGGTTGATTGGTTATCAcc                                                                                                                                                                                                                          <  1:258533/239‑1 (MQ=255)
                     tGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCCGAGTTTTCTGCCCTTTATTATCGGGTTTGTTGGGTTTCAccccg                                                                                                                                                                                                                       >  2:485413/1‑237 (MQ=255)
                      ggtgggtttttgaatttactcatGCCGGCAATCGTTTTCCTGGGTGCTGCTTCGTTGGCAACCTTCGCCACGGTGGGGGTTGGGATTAGGGTCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGa                                                                                                                                                                                                                      <  1:500752/220‑1 (MQ=255)
                       acaggtATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGa                                                                                                                                                                                                                      <  2:77233/234‑1 (MQ=255)
                          tgtATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTTGCCTTTATTTTCCGGTTGGTTGGGTTTTAccccggtaac                                                                                                                                                                                                                  >  2:295009/1‑232 (MQ=255)
                            tATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATAcgc                                                                                                                                                                                                                 <  2:856694/238‑1 (MQ=255)
                               cTGAACTGCATCATGCGGGAAAACGGCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTAt                                                                                                                                                                                                             <  2:897848/239‑1 (MQ=255)
                                   acttcctcttGCCGGAAATCTTTTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  1:258795/230‑1 (MQ=255)
                                     tACATCATGCCGGAAAACGTTCTCTTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCt                                                                                                                                                                                                       <  2:882051/239‑1 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGACGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCCGATTTTCTTTACCTTTTTTTTTCGGGTGGTTCGTTTTTCCCCGGAAACACGTATTCCGCTg                                                                                                                                                                                                      >  1:243630/1‑239 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAGGTGCCGGTGGGGGTAGCAACGACCAACGGCGTGCGTGTTTGCTTTCTCTTTTTTATTGGGGTGATTTGTTATCAACCCATTACGCGGTtctccgtt                                                                                                                                                                                                      >  1:185694/1‑235 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTTCTTCTCTTTTTTTTCGGGTTGATTGGTTTTCACCCGGGTACGCGTAtttc                                                                                                                                                                                                          >  2:258795/1‑235 (MQ=255)
                                        atcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCCCTTTTTTATCCGGTTGTTTGGTTTTAACCCGGGTAACCGTATTTCGCTGTa                                                                                                                                                                                                    >  2:273707/1‑239 (MQ=255)
                                           atGCCGGAAAACGTTTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGt                                                                                                                                                                                                 <  2:280745/239‑1 (MQ=255)
                                            tGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCGTATTTTTCTTCTCTTTTTTTTTGGGGTGGTTGTGTTTCCCCCGGGTaacactactttcctcttattg                                                                                                                                                                                                >  1:89666/1‑219 (MQ=255)
                                                    accggtttCCTGGTGATCGCTTCCCTTGCAACCGTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  2:564820/233‑1 (MQ=255)
                                                    aaCGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  2:693174/239‑1 (MQ=255)
                                                     atgTCTTCCGTGTTTTCGGTTGGCTGGACTCCTTCGCCAAGGTGTGGGTGGGGATTAGGAGTGTGGTTGTGTAAATTGCCTTCCGTCGCCGGTTGCCGTCAGGAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTGATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgcg                                                                                                                                                                                       <  1:271338/237‑1 (MQ=255)
                                                            ccTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttg                                                                                                                                                                                <  2:518474/239‑1 (MQ=255)
                                                              tGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGCCGGTTTCGCGTGGAttgtt                                                                                                                                                                              <  2:617809/239‑1 (MQ=255)
                                                                   ttttgcttcgttgtccgttttCGCCGTGGTGTGGGTGTGGATTAGTATCTTTCTATCTGAAATTGCACTCCGTCGCCGTTTGACCCCAGAAGAAGTTAAGGGGCTAAAATTTAAATTGCCGGGTGGGGTAGCAATGATCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTTTCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                          <  1:548102/220‑1 (MQ=255)
                                                                    aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                          <  2:770978/238‑1 (MQ=255)
                                                                         ttcgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGa                                                                                                                                                                    <  2:788893/239‑1 (MQ=255)
                                                                                      ttCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTAt                                                                                                                                                                                                                              <  2:281168/168‑1 (MQ=255)
                                                                                        cGCCACGGTGTGGGTGTGGATTATGATCCTGCTGGCGCAAATTGCTTTCCGTCGCCGTTTTCCGCCAGAAGCAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCACCGACCATCGGCGGGCTGAGTTTCCTTCTCTTTATTTTCGCGCTGTATGGCTATGCCCCGCTGCTCCGTGACTAGCTGCATGTTGGGTCCGCGCGGGTTGTTGTGcgtgtgactggaggagtgtggt                                                                                                                                                     >  1:845476/1‑218 (MQ=255)
                                                                                          ccgcGGTGTGGGTTGGTTTTTTGTTCCTGCTTTAGAAAATGTCCTCCCGGCGCCGTTTGCCGCCGGAAGAAGTTAGGGCGATTAAATTTAAAGGGCCGGGGGGGGTAGCAACGACCGTCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTaaa                                                                                                                                                   <  1:107155/236‑1 (MQ=255)
                                                                                                 gtgtgtgtgaTTATGATCCTCGTGCCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGTCTGATTTCCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                           <  2:768118/239‑1 (MQ=255)
                                                                                                  ggggtgttgATTATGAGCCTGTGTGCGCAAATTGCCTGCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGGGCTGGAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGTTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                           <  1:289594/238‑1 (MQ=255)
                                                                                                          tgtttagggacctggttgccgaaaTTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTTGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCGCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGt                                                                                                                                    <  2:752973/218‑1 (MQ=255)
                                                                                                                     tgctgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGTGATTTTTCTGCTCTTTTTTATCGGGTGGTTTTGTTTTCACCCGGGTACGCGGATCTCGCAGTATGTCGGTTTTGCGTGGGTGTTTGTGTTGTTGGTTTGCTGGCGGTTTAAAAGCCGCCCCGACGGGCACCTGGCGGaa                                                                                                                        >  1:711465/1‑239 (MQ=255)
                                                                                                                     tgctgTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGaa                                                                                                                        <  2:116414/239‑1 (MQ=255)
                                                                                                                             ggaaaTTGCCTTCCGTCGCCGTTTGCCGCCAGAAAAAGTTAAGGCGCTGAAATTTAAAGGGCCGGGTGGGGTAGCAACGACCATCGGCGGGTTTATTTTCCTGCTTTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGt                                                                                                                 <  1:634403/237‑1 (MQ=255)
                                                                                                                                       ttCCGTCTCCGTTGGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGCCCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCg                                                                                                       <  2:645410/239‑1 (MQ=255)
                                                                                                                                       ttCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCg                                                                                                       <  2:388454/239‑1 (MQ=255)
                                                                                                                                               ccGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGGTACCCGTATCTCGCCGTATTTCGGTTTCGCGTGGATTTTTGTTCTGTTGATTTGCTGGAATTTCAAAAGCCGCCACGATCGTCAGCTGGCTGCAAAACACGAATCCTTTCGTTCTggta                                                                                               >  2:17120/1‑236 (MQ=255)
                                                                                                                                                        gccAAAAGAAGGTAAGGCGTTTAAATTTAAAGTGGCGGGTGGGGTTGCCACGCACCTGGGCGGGATGATTTTCCGGCTCTTTATTTACGGTTTTATTGGTTAACCCCCGGATACGCGTATCTCGCTGTTTGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCa                                                                                     <  1:279098/239‑1 (MQ=255)
                                                                                                                                                              agaagTTAAGGCGTTGAAATTTAAAGTTCCGGGGGGGGTAGCAACGACCATCGGCGGGTTGATTTTCCTGTTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCaggagg                                                                                <  2:291740/238‑1 (MQ=255)
                                                                                                                                                                    taagggctttaaaTTTAAAGGGCCGGGTGGGGGAGAAACGTCCCGGGGGGGGGTTGTTTTCCTGCTCTTTATTTGCGGGTCGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATc                                                                           <  1:774589/227‑1 (MQ=255)
                                                                                                                                                                                         gCCGGGGGGGGTAGCAACGGCCATCGGGGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCttacagaacgca                                                     <  2:89666/238‑12 (MQ=255)
                                                                                                                                                                                           cGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTGCTCTTTTTTATCGGGTTGATTGGTTATTACCCCGGAACGCGTATCTCGCCGTATGTCGGGTTTGCGGGGGTTGTTTTGCTGTTTGTTGGCTGGGGGTTTAAAAACCGCCACCATCGCCAGCCGGCTGAAAAACCGAAATCCTTGGGGTCTTGAAACCGGTTTCGGGGGGAAACCCTATACGCGTATTCCGCAAAttt                                                  >  1:784641/1‑237 (MQ=255)
                                                                                                                                                                                                  ggTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTGCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGccc                                           <  2:189523/238‑1 (MQ=255)
                                                                                                                                                                                                     ggCAACGCCCATTGGCGGGCTTTTTTTCCTTGTCGTTTTTGTCGGTTTTATTGGTTATCGCCCGGTTACGCGTATCCCGCTGTATGTCGGTCTCGCGTGGATTGTTGTGCTGTTGATGGGCGTGATGTTTAACCGCCGCCCCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCg                                        <  1:202720/237‑1 (MQ=255)
                                                                                                                                                                                                                      ggcttgttttCCTGCTCTTTATTATCGGGTTTATTGGGTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTTGCGTGGGTTTTTGTTCTGTTGATTGGCCGGCTGTTTAAACGCCCCCCCCACCGTCAGCTGGCCGGAAACCGGAAATCCCTCCGTTCTGTAAGCCGGTTCCAGAGGGAATACTCGTCCGGCTAGTCCCAAAAGTTTCCCCCCCCCCCCCCCCCCCCCCaaaa                        >  1:712534/7‑237 (MQ=255)
                                                                                                                                                                                                                      ggCTGTTTTTCCTGTGTTTTTTTTGAGGGTTGATTGGTTATCACCCGGATCGGCGTCTGTCGTTGTATGTCGGTTTCGCGGCGGTTGTTGGACTGTTGATGTGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                      <  1:681754/239‑1 (MQ=255)
                                                                                                                                                                                                                          ttttttCCTTTTTTTTTTTTGCGGGTTTGTTGGTTACCCCCCGGTTACGCGTATCTCGCTGTATGTCGGTTTCCGGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCtt                  <  1:850064/237‑1 (MQ=255)
                                                                                                                                                                                                                          gATTTTCCTGCTTTTTATGATCGGGTTGATTGGTTATCACCCGGATCCGCGTTTCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCtt                  <  2:276106/239‑1 (MQ=255)
                                                                                                                                                                                                                               tCCTGCTCTTTATTATCGTATTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgat             <  2:185694/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tctTTATTATCGGGTTGATTGGTTATCCCCCGGATTCGCGTATCTCGCTGGATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  2:486223/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tctTGATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  2:69090/239‑1 (MQ=255)
                                                                                                                                                                                                                                       tttATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGAtgatgatgg     >  2:131964/1‑236 (MQ=255)
                                                                                                                                                                                                                                          agtatCGGGTTGAATGGTTAGCACCCGGATACGCGTATCCCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc  <  2:767282/237‑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 < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: