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

GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  NC_000913/421996‑422466
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
gttgaCCACCGCGTGTCTGTTCGGGGTGTATCTTAACCACATCATGCCGGAAAACTTCTTCCG‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATc                                                                                                                                                                                                                                             <  2:457881/236‑1 (MQ=255)
 gtggcttccggggtgttttttGCGGTGTTCTTGAACTCCATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCg                                                                                                                                                                                                                                            <  1:388554/220‑1 (MQ=255)
   gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTCAGGCGCTGAAATTTAAAGTGCCGGGTGGGTGAGCAACGACCATCGGCGGGCCGCTTTTCTCGCTCTTTATTTTCcggg                                                                                                                                                                                                                                         >  2:477776/1‑238 (MQ=255)
        ccGCGCTGCTGTTTGAGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACATTCGACACGGCGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTACGTCGCCGGTTGCCGCCTGACGAGGTTAAGGAGCTGAAATTTAAATTGACGGGGGGGGGCGCACAGCCACTCGGCCGGCGGTGCTTTCTGCTCTTTTTTATTggagggaga                                                                                                                                                                                                                                    >  1:762078/1‑232 (MQ=255)
         cgcgCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGAAGCAACGACCATCGCCGGGCGGATTTTTCTGCCCCTTATTTCTGGgttggttg                                                                                                                                                                                                                                   >  2:850064/1‑239 (MQ=255)
             ctgctgTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCGGCCAGAAGAAGTTAAGGGGCTGAAATTTAAAGTGCCGGGGGGGGAAGCACCGACCATCGGCGGGGGGTTTTTCTTCCTTTGTTTTTTAGGGTGGATGGGttg                                                                                                                                                                                                                               >  1:693174/1‑238 (MQ=255)
                 tttttGGGGTGTATCTGAACCTAATCATTGCGGAAAACGGCTTCCT‑GGTGATCGCTTCGGTGGCAACCTTCGCCACGGGGTGGGTGTGGATTATGATTCTGGTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc                                                                                                                                                                                                                           <  1:212803/237‑1 (MQ=255)
                 tgTTGGCGGTTTTTCTGATCTTGAGCATCCCGGAAATCGTCTTCCT‑GGTGATGGCTGCGCTGGCAACCTGCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCAc                                                                                                                                                                                                                           <  1:502702/239‑1 (MQ=255)
                  gTTTGTGGTGTGTTTGAACCACACTAGGCCGGAAAACGTCTTCCT‑GGGGATGGCTTCGCGGGCAACCTTTGCCAGGGTTTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCATTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTAGTATCGGGTTGATTGGTTATCAcc                                                                                                                                                                                                                          <  1:258533/239‑1 (MQ=255)
                     tGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCCGAGTTTTCTGCCCTTTATTATCGGGTTTGTTGGGTTTCAccccg                                                                                                                                                                                                                       >  2:485413/1‑237 (MQ=255)
                      ggtgggtttttgaatttactcatGCCGGCAATCGTTTTCCTGGGTGCT‑GCTTCGTTGGCAACCTTCGCCACGGTGGGGGTTGGGATTAGGGTCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGa                                                                                                                                                                                                                      <  1:500752/220‑1 (MQ=255)
                       acaggtATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGa                                                                                                                                                                                                                      <  2:77233/234‑1 (MQ=255)
                          tgtATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTTGCCTTTATTTTCCGGTTGGTTGGGTTTTAccccggtaac                                                                                                                                                                                                                  >  2:295009/1‑232 (MQ=255)
                            tATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATAcgc                                                                                                                                                                                                                 <  2:856694/238‑1 (MQ=255)
                               cTGAACTGCATCATGCGGGAAAACGGCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTAt                                                                                                                                                                                                             <  2:897848/239‑1 (MQ=255)
                                   acttcctcttGCCGGAAATCTTTTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                         <  1:258795/230‑1 (MQ=255)
                                     tACATCATGCCGGAAAACGTTCTCTT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCt                                                                                                                                                                                                       <  2:882051/239‑1 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGACGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCCGATTTTCTTTACCTTTTTTTTTCGGGTGGTTCGTTTTTCCCCGGAAACACGTATTCCGCTg                                                                                                                                                                                                      >  1:243630/1‑239 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAGGTGCCGGTGGGGGTAGCAACGACCAACGGCGTGCGTGTTTGCTTTCTCTTTTTTATTGGGGTGATTTGTTATCAACCCATTACGCGGTtctccgtt                                                                                                                                                                                                      >  1:185694/1‑235 (MQ=255)
                                      aCATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTTCTTCTCTTTTTTTTCGGGTTGATTGGTTTTCACCCGGGTACGCGTAtttc                                                                                                                                                                                                          >  2:258795/1‑235 (MQ=255)
                                        atcatGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCCCTTTTTTATCCGGTTGTTTGGTTTTAACCCGGGTAACCGTATTTCGCTGTa                                                                                                                                                                                                    >  2:273707/1‑239 (MQ=255)
                                           atGCCGGAAAACGTTTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGt                                                                                                                                                                                                 <  2:280745/239‑1 (MQ=255)
                                            tGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCGTATTTTTCTTCTCTTTTTTTTTGGGGTGGTTGTGTTTCCCCCGGGTaacactactttcctcttattg                                                                                                                                                                                                >  1:89666/1‑219 (MQ=255)
                                                    accggtttCCT‑GGTGATCGCTTCCCTTGCAACCGTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  2:564820/233‑1 (MQ=255)
                                                    aaCGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                        <  2:693174/239‑1 (MQ=255)
                                                     atgTCTTCCG‑TGTTTTCGGTTGGCTGGACTCCTTCGCCAAGGTGTGGGTGGGGATTAGGAGTGTGGTTG‑TGTAAATTGCCTTCCGTCGCCGGTTGCCGTCAGGAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTGATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgcg                                                                                                                                                                                       <  1:271338/237‑1 (MQ=255)
                                                            ccT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGAttg                                                                                                                                                                                <  2:518474/239‑1 (MQ=255)
                                                              t‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGCCGGTTTCGCGTGGAttgtt                                                                                                                                                                              <  2:617809/239‑1 (MQ=255)
                                                                   ttttgcttcgttgtccgttttCGCCGTGGTGTGGGTGTGGATTAGTATCTTTCTAT‑CTGAAATTGCACTCCGTCGCCGTTTGACCCCAGAAGAAGTTAAGGGGCTAAAATTTAAATTGCCGGGTGGGGTAGCAATGATCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTTTCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                          <  1:548102/220‑1 (MQ=255)
                                                                    aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                          <  2:770978/238‑1 (MQ=255)
                                                                         ttcgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGa                                                                                                                                                                    <  2:788893/239‑1 (MQ=255)
                                                                                      ttCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTAt                                                                                                                                                                                                                              <  2:281168/168‑1 (MQ=255)
                                                                                        cGCCACGGTGTGGGTGTGGATTATGATCCTGCTGG‑CGCAAATTGCTTTCCGTCGCCGTTTTCCGCCAGAAGCAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCACCGACCATCGGCGGGCTGAGTTTCCTTCTCTTTATTTTCGCGCTGTATGGCTATGCCCCGCTGCTCCGTGACTAGCTGCATGTTGGGTCCGCGCGGGTTGTTGTGcgtgtgactggaggagtgtggt                                                                                                                                                     >  1:845476/1‑218 (MQ=255)
                                                                                          ccgcGGTGTGGGTTGGTTTTTTGTTCCTGCTTT‑AGAAAATGTCCTCCCGGCGCCGTTTGCCGCCGGAAGAAGTTAGGGCGATTAAATTTAAAGGGCCGGGGGGGGTAGCAACGACCGTCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTaaa                                                                                                                                                   <  1:107155/236‑1 (MQ=255)
                                                                                                 gtgtgtgt‑gaTTATGATCCTCGTGC‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGTCTGATTTCCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                           <  2:768118/239‑1 (MQ=255)
                                                                                                  ggggtgttgATTATGAGCCTG‑TGTGCGCAAATTGCCTGCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGGGCTGGAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGTTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                           <  1:289594/238‑1 (MQ=255)
                                                                                                          tgtttagggacctggttgccgaaaTTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTTGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCGCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGt                                                                                                                                    <  2:752973/218‑1 (MQ=255)
                                                                                                                     tgctgT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGTGATTTTTCTGCTCTTTTTTATCGGGTGGTTTTGTTTTCACCCGGGTACGCGGATCTCGCAGTATGTCGGTTTTGCGTGGGTGTTTGTGTTGTTGGTTTGCTGGCGGTTTAAAAGCCGCCCCGACGGGCACCTGGCGGaa                                                                                                                        >  1:711465/1‑239 (MQ=255)
                                                                                                                     tgctgT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGaa                                                                                                                        <  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)
                                                                                                                                                        gccAAAAGAAGGTAAGGCGTTTAAATTTAAAGTGGCGGGTGGGGTTGCCACGCACCTGGGCGGGATGATTTTCCGGCTCTTTATTTACGGTTTTATTGGTTAACCCCCGGATACGCGTATCTCGCTGTTTGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCa                                                                                     <  1:279098/239‑1 (MQ=255)
                                                                                                                                                              agaagTTAAGGCGTTGAAATTTAAAGTTCCGGGGGGGGTAGCAACGACCATCGGCGGGTTGATTTTCCTGTTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCaggagg                                                                                <  2:291740/238‑1 (MQ=255)
                                                                                                                                                                    taagggctttaaaTTTAAAGGGCCGGGTGGGGGAGAAACGTCCCGGGGGGGGGTTGTTTTCCTGCTCTTTATTTGCGGGTCGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATc                                                                           <  1:774589/227‑1 (MQ=255)
                                                                                                                                                                                         gCCGGGGGGGGTAGCAACGGCCATCGGGGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCttacagaacgca                                                     <  2:89666/238‑12 (MQ=255)
                                                                                                                                                                                           cGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTGCTCTTTTTTATCGGGTTGATTGGTTATTACCCCGGAACGCGTATCTCGCCGTATGTCGGGTTTGCGGGGGTTGTTTTGCTGTTTGTTGGCTGGGGGTTTAAAAACCGCCACCATCGCCAGCCGGCTGAAAAACCGAAATCCTTGGGGTCTTGAAACCGGTTTCGGGGGGAAACCCT‑ATACGCGTA‑TTCCGCAAAttt                                                  >  1:784641/1‑237 (MQ=255)
                                                                                                                                                                                                  ggTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTGCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGccc                                           <  2:189523/238‑1 (MQ=255)
                                                                                                                                                                                                     ggCAACGCCCATTGGCGGGCTTTTTTTCCTTGTCGTTTTTGTCGGTTTTATTGGTTATCGCCCGGTTACGCGTATCCCGCTGTATGTCGGTCTCGCGTGGATTGTTGTGCTGTTGATGGGCGTGATGTTTAACCGCCGCCCCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCg                                        <  1:202720/237‑1 (MQ=255)
                                                                                                                                                                                                                      ggcttgttttCCTGCTCTTTATTATCGGGTTTATTGGGTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTTGCGTGGGTTTTTGTTCTGTTGATTGGCCGGCTGTTTAAACGCCCCCCCCACCGTCAGCTGGCCGGAAACCGGAAATCCCTCCGTTCTGTAAGCCGG‑TTCCAGAGGGAATACTCGTCCGGCTAGTCCCAAAAGTTTCCCCCCCCCCCCCCCCCCCCCCaaaa                        >  1:712534/7‑237 (MQ=255)
                                                                                                                                                                                                                      ggCTGTTTTTCCTGTGTTTTTTTTGAGGGTTGATTGGTTATCACCCGGATCGGCGTCTGTCGTTGTATGTCGGTTTCGCGGCGGTTGTTGGACTGTTGATGTGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                      <  1:681754/239‑1 (MQ=255)
                                                                                                                                                                                                                          ttttttCCTTTTTTTTTTTTGCGGGTTTGTTGGTTACCCCCCGGTTACGCGTATCTCGCTGTATGTCGGTTTCCGGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCtt                  <  1:850064/237‑1 (MQ=255)
                                                                                                                                                                                                                          gATTTTCCTGCTTTTTATGATCGGGTTGATTGGTTATCACCCGGATCCGCGTTTCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCtt                  <  2:276106/239‑1 (MQ=255)
                                                                                                                                                                                                                               tCCTGCTCTTTATTATCGTATTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgat             <  2:185694/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tctTTATTATCGGGTTGATTGGTTATCCCCCGGATTCGCGTATCTCGCTGGATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  2:486223/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tctTGATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg       <  2:69090/239‑1 (MQ=255)
                                                                                                                                                                                                                                       tttATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGAtgatgatgg     >  2:131964/1‑236 (MQ=255)
                                                                                                                                                                                                                                          agtatCGGGTTGAATGGTTAGCACCCGGATACGCGTATCCCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc  <  2:767282/237‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
GATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCT‑GGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGT‑CGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTGCTCTTTATTATCGGGTTGATTGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGG‑TTCAGGAGGAAATCCTGATCCGGCTA‑TTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATC  >  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: ‑