Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,213 C→T 15.3% 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,2130CT15.3% 326.6 / 40.3 132T121I (ACC→ATC) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base C (57/53);  new base T (13/7);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 142 total aligned reads displayed.

GTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCG  >  NC_000913/1327114‑1327312
                                                                                                   |                                                                                                   
gTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTcaccac                                                                                                     <  2:1610358/100‑1 (MQ=255)
 tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                    >  1:2112010/1‑100 (MQ=255)
 tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                    <  1:3642938/100‑1 (MQ=255)
     aCGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATg                                                                                                <  1:3199314/100‑1 (MQ=255)
     aCGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATg                                                                                                <  2:2401608/100‑1 (MQ=255)
       ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                              <  2:2328476/100‑1 (MQ=255)
        gtgtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGgtggtg                                                                                             <  1:2492415/100‑1 (MQ=255)
        gtgtTTGACCGTCGGCCAAAACTGCGTGGCGCACGCTTGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTg                                                                                             <  2:1746067/100‑1 (MQ=255)
        gtgtTTGAACGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCCCCACCGATGGTg                                                                                             <  2:1331132/100‑1 (MQ=255)
          gtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGaa                                                                                           >  1:1389537/1‑100 (MQ=255)
            ttGACCGTCTGCCGAAACTGCGTGGCGCACGCTGGATTGCCGGGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACt                                                                                         >  1:2673125/1‑100 (MQ=255)
             tGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcag                                                                                        <  2:122327/100‑3 (MQ=255)
             tGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg                                                                                        <  2:2199447/100‑1 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                       <  1:1524377/100‑4 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                       <  1:937292/100‑1 (MQ=255)
              gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                       >  2:3535433/1‑100 (MQ=255)
               aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagca                                                                                      <  2:2057164/100‑5 (MQ=255)
                ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcag                                                                                     <  1:2897818/100‑6 (MQ=255)
                ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCAGCGAGGGTGAACTGGc                                                                                      >  1:178360/1‑99 (MQ=255)
                  gTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaaa                                                                                   >  2:652636/1‑100 (MQ=255)
                    cTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAcc                                                                                 <  1:3425237/100‑1 (MQ=255)
                     tGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGCACTGGCAAAc                                                                                  >  2:1923845/1‑98 (MQ=255)
                     tGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCg                                                                                >  1:3435926/1‑100 (MQ=255)
                      gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCg                                                                                >  2:3196497/1‑99 (MQ=255)
                      gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                               >  1:1023670/1‑100 (MQ=255)
                      gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                               <  2:3310391/100‑1 (MQ=255)
                      gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                               >  1:3009629/1‑100 (MQ=255)
                       ccAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGtt                                                                              >  1:66183/1‑100 (MQ=255)
                        cAAAACTGCGTGGGGCACGCTGGATTGCCGTGGGTCGCGTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                             >  2:2434028/1‑100 (MQ=255)
                        cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                             <  1:1028257/100‑1 (MQ=255)
                        cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                             <  2:2612535/100‑1 (MQ=255)
                         aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                             <  2:3412387/99‑1 (MQ=255)
                          aaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTaa                                                                           >  2:2528539/1‑100 (MQ=255)
                           aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAAt                                                                          <  1:3097317/100‑1 (MQ=255)
                            aCTGCGTGGCGCACGCTGGATTGCCGTGTGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTaacg                                                                         <  1:2700877/100‑3 (MQ=255)
                            aCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                         <  1:676683/100‑1 (MQ=255)
                             cTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGc                                                                        >  1:541665/1‑100 (MQ=255)
                              tGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCa                                                                       <  1:1941671/100‑1 (MQ=255)
                              tGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCa                                                                       >  1:1614986/1‑100 (MQ=255)
                               gCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAc                                                                      >  1:3708552/1‑100 (MQ=255)
                                 gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAccc                                                                    <  2:591826/100‑1 (MQ=255)
                                 gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAccc                                                                    >  1:405851/1‑100 (MQ=255)
                                 gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAccc                                                                    <  2:2397500/100‑1 (MQ=255)
                                   ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCaa                                                                  >  1:2152872/1‑100 (MQ=255)
                                     cgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGc                                                                >  1:1065762/1‑100 (MQ=255)
                                         cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                            <  2:2013808/100‑1 (MQ=255)
                                            tGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAg                                                         <  2:1795839/100‑1 (MQ=255)
                                             ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                        >  2:3076995/1‑100 (MQ=255)
                                             ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                        <  2:3708524/100‑1 (MQ=255)
                                             ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTCGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                        <  1:3748021/100‑1 (MQ=255)
                                                  gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaat                                                   <  2:1023664/100‑2 (MQ=255)
                                                  gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                   >  2:3636621/1‑100 (MQ=255)
                                                  gCCGTGGGGCGCTTGGACGTTAATACCTGTGGTCTGCTGCGGTTCCCCACGGATGGTGCGCTGGCACCCCGGTTGATGCCCCCATGGCGGGAAGTTGGAc                                                   >  2:377799/1‑100 (MQ=255)
                                                   ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                  <  1:323669/100‑1 (MQ=255)
                                                     gTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGAGAAGTTGAACGTg                                                >  2:1208365/1‑100 (MQ=255)
                                                       gggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGaa                                              >  2:2894651/1‑100 (MQ=255)
                                                        ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAat                                             >  1:3787964/1‑100 (MQ=255)
                                                         gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                            >  2:2568750/1‑100 (MQ=255)
                                                         gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                            >  1:1746330/1‑100 (MQ=255)
                                                                ggACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                     <  1:3344161/100‑1 (MQ=255)
                                                                ggACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                     <  1:1972315/100‑1 (MQ=255)
                                                                  aCGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                   <  2:88939/100‑1 (MQ=255)
                                                                   cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgt                                      >  2:2267358/1‑96 (MQ=255)
                                                                    gTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcgtg                                 >  1:3649341/1‑100 (MQ=255)
                                                                       aaTACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAt                              >  2:2670298/1‑100 (MQ=255)
                                                                         ttccTGTGGCCTGCTGCTGTTCACCACCGATGGTGAACTGGAAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAttt                            <  1:2173553/98‑1 (MQ=255)
                                                                         tACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAttt                            >  2:1927599/1‑100 (MQ=255)
                                                                         tACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAttt                            <  1:1226425/100‑1 (MQ=255)
                                                                           ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                          <  2:647245/100‑1 (MQ=255)
                                                                           ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                          <  2:1132493/100‑1 (MQ=255)
                                                                            cTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGt                         >  2:445390/1‑100 (MQ=255)
                                                                             tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTc                        <  1:3063536/100‑1 (MQ=255)
                                                                             tgtgGTCTGCTGCTGTACACCAACGATGGTGAACTGGCAAACCGTTTAATGCACCCAACCCTGTAAGTTTAACGTGAATATTCCGTGCGGGTATGTGGTc                        >  2:2297545/1‑100 (MQ=255)
                                                                              gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCa                       <  2:1067616/100‑1 (MQ=255)
                                                                                ggTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                     <  2:870969/100‑1 (MQ=255)
                                                                                ggTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                     <  2:66183/100‑1 (MQ=255)
                                                                                  tCTGCTGCTGTTCACCACCGATGGTGAACGGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGACCGTGAATATGCCGTGCGTGTATTTGGTCAGGtt                   <  2:2369534/100‑1 (MQ=255)
                                                                                   gtttgccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                  >  1:3908227/9‑100 (MQ=38)
                                                                                    tgctgctgTTCACCACCGGTGGTGAACTGGCAAACCGGTTAATGCACCCAAGCCGTGAAGTTGAACGTGGATATGCCGTGCGAGTATTTGGTCAGGTTg                  >  2:1389510/1‑99 (MQ=255)
                                                                                    tgctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTga                 >  2:736731/1‑100 (MQ=255)
                                                                                      tgccagTTCACCATCGGTGGTGAACTGTCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCATGCGTGTATTTGGTCAGGTTgacg               >  2:1636295/6‑100 (MQ=255)
                                                                                       tgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacga              >  1:3037173/1‑100 (MQ=255)
                                                                                       gccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacga              >  1:2717035/5‑100 (MQ=255)
                                                                                        gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac             >  1:1707819/1‑100 (MQ=255)
                                                                                        ccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTTGTCAGGTTgacgac             >  1:884462/4‑100 (MQ=255)
                                                                                         ctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg            >  2:3823653/1‑100 (MQ=255)
                                                                                         cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg            >  1:2423487/3‑100 (MQ=255)
                                                                                         cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg            >  2:3433316/3‑100 (MQ=255)
                                                                                         cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg            >  1:2621566/3‑100 (MQ=255)
                                                                                         cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg            >  1:3862268/3‑100 (MQ=255)
                                                                                          tgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc           <  2:2055597/100‑1 (MQ=255)
                                                                                          agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc           >  2:1538209/2‑100 (MQ=255)
                                                                                          agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc           >  1:1102719/2‑100 (MQ=255)
                                                                                           gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGacgacg            <  1:1586911/98‑1 (MQ=255)
                                                                                           gTTCACCACCGATGGTGAACTGGCAAACCGCTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg          >  2:2824594/1‑100 (MQ=255)
                                                                                             tCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaa        >  1:2245988/1‑100 (MQ=255)
                                                                                                ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACt     >  2:3905801/1‑100 (MQ=255)
                                                                                                ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACt     >  1:751785/1‑100 (MQ=255)
                                                                                                  accGATGGTGAACTGGCAAACCGATTAATGCACCCAAGCAGGGAAGTTGAACGTGTATATGACATGGGTGTAAGTGGTGGGGTTGAGTACGCGAAATTg    >  2:1307355/1‑99 (MQ=255)
                                                                                                   ccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg  >  2:1515329/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCG  >  NC_000913/1327114‑1327312

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