Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,213 C→T 17.1% T121I (ACC→ATC)  rluB → 23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,327,2130CT17.1% 270.0 / 43.0 119T121I (ACC→ATC) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base C (56/42);  new base T (11/9);  total (67/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 130 total aligned reads displayed.

GTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGC  >  NC_000913/1327114‑1327311
                                                                                                   |                                                                                                  
gTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTcaccac                                                                                                    >  2:734038/1‑100 (MQ=255)
 tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                   >  2:202410/1‑100 (MQ=255)
 tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                   <  1:340371/100‑1 (MQ=255)
 tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                   >  1:1992406/1‑100 (MQ=255)
 tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                   <  2:1053359/100‑1 (MQ=255)
   cAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                 >  2:1347724/1‑100 (MQ=255)
    aaCGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                >  2:1195822/1‑100 (MQ=255)
    aaCGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                >  1:2985226/1‑100 (MQ=255)
     aCGGTGTTTGTCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTTATACCTGTGGTCTGCTGCTGTTCACCACCGATg                                                                                               <  2:2505030/100‑1 (MQ=255)
      cGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTgg                                                                                              <  1:2176151/100‑1 (MQ=255)
       ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                              <  2:2097758/99‑1 (MQ=255)
       ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                             >  1:3715185/1‑100 (MQ=255)
       ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                             >  2:2545472/1‑100 (MQ=255)
       ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                             <  2:622631/100‑1 (MQ=255)
          gtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATAGCTGTGGTCTGCTGCTGTTAACCACCGATGGTGaa                                                                                          >  2:218264/1‑100 (MQ=255)
           tttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAAc                                                                                         <  2:784566/100‑1 (MQ=255)
           tttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAAc                                                                                         <  2:239213/100‑1 (MQ=255)
           tttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAAc                                                                                         <  2:239639/100‑1 (MQ=255)
             tGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg                                                                                       >  2:2883357/1‑100 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  1:3362884/100‑4 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                      <  2:2583602/100‑4 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                      <  1:1781246/100‑1 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCAGTTCACCATCGGTGGTGAAcagc                                                                                      <  1:2440184/100‑4 (MQ=255)
               aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                     <  1:3038445/100‑1 (MQ=255)
               aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                     <  1:3110240/100‑1 (MQ=255)
               aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCAGTTCACCATCGGTGGTGAAcagca                                                                                     <  1:3000842/100‑5 (MQ=38)
                 cGTCTGCCAAACCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagc                                                                                   <  1:3044284/100‑7 (MQ=38)
                 cGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagc                                                                                   <  1:2870814/100‑7 (MQ=39)
                  gTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaaa                                                                                  >  1:154604/1‑100 (MQ=255)
                  gTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaa                                                                                   >  2:3033442/1‑99 (MQ=255)
                   tCTGCCAAAGCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAc                                                                                 <  2:3106977/100‑1 (MQ=255)
                   tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagcag                                                                                 <  1:876243/100‑9 (MQ=38)
                    cTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAcc                                                                                <  1:2051306/100‑1 (MQ=255)
                       ccAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGtt                                                                             <  1:943777/100‑1 (MQ=255)
                        cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                            >  2:909974/1‑100 (MQ=255)
                         aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTa                                                                           >  1:2606750/1‑100 (MQ=255)
                          aaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTAGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTaa                                                                          <  1:2204581/100‑1 (MQ=255)
                           aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTa                                                                           <  2:3198065/98‑1 (MQ=255)
                            aCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                        <  1:619834/100‑1 (MQ=255)
                               gCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAc                                                                     <  2:628531/100‑1 (MQ=255)
                                cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                    >  2:607575/1‑100 (MQ=255)
                                cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAc                                                                     >  2:1889257/1‑99 (MQ=255)
                                cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTACACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                    <  2:1513130/100‑1 (MQ=255)
                                 gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAccc                                                                   >  1:1353627/1‑100 (MQ=255)
                                 atctcgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAccc                                                                   <  1:1723259/96‑1 (MQ=255)
                                  tGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCa                                                                  >  2:1262843/1‑100 (MQ=255)
                                   ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCaa                                                                 >  1:128802/1‑100 (MQ=255)
                                     cgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCagat                                                               >  2:1723249/1‑97 (MQ=255)
                                       cACGCTGGATTGCCGTGGGTCGCTTGGCCGTTTATACCTGTCGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAGACCGTTTAATGCACCCCAGCCg                                                             >  2:920378/1‑100 (MQ=255)
                                         cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                           <  1:571696/100‑1 (MQ=255)
                                         cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                           <  1:1979337/100‑1 (MQ=255)
                                          gCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                          <  2:950919/100‑1 (MQ=255)
                                             ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                       >  1:3637731/1‑100 (MQ=255)
                                             ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                       >  1:3422548/1‑100 (MQ=255)
                                                 tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaa                                                   <  2:154604/100‑1 (MQ=255)
                                                 tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaa                                                   <  1:2095317/100‑1 (MQ=255)
                                                  gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                  >  2:351421/1‑100 (MQ=255)
                                                   ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                  <  2:3092972/99‑1 (MQ=255)
                                                    cGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGAGAAGTTCAACGt                                                >  2:3604196/1‑100 (MQ=255)
                                                    cGTGGGTCGCTTGGACGTTAATACATGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGt                                                >  2:2062007/1‑100 (MQ=255)
                                                         gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                           >  1:2895997/1‑100 (MQ=255)
                                                         gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                           >  2:399223/1‑100 (MQ=255)
                                                           cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATg                                         <  2:2895976/100‑1 (MQ=255)
                                                             cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGcc                                       >  2:2064198/1‑100 (MQ=255)
                                                              ttGGACGTTAATACCTGTGGTCTGCTGCTGTTCTCCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcg                                      <  1:1445509/100‑1 (MQ=255)
                                                                ggACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                    >  1:1793012/1‑100 (MQ=255)
                                                                 gACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgc                                   >  2:3379122/1‑100 (MQ=255)
                                                                   cGTTAATACCTGTGGGCTGCTGCTGTTAACAACCGATGGTGAACTGGAAAACCGTTTAATGCACCCAAACCGTGAAGTTGAACGTGAATATGCcgtgcgt                                 >  2:2389503/1‑100 (MQ=255)
                                                                          aCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                          <  1:2270576/100‑1 (MQ=255)
                                                                          aCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                          >  1:672780/1‑100 (MQ=255)
                                                                          aCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                          >  1:1329989/1‑100 (MQ=255)
                                                                           ccTGTGGTCTGCTGCTGTTCACCTCCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                         >  2:254525/1‑100 (MQ=255)
                                                                           ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGCACTGTCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGAGAATATGCCGTGCGCGGAtt                            >  2:2362600/1‑97 (MQ=255)
                                                                           ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                         <  2:2148709/100‑1 (MQ=255)
                                                                           ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                         >  1:3204153/1‑100 (MQ=255)
                                                                              gtgGTCTGCTGTTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCa                      <  2:1395989/100‑1 (MQ=255)
                                                                              gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGGGAATATGCCGTGCGTGTATTTGGTCa                      >  2:522100/1‑100 (MQ=255)
                                                                               tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAg                     >  1:2849118/1‑100 (MQ=255)
                                                                                ggTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                    <  1:2879897/100‑1 (MQ=255)
                                                                                 gTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGt                   <  2:1087715/100‑1 (MQ=255)
                                                                                 cggtttgccagTTCACCATCGGTGGTGAACAGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGt                   >  1:1191791/11‑100 (MQ=37)
                                                                                   ctgctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                 <  1:3401874/100‑1 (MQ=255)
                                                                                    tgctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTga                >  2:1904858/1‑100 (MQ=255)
                                                                                     gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac               <  2:1183268/100‑1 (MQ=255)
                                                                                     gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac               >  2:2854011/1‑100 (MQ=255)
                                                                                     gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac               <  1:3347862/100‑1 (MQ=255)
                                                                                      ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg              >  1:3005313/1‑100 (MQ=255)
                                                                                      ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg              >  2:3781750/1‑100 (MQ=255)
                                                                                       gccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacga             >  2:3723581/5‑100 (MQ=255)
                                                                                       gccagTTCACCATCGGTGGTGAACAGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacga             >  1:470130/5‑100 (MQ=38)
                                                                                        ccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac            >  2:3345045/4‑100 (MQ=255)
                                                                                         cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGaagacg           >  2:331789/3‑100 (MQ=255)
                                                                                          tgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc          >  1:2351457/1‑100 (MQ=255)
                                                                                          agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc          >  1:1271358/2‑100 (MQ=255)
                                                                                           gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg         <  2:1793002/100‑1 (MQ=255)
                                                                                              caccacCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaaa      >  2:1667403/1‑100 (MQ=255)
                                                                                               accaccGATGGTGAACTGGCAAACCGTTTACTGCACCCAAGCCGTGAAGTTGACCGTGAATATGCCGTGCGTGTACTTGGTCAGGTTGCCGACNCGaaa      >  1:3155277/1‑99 (MQ=255)
                                                                                                ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGCAACt    >  2:1516496/1‑100 (MQ=255)
                                                                                                ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACt    >  2:3816126/1‑100 (MQ=255)
                                                                                                  accGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc  >  2:731882/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGC  >  NC_000913/1327114‑1327311

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