Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,152,110 A→G 5.2% T251A (ACC→GCC)  rluC → ribosomal large subunit pseudouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,152,1100AG 5.2% 231.8 / 2.2 76T251A (ACC→GCC) rluCribosomal large subunit pseudouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (34/38);  new base G (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

TGCGTTCCGGCGAGCGCATGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTC  >  NC_002947/2151964‑2152259
                                                                                                                                                  |                                                                                                                                                     
tGCGTTCCGGCGAGCGCATGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACAccc                                                                                                                                                    <  2:1195764/150‑1 (MQ=255)
        ggCGAGCGCATGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGcc                                                                                                                                            <  2:283409/150‑1 (MQ=255)
         atctgcgcATGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCg                                                                                                                                           <  1:1168964/146‑1 (MQ=255)
              cgcATGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATACGAGTGAGCACAATGCCTGCCgggggt                                                                                                                                      >  2:1168964/1‑147 (MQ=255)
                  tGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGa                                                                                                                                  <  1:817925/150‑1 (MQ=255)
                  tGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGa                                                                                                                                  <  2:1519251/150‑1 (MQ=255)
                  tGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGa                                                                                                                                  <  1:150806/150‑1 (MQ=255)
                     tCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCg                                                                                                                               >  1:518492/1‑150 (MQ=255)
                     tCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTACGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCAACGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCg                                                                                                                               >  1:1567526/1‑150 (MQ=255)
                       gAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCt                                                                                                                             >  2:1386079/1‑150 (MQ=255)
                         aGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTgg                                                                                                                           <  1:391846/150‑1 (MQ=255)
                            gAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCGTCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGa                                                                                                                        <  2:1442516/150‑1 (MQ=255)
                             aaCGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGAc                                                                                                                       >  1:279517/1‑150 (MQ=255)
                              aCGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGCTGACTTTGGCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACa                                                                                                                      <  1:412777/150‑1 (MQ=255)
                               cGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAg                                                                                                                     <  2:518492/150‑1 (MQ=255)
                                         aaGGAGGCGCTGACCCTGTTTCGCTTCCTGCGCCGTTTTGGTGACTTTGCCACCATGGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGc                                                                                                           >  2:229998/1‑150 (MQ=255)
                                         aaGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACGCCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGc                                                                                                           >  1:838497/1‑150 (MQ=255)
                                         aaGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGc                                                                                                           <  1:74676/150‑1 (MQ=255)
                                                     aCCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCg                                                                                                                                           <  1:414386/106‑1 (MQ=255)
                                                     aCCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCg                                                                                                                                           >  2:414386/1‑106 (MQ=255)
                                                      cccTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCa                                                                                              >  1:513706/1‑150 (MQ=255)
                                                       ccTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAg                                                                                             >  1:1063485/1‑150 (MQ=255)
                                                          gTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTccgccg                                                                                          >  1:980980/1‑150 (MQ=255)
                                                           tttCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCgg                                                                                         <  2:980980/150‑1 (MQ=255)
                                                               gcgTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCAGATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGaaa                                                                                     >  1:1440291/1‑150 (MQ=255)
                                                                    cTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTcgc                                                                                <  1:1386079/150‑1 (MQ=255)
                                                                     tGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTcgcg                                                                               >  2:1602645/1‑150 (MQ=255)
                                                                      gcgcCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGa                                                                              >  1:329199/1‑150 (MQ=255)
                                                                      gcgcCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGa                                                                              <  2:279517/150‑1 (MQ=255)
                                                                          cGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCa                                                                                              >  1:450404/1‑130 (MQ=255)
                                                                          cGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCa                                                                                              <  2:450404/130‑1 (MQ=255)
                                                                            ttttGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGc                                                                                                                    >  2:944944/1‑106 (MQ=255)
                                                                            ttttGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGc                                                                                                                    <  1:944944/106‑1 (MQ=255)
                                                                             ggtggcaaCTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGggc                                                                       >  2:1241191/8‑150 (MQ=255)
                                                                              ttGGTGACTTTGCCACCATCGTCGACGAGAGCTCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGggcg                                                                      <  2:1440291/150‑1 (MQ=255)
                                                                                 gTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCa                                                                   >  2:1364672/1‑150 (MQ=255)
                                                                                          ccaccaTCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCt                                                                                                                             >  1:1061774/1‑83 (MQ=255)
                                                                                          ccaccaTCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCt                                                                                                                             <  2:1061774/83‑1 (MQ=255)
                                                                                          ccaccaTCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGt                                                          >  1:260984/1‑150 (MQ=255)
                                                                                           caccaTCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGtt                                                         >  2:884563/1‑150 (MQ=255)
                                                                                                 cgtcgAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCa                                                   <  2:566768/150‑1 (MQ=255)
                                                                                                   tcgAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCgg                                                                                                                                          <  1:1454915/61‑1 (MQ=255)
                                                                                                   tcgAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCgg                                                                                                                                          >  2:1454915/1‑61 (MQ=255)
                                                                                                      aGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGc                                                                                                                          <  1:668160/74‑1 (MQ=255)
                                                                                                      aGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGc                                                                                                                          >  2:668160/1‑74 (MQ=255)
                                                                                                          gcgcCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTAcgc                                          <  1:1241191/150‑1 (MQ=255)
                                                                                                          gcgcCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTAcgc                                          <  1:1602645/150‑1 (MQ=255)
                                                                                                             cccGATCACCGGTCGTACCCATCAGATCCGCGTGCACGCCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGtt                                       <  2:838497/150‑1 (MQ=255)
                                                                                                              ccGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTg                                      <  1:229998/150‑1 (MQ=255)
                                                                                                              ccGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTg                                      >  1:137479/1‑150 (MQ=255)
                                                                                                              ccGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTg                                      >  1:992040/1‑150 (MQ=255)
                                                                                                               cGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGa                                     <  1:1364672/150‑1 (MQ=255)
                                                                                                               cGAACACCGGTCGTACCGAGCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGa                                     <  2:260984/150‑1 (MQ=255)
                                                                                                                gATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGAc                                    >  1:1553697/1‑150 (MQ=255)
                                                                                                                 aTCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTg                                                                          <  2:873097/111‑1 (MQ=255)
                                                                                                                 aTCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTg                                                                          >  1:873097/1‑111 (MQ=255)
                                                                                                                  tCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAAt                                                                                    <  1:397207/100‑1 (MQ=255)
                                                                                                                  tCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAAt                                                                                    >  2:397207/1‑100 (MQ=255)
                                                                                                                     ccGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACtt                                                                                                <  1:112948/85‑1 (MQ=255)
                                                                                                                     ccGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACtt                                                                                                >  2:112948/1‑85 (MQ=255)
                                                                                                                            tACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGcc                        <  2:574751/150‑1 (MQ=255)
                                                                                                                             aCCCATCAGATCCGCGTGCACGCCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCg                                                                                          >  2:505670/1‑83 (MQ=255)
                                                                                                                             aCCCATCAGATCCGCGTGCACGCCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCg                                                                                          <  1:505670/83‑1 (MQ=255)
                                                                                                                                cATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGcgacga                                                                                                      >  2:778287/1‑68 (MQ=255)
                                                                                                                                cATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGcgacga                                                                                                      <  1:778287/68‑1 (MQ=255)
                                                                                                                                 aTCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGtt                                       >  2:634697/1‑130 (MQ=255)
                                                                                                                                 aTCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGtt                                       <  1:634697/130‑1 (MQ=255)
                                                                                                                                 aTCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGa                                     >  2:490154/1‑132 (MQ=255)
                                                                                                                                 aTCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGa                                     <  1:490154/132‑1 (MQ=255)
                                                                                                                                  tCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGAcgg                  >  2:780445/1‑150 (MQ=255)
                                                                                                                                      aTCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAAt                                                                                    >  2:821316/1‑80 (MQ=255)
                                                                                                                                      aTCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAAt                                                                                    <  1:821316/80‑1 (MQ=255)
                                                                                                                                        ccGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGa            <  2:329199/150‑1 (MQ=255)
                                                                                                                                          gcgTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAAt          >  1:1061584/1‑150 (MQ=255)
                                                                                                                                              gCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGaa      <  1:884563/150‑1 (MQ=255)
                                                                                                                                              gCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGaa      <  2:1061584/150‑1 (MQ=255)
                                                                                                                                                  acCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGa                                                                              >  2:25527/1‑74 (MQ=255)
                                                                                                                                                  acCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGa                                                                              <  1:25527/74‑1 (MQ=255)
                                                                                                                                                  acCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTc  >  1:624699/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCGTTCCGGCGAGCGCATGGTCGAAGTGAACGAGGAAGGCAAGGAGGCGCTGACCCTGTTTCGCGTCCTGCGCCGTTTTGGTGACTTTGCCACCATCGTCGAGGCGCGCCCGATCACCGGTCGTACCCATCAGATCCGCGTGCACACCCTGCATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTC  >  NC_002947/2151964‑2152259

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