Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 254,233 A→T 18.3% intergenic (+31/+26) prfH → / ← pepD pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913254,2330AT18.3% 379.2 / 75.1 163intergenic (+31/+26)prfH/pepDpseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (67/66);  new base T (12/18);  total (79/84)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 170 total aligned reads displayed.

TTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTT  >  NC_000913/254134‑254332
                                                                                                   |                                                                                                   
ttCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                     >  2:1881833/1‑100 (MQ=255)
ttCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                     >  1:1350044/1‑100 (MQ=255)
 tCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAAc                                                                                                    <  2:2811984/100‑1 (MQ=255)
  ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATGATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                                   >  2:2458743/1‑100 (MQ=255)
  ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                                   >  1:3483196/1‑100 (MQ=255)
  ccATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACg                                                                                                   <  1:1371330/100‑1 (MQ=255)
   cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                                  >  1:2644234/1‑100 (MQ=255)
   cATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGc                                                                                                  <  2:2664807/100‑1 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                              >  1:1024840/1‑100 (MQ=255)
       aCCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                              >  2:3470128/1‑100 (MQ=255)
         cAGATTGAACGTGGCGGCCTGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGGAGGTGGAAACTGCCGGAAGGCGATTAAACGc                                                                                                  >  2:1105498/1‑94 (MQ=255)
          aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCgg                                                                                           <  1:2959729/100‑1 (MQ=38)
          aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCgg                                                                                           <  1:2021556/100‑1 (MQ=255)
           gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                          <  2:38511/100‑1 (MQ=38)
           gATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                          <  1:1592866/100‑1 (MQ=255)
            aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgca                                                                                         <  1:3340450/100‑1 (MQ=38)
            aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgca                                                                                         <  1:822508/100‑1 (MQ=38)
            aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgc                                                                                          <  1:512665/99‑1 (MQ=38)
            aTTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgca                                                                                         <  2:647182/100‑1 (MQ=255)
             ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcag                                                                                        <  2:1048424/100‑1 (MQ=38)
              tGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagc                                                                                       <  1:860672/100‑1 (MQ=255)
                aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                     <  1:3570603/100‑1 (MQ=255)
                  cGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaat                                                                                   >  2:2413090/1‑100 (MQ=255)
                    tGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGAGTAAACGCCATCCGGCAGcaaatca                                                                                 >  2:1731521/1‑100 (MQ=255)
                     ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTCCGgcagtttcatcc                                                                                <  1:1761866/100‑9 (MQ=21)
                     ggCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa                                                                                <  1:772232/100‑1 (MQ=255)
                      gCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa                                                                               <  2:895661/100‑1 (MQ=255)
                      gCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa                                                                               <  1:341313/100‑1 (MQ=255)
                      gCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa                                                                               <  2:3185658/100‑1 (MQ=255)
                        aaCCCGCGACGCACATTTACAGGGATGGCTATTATCGAATGATAATGAAGGATGAAACTGCCATATGGCGTTTAATCGCCTTCCGgcagtttcatc                                                                                 <  2:150263/96‑8 (MQ=21)
                           ccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAAtt                                                                          <  1:2898051/100‑1 (MQ=255)
                           ccGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAAtt                                                                          <  1:1133674/100‑1 (MQ=255)
                            cgcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                                         >  2:2579489/1‑100 (MQ=255)
                            cgcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                                         >  1:96077/1‑100 (MQ=255)
                               gACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACtt                                                                      <  2:210198/100‑1 (MQ=255)
                                aCGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTc                                                                     >  2:3608726/1‑100 (MQ=255)
                                   gACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGcc                                                                  >  1:3205662/1‑100 (MQ=255)
                                    aCATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCg                                                                 <  2:1261183/100‑1 (MQ=255)
                                    aCATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCg                                                                 <  2:3579814/100‑1 (MQ=255)
                                     cATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCgg                                                                >  2:1101718/1‑100 (MQ=255)
                                     cATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCgg                                                                <  1:1166188/100‑1 (MQ=255)
                                         tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtt                                                            >  1:787874/1‑100 (MQ=255)
                                           cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttc                                                          >  2:1048508/1‑100 (MQ=255)
                                           cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttc                                                          >  2:2789737/1‑100 (MQ=255)
                                                aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTca                                                     >  1:1680931/1‑100 (MQ=255)
                                                  ggCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                                   <  1:3501557/100‑1 (MQ=255)
                                                    cTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                                 >  1:3489012/1‑100 (MQ=255)
                                                     ttttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt                                                <  2:1153716/100‑1 (MQ=255)
                                                       ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                              <  2:2098641/100‑1 (MQ=255)
                                                        tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                                             >  1:1721547/1‑100 (MQ=255)
                                                         atccttcgATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                            >  2:434685/8‑100 (MQ=255)
                                                         aTCGAAGTATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                            <  1:3168661/100‑1 (MQ=255)
                                                         aTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                                            <  1:2726625/100‑1 (MQ=255)
                                                                 aTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgt                                    <  1:16509/100‑1 (MQ=255)
                                                                    aTGAAGGATGAAACTGCGGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTAAGCAGTTCAGTCAGCAgtgt                                  >  2:2013205/1‑99 (MQ=255)
                                                                    aTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                                 >  1:3130888/1‑100 (MQ=255)
                                                                    aTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCAGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTTAGTCAGCAGTGTc                                 >  2:3013944/1‑100 (MQ=255)
                                                                      gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                               >  2:580195/1‑100 (MQ=255)
                                                                      gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                               >  2:2481702/1‑100 (MQ=255)
                                                                       atttgatttgcTGCCGGATGGCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                              <  2:3331015/90‑1 (MQ=21)
                                                                       aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg                              <  2:3205637/100‑1 (MQ=255)
                                                                        aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt                             >  1:221550/1‑100 (MQ=255)
                                                                         ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                            >  1:2423736/1‑100 (MQ=255)
                                                                          gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTaa                           <  1:28642/100‑1 (MQ=255)
                                                                           aTGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt                          >  1:1052484/1‑100 (MQ=255)
                                                                               aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc                      <  2:2347290/100‑1 (MQ=255)
                                                                                aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt                     <  1:2765/100‑1 (MQ=255)
                                                                                 cTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa                    >  1:2339014/1‑100 (MQ=38)
                                                                                 cTGCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCAAGTAATGACCt                     >  1:1209099/1‑99 (MQ=37)
                                                                                 cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa                    <  1:606359/100‑1 (MQ=255)
                                                                                   gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCGGCAGTTCAGTCAGCAGTGTCCAGTAATGACCtc                    >  2:2410154/1‑97 (MQ=37)
                                                                                   gCCGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg                  >  2:1553879/1‑100 (MQ=38)
                                                                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg                  >  1:1027385/1‑100 (MQ=255)
                                                                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg                  <  1:580197/100‑1 (MQ=255)
                                                                                    ccGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc                 >  1:1828063/1‑100 (MQ=38)
                                                                                    ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc                 >  1:1964875/1‑100 (MQ=255)
                                                                                    ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc                 <  2:571557/100‑1 (MQ=255)
                                                                                     cGGATGGCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCt                >  1:2951448/1‑100 (MQ=38)
                                                                                        aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc             <  2:1024832/100‑1 (MQ=255)
                                                                                         aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg            <  1:154149/100‑1 (MQ=255)
                                                                                         aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg            <  1:702118/100‑1 (MQ=255)
                                                                                            cgggtAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTAGAtg         <  2:1189940/96‑1 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg         >  2:789349/1‑100 (MQ=255)
                                                                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg         <  1:3487240/100‑1 (MQ=255)
                                                                                             gATTAAACTCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg         <  2:1614158/99‑1 (MQ=255)
                                                                                             gATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCGGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg         <  2:2912531/99‑1 (MQ=255)
                                                                                             gATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgt        >  1:2224578/1‑100 (MQ=255)
                                                                                             gATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgt        >  1:209387/1‑100 (MQ=255)
                                                                                             gATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgt        >  1:3350768/1‑100 (MQ=255)
                                                                                             gATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgt        >  1:1838980/1‑100 (MQ=255)
                                                                                              aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg       >  2:3499238/1‑100 (MQ=255)
                                                                                              aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg       <  1:3607282/100‑1 (MQ=255)
                                                                                              aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgtg       <  1:3386138/100‑1 (MQ=255)
                                                                                               ttAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGa      <  1:1974299/100‑1 (MQ=255)
                                                                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa     >  1:220784/1‑100 (MQ=255)
                                                                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa     >  2:2188695/1‑100 (MQ=255)
                                                                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa     <  2:787869/100‑1 (MQ=255)
                                                                                                 aaaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAAc    >  2:37627/1‑100 (MQ=255)
                                                                                                 aaaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGCCCAGTAATGACCTACGCTTTCGATGTGAAc    >  2:890123/1‑100 (MQ=255)
                                                                                                   aCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACtt  >  1:3015488/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTT  >  NC_000913/254134‑254332

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: