Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,213 C→T 5.7% 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,2130CT 5.7% 286.5 / 4.3 91T121I (ACC→ATC) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base T (2/3);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 107 total aligned reads displayed.

CTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAA  >  NC_000913/1327064‑1327362
                                                                                                                                                     |                                                                                                                                                     
cTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTcaccac                                                                                                                                                       >  1:73596/1‑150 (MQ=255)
 tATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGCTAATACCTGTGGTCTGCTGCTGTTCaccacc                                                                                                                                                      >  1:814126/1‑150 (MQ=255)
     aCAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATg                                                                                                                                                  <  1:1260554/150‑1 (MQ=255)
      cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAATACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACTTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                                                                                 >  1:227359/1‑150 (MQ=255)
      cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGTCGCAAGCTTGATTTCCGTGGGTGGCTTGGACGTTAATACATGTGGTGTGCTGCTGTTCACCAGCGATgg                                                                                                                                                 >  2:1212694/1‑150 (MQ=255)
      cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                                                                                 >  1:1631957/1‑150 (MQ=255)
      cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAAAACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                                                                                 >  1:460912/1‑150 (MQ=255)
        aaCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTg                                                                                                                                               <  1:1309916/150‑1 (MQ=255)
        aaCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGACGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTg                                                                                                                                               >  2:27400/1‑150 (MQ=255)
          ccGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGaa                                                                                                                                             <  1:1232120/150‑1 (MQ=255)
           cGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCTCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAAc                                                                                                                                            <  2:210382/150‑1 (MQ=255)
           cGGAAGGTGAGTTGAGTACCCGTAACGATCCGGACGGACGTCCAACGGTGTTTGACCGTCTGCAAAAACTGCGTGGCGCACGCTGTATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAc                                                                                                                                            <  1:1218316/150‑1 (MQ=255)
            ggAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAca                                                                                                                                           <  2:1016059/150‑2 (MQ=255)
             gAAGGTGAGTGGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTg                                                                                                                                          <  1:874071/150‑1 (MQ=255)
               aGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                                                                        >  2:75183/1‑150 (MQ=255)
                ggTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCTCCGATGGTGAACTGGCa                                                                                                                                       >  1:1028491/1‑150 (MQ=255)
                 gTGAGTTGTGTACCCGTAACGAGCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagc                                                                                                                                      <  2:1845542/150‑7 (MQ=255)
                  tGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCGGTTCACCACCGATGGTGAACTGGCaaa                                                                                                                                     <  1:1051783/150‑1 (MQ=255)
                   gAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGGTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAc                                                                                                                                    <  2:1385482/150‑1 (MQ=255)
                    aGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGTACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAcc                                                                                                                                   <  1:1137579/150‑1 (MQ=255)
                     gTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCg                                                                                                                                  >  1:1273877/1‑150 (MQ=255)
                     gTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCg                                                                                                                                  >  1:642339/1‑150 (MQ=255)
                     gTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCg                                                                                                                                  <  1:1852249/150‑1 (MQ=255)
                       tgtgTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGtt                                                                                                                                >  1:1590292/1‑150 (MQ=255)
                        gtgtACCCGTCACGATCCGGAAGGACGTCCAACGGGGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                                                                               <  2:580910/150‑1 (MQ=255)
                        gtgtACCCGTAACGATCCGGAAGGACGTCCAACGTTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                                                                               <  2:88587/150‑1 (MQ=255)
                        gtgtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                                                                               <  1:1495533/150‑1 (MQ=255)
                         tgtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTa                                                                                                                              >  1:660437/1‑150 (MQ=255)
                                       tCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCg                                                                                                                <  2:415662/150‑1 (MQ=255)
                                       gccGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCg                                                                                                                <  1:1643034/149‑1 (MQ=255)
                                         cGGAAGGACGTCAAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGGACGCTGGATTGCCGTGGGGCGCTTGGACGTTATTACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTg                                                                                                              <  1:443961/150‑1 (MQ=255)
                                          ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACAACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                             >  2:1556645/1‑150 (MQ=255)
                                          ggAAGGACGTCCAACGGTGTTTGACCGGCTGCCAAAACTGCGGGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                             <  1:219839/150‑1 (MQ=255)
                                             aGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                                                                          <  1:1496149/150‑1 (MQ=255)
                                             aGGACGTCCAACGGGGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGGCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                                                                          <  2:1681898/150‑1 (MQ=255)
                                               gACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCCCCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTg                                                                                                        <  2:702096/150‑1 (MQ=255)
                                               gACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTg                                                                                                        >  2:1370346/1‑150 (MQ=255)
                                                aCGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGa                                                                                                       <  1:383265/150‑1 (MQ=255)
                                                  gTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                                                                     <  1:1099420/150‑1 (MQ=255)
                                                  gTCCAACGGAGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGTTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAAc                                                                                                     <  1:263787/150‑1 (MQ=255)
                                                   tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                    >  1:880237/1‑150 (MQ=255)
                                                   tCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACg                                                                                                    >  1:1393465/1‑150 (MQ=255)
                                                     cAACGGTGTTTGCCCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTg                                                                                                  >  2:1489163/1‑150 (MQ=255)
                                                     cAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTg                                                                                                  >  2:886277/1‑150 (MQ=255)
                                                      aaCGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCCCCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGa                                                                                                 >  2:862652/1‑150 (MQ=255)
                                                              ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGGGGGGCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACGGGCAAACCGTTTAATGCACCCAAGCCGGGAAGTTGACGGTGAATTTGCcg                                                                                         >  2:302092/1‑150 (MQ=255)
                                                                gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                       >  1:806607/1‑150 (MQ=255)
                                                                gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGACAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                       <  1:1661743/150‑1 (MQ=255)
                                                                        gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                                                                               >  2:1047750/1‑150 (MQ=255)
                                                                        gCCAAAACTGCGGGGCGCACGCTGGATTGCCGGGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTAAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                                                                               <  2:86790/150‑1 (MQ=255)
                                                                          cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                                                                             >  1:1633262/1‑150 (MQ=255)
                                                                                tGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                                                                       <  2:1628882/150‑1 (MQ=255)
                                                                                tGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAgg                                                                       >  1:1264130/1‑150 (MQ=255)
                                                                                  cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACAGGAAAACCGTTTAATGCACACAAGACTTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGtt                                                                     >  2:1097910/1‑150 (MQ=255)
                                                                                     ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCAGTGAAGTTGAGCGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac                                                                  >  1:1289937/1‑150 (MQ=255)
                                                                                      gctcACGCTGGATTGCCGTGGGTCGCTTGTATGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg                                                                 <  1:303912/147‑1 (MQ=255)
                                                                                      gcgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacg                                                                 >  2:877819/1‑150 (MQ=255)
                                                                                        gcACGCTGGATTGCCGTGGGTCGCTGGGACGTTAATACCTCTGGTCTGCTGCTGTTCACCACCGAGGGTGAACTGGCTAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac                                                               <  1:304939/150‑1 (MQ=255)
                                                                                           cGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAACCCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgcg                                                            >  2:101161/1‑150 (MQ=255)
                                                                                             cTGGATTGCCGTGGGTCGCTTGGAAGTTAATACCTGTGGTCTGCTGCTGTTCTCCACCGATGGTGAACTGGCTAACCGTTTAATGCCCCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTCTATTTGGTCAGCTTGACGACCCGaa                                                          >  2:319980/1‑150 (MQ=255)
                                                                                                  ttGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc                                                     >  2:1875507/1‑150 (MQ=255)
                                                                                                       gTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                <  2:806607/150‑1 (MQ=255)
                                                                                                       gTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                <  1:406324/150‑1 (MQ=255)
                                                                                                         gggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAttt                                              <  2:1590292/150‑1 (MQ=255)
                                                                                                         gggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAAGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAttt                                              <  1:583089/150‑1 (MQ=255)
                                                                                                          agTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                             >  1:991827/2‑150 (MQ=255)
                                                                                                            tCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAg                                           <  1:620532/150‑1 (MQ=255)
                                                                                                            gcGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAg                                           <  1:1662125/149‑1 (MQ=255)
                                                                                                              gCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTc                                         <  1:870246/150‑1 (MQ=255)
                                                                                                               cTTGGACGTTAATCCCTGTGGTCTGCTGCTGATCCCCCCCGATGGTGAAATCGCAAACCGTTTAATGCACCCTAGCAGTGAAGGTCAACGTGAATATGCCGTGCGAGTATTTGGTCAGGTTGCCGACGCGAAACTGCGTGATTTGAGTCg                                        >  2:789660/1‑150 (MQ=255)
                                                                                                                  ggACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTGGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTgg                                     <  1:512596/150‑1 (MQ=255)
                                                                                                                   gACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGc                                    >  1:423856/1‑150 (MQ=255)
                                                                                                                     cGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGGGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGt                                  <  1:217081/150‑1 (MQ=255)
                                                                                                                      gTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                 <  2:1287986/150‑1 (MQ=255)
                                                                                                                       ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGc                                <  1:1737514/150‑1 (MQ=255)
                                                                                                                       ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGc                                <  2:1631957/150‑1 (MQ=255)
                                                                                                                       ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCCTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGc                                >  2:1244967/1‑150 (MQ=255)
                                                                                                                       ttAATACCTGTGGTCTGCGGCTGTTCACCAGCGATGGTGAACTGGCATACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGGGCGTGTATTTGGGCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGc                                <  1:141839/150‑1 (MQ=255)
                                                                                                                           aaCCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGtt                            <  2:73596/149‑1 (MQ=255)
                                                                                                                              cTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGa                         >  2:1783792/1‑150 (MQ=255)
                                                                                                                              cTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGa                         <  1:47248/150‑1 (MQ=255)
                                                                                                                               tgtgGTCTGCTGCTGTTCACAACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATAAGCCGGGCGTGTATTTGGTCAGGTTGACGCCGCGAAACTGGGTGATTTGAGTCGTGGCGTGCAGGTGGaa                        >  2:1244519/1‑150 (MQ=255)
                                                                                                                                  ggTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGAt                     <  2:814126/150‑1 (MQ=255)
                                                                                                                                   ggctgctgctgTGCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATg                    <  2:1426174/148‑1 (MQ=255)
                                                                                                                                    tCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATgg                   >  2:788308/1‑150 (MQ=255)
                                                                                                                                       gctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTTCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTcc                >  2:138813/1‑150 (MQ=255)
                                                                                                                                         tgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAGCGTGACTGTGCCGTGCGTGTATTTGGGCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGGTGGTAGATGGTCCgg              >  1:1209488/1‑150 (MQ=255)
                                                                                                                                           cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCCCCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGTAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCa            >  1:1146446/3‑150 (MQ=255)
                                                                                                                                                  ccaccGATGGTGAACTGGCAAACCGTTTATTGCACCCACGCCGTGAAGTTGAACGTGAATATGCCGGGCGTGTAGTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCttttt     >  2:1512746/1‑149 (MQ=255)
                                                                                                                                                  ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTa     >  1:1281904/1‑150 (MQ=255)
                                                                                                                                                  ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTa     >  1:686924/1‑150 (MQ=255)
                                                                                                                                                  ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTa     >  1:583870/1‑150 (MQ=255)
                                                                                                                                                   caccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATGTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaa    >  1:1757917/1‑150 (MQ=255)
                                                                                                                                                   caccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGCGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaa    <  2:1857831/150‑1 (MQ=255)
                                                                                                                                                   caccGATGGTGAACTGGCAAACCGTTTAAAGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGGCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTTGAAGATGGTCCGGCAGCTTTTaa    <  1:1496268/150‑1 (MQ=255)
                                                                                                                                                    accGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGGCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaa   <  2:1633262/150‑1 (MQ=255)
                                                                                                                                                    accGATGGTGAACTGGCAAAACGTTTAATGCACCCAATCCGTGAATTTGATCGTGAATATGCCGTGCGTGTATATGGTCAGGTTTACTACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTGGGAAGATGGTCCGGCAGCTTTTaaa   <  2:1381291/150‑1 (MQ=255)
                                                                                                                                                    accGAGGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaa   <  1:1068150/150‑1 (MQ=255)
                                                                                                                                                     ccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaaa  <  2:660437/150‑1 (MQ=255)
                                                                                                                                                     ccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaaa  >  2:390154/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTATTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAA  >  NC_000913/1327064‑1327362

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