Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,216 A→G 16.4% D122G (GAT→GGT)  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,2160AG16.4% 283.6 / 48.3 128D122G (GAT→GGT) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (62/45);  new base G (11/10);  total (73/55)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 129 total aligned reads displayed.

CAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTG  >  NC_000913/1327117‑1327314
                                                                                                   |                                                                                                  
cAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                    >  2:1347724/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)
    ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                                >  1:3715185/1‑100 (MQ=255)
    ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                                <  2:622631/100‑1 (MQ=255)
    ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                                >  2:2545472/1‑100 (MQ=255)
       gtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATAGCTGTGGTCTGCTGCTGTTAACCACCGATGGTGaa                                                                                             >  2:218264/1‑100 (MQ=255)
        tttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAAc                                                                                            <  2:239213/100‑1 (MQ=255)
        tttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAAc                                                                                            <  2:784566/100‑1 (MQ=255)
        tttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAGGGTGAAc                                                                                            <  2:1998028/100‑1 (MQ=255)
          tGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg                                                                                          >  2:2883357/1‑100 (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)
            aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                        <  1:1347733/100‑1 (MQ=255)
            aCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCAGTTCACCATCGGTGGTGAAcagca                                                                                        <  1:3000842/100‑5 (MQ=38)
             ccGTCTGCCAATACTGCGTGGCGCACGCTGGGTTTCCGTGGGTCGGTTGGACGTTAAGACCTGTGGTCTGCTGCTGTTCACCAGCGGTGGTGAAcagc                                                                                         <  2:2479214/98‑4 (MQ=255)
             ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcag                                                                                       <  2:2098125/100‑6 (MQ=255)
             ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCa                                                                                       >  1:2133929/1‑100 (MQ=255)
              cGTCTGCCAAACCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagc                                                                                      <  1:3044284/100‑7 (MQ=38)
              cGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagc                                                                                      <  1:2870814/100‑7 (MQ=39)
               gTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCaa                                                                                      >  2:3033442/1‑99 (MQ=255)
                tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagcag                                                                                    <  1:876243/100‑9 (MQ=38)
                 cTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAcc                                                                                   <  1:2051306/100‑1 (MQ=255)
                 cTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGAAGTTAATACCTGTGGTCTGCTGCTGTTCACCACAGATGGTGAACTAGCAAAc                                                                                    >  1:2961501/1‑99 (MQ=255)
                    ccAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGtt                                                                                <  1:943777/100‑1 (MQ=255)
                     cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                               >  2:909974/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)
                         aCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGAGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                           >  1:2336604/1‑100 (MQ=255)
                             cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAcc                                                                       >  2:607575/1‑100 (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:3383429/1‑100 (MQ=255)
                                ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCaa                                                                    >  1:128802/1‑100 (MQ=255)
                                ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCaa                                                                    >  1:1849521/1‑100 (MQ=255)
                                  cgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCagat                                                                  >  2:1723249/1‑97 (MQ=255)
                                    cACGCTGGATTGCCGTGGGTCGCTTGGCCGTTTATACCTGTCGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAGACCGTTTAATGCACCCCAGCCg                                                                >  2:920378/1‑100 (MQ=255)
                                      cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                              >  1:1613738/1‑100 (MQ=255)
                                      cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                              <  1:1979337/100‑1 (MQ=255)
                                       gCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                             <  2:950919/100‑1 (MQ=255)
                                        cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCATGCCGTGaa                                                            >  2:4547/1‑100 (MQ=255)
                                        cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGaa                                                            >  2:3480265/1‑100 (MQ=255)
                                        cTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                             <  2:674801/99‑1 (MQ=255)
                                          ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                          >  1:3637731/1‑100 (MQ=255)
                                          ggATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                          >  1:3422548/1‑100 (MQ=255)
                                              tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGaa                                                      <  2:154604/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)
                                                     ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAat                                               >  1:2760584/1‑100 (MQ=255)
                                                          cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGcc                                          >  2:2064198/1‑100 (MQ=255)
                                                           ttGGACGTTAATACCTGTGGTCTGCTGCTGTTCTCCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcg                                         <  1:1445509/100‑1 (MQ=255)
                                                                cGTTAATACCTGTGGGCTGCTGCTGTTAACAACCGATGGTGAACTGGAAAACCGTTTAATGCACCCAAACCGTGAAGTTGAACGTGAATATGCcgtgcgt                                    >  2:2389503/1‑100 (MQ=255)
                                                                   tAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTa                                 >  1:529460/1‑100 (MQ=255)
                                                                       aCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                             >  1:672780/1‑100 (MQ=255)
                                                                        ccTGTGGTCTGCTGCTGTTCACCTCCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                            >  2:254525/1‑100 (MQ=255)
                                                                        ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                            >  1:3204153/1‑100 (MQ=255)
                                                                        ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                            <  2:2148709/100‑1 (MQ=255)
                                                                        ccTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                             >  2:1535073/1‑99 (MQ=255)
                                                                          tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTc                          <  1:3033466/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:1263947/1‑100 (MQ=255)
                                                                                  gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac                  >  2:2854011/1‑100 (MQ=255)
                                                                                  gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac                  <  1:3347862/100‑1 (MQ=255)
                                                                                  gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac                  <  2:1183268/100‑1 (MQ=255)
                                                                                   tgccagTTCACCATCGGTGGTGAACAGGCAAACCGTTTAATGCACCCCAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgccg                 >  2:2378592/6‑100 (MQ=37)
                                                                                   ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg                 >  2:3781750/1‑100 (MQ=255)
                                                                                    gccagTTCACCATCGGTGGTGAACAGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacga                >  1:470130/5‑100 (MQ=38)
                                                                                     gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac               <  2:2665740/100‑1 (MQ=255)
                                                                                     ccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac               >  2:3345045/4‑100 (MQ=255)
                                                                                      cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGaagacg              >  2:331789/3‑100 (MQ=255)
                                                                                      cagTTCACCATCGGTGGTGAACTGGCAAACCGCTTAATGCACCCAACCCGTGAAGTTGAACGAGAATATGCCGTGCGTGTATTTGGTCAGGTAgacg                 >  2:1728604/3‑97 (MQ=255)
                                                                                       tgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  1:2351457/1‑100 (MQ=255)
                                                                                       tgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  2:3767276/1‑100 (MQ=255)
                                                                                       agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  2:2945173/2‑100 (MQ=255)
                                                                                       agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  1:1271358/2‑100 (MQ=255)
                                                                                       agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc             >  2:1155266/2‑100 (MQ=255)
                                                                                        gTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg            <  2:1793002/100‑1 (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)
                                                                                                 cGATGGTGAACTGGCAAACGGTTTAATGCACACAAGACGTGAAGTTGCACGTGAATATGCCGTGCGCGTATTTGGTCAGGTTGTCGACGCGAAGCTg      >  2:2970819/1‑97 (MQ=255)
                                                                                                 cGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt   <  2:3383403/100‑1 (MQ=255)
                                                                                                 cGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt   <  1:2206909/100‑1 (MQ=255)
                                                                                                 cGATGGTGAACGGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt   <  1:1251203/100‑1 (MQ=255)
                                                                                                  gATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTg  <  1:1716238/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTG  >  NC_000913/1327117‑1327314

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: