Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 970,240 T→G 7.6% F29C (TTC→TGC) ‡ tgt → queuine tRNA‑ribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947970,2400TG 7.6% 258.8 / 9.3 92F29C (TTC→TGC) ‡tgtqueuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (43/42);  new base G (2/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCACGATT  >  NC_002947/970097‑970389
                                                                                                                                               |                                                                                                                                                     
tcaCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGcc                                                                                                                                                 >  1:49637/1‑150 (MQ=255)
tcaCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGcc                                                                                                                                                 >  2:259121/1‑150 (MQ=255)
 caCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCg                                                                                                                                                <  2:1369774/150‑1 (MQ=255)
  aCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCgg                                                                                                                                               >  1:1339864/1‑150 (MQ=255)
  aCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACACCGGCATGAATGCCgg                                                                                                                                               <  1:1254838/150‑1 (MQ=255)
   cccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGt                                                                                                                                              >  2:1914607/1‑150 (MQ=255)
   cccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCggg                                                                                                                                              <  1:1136551/150‑2 (MQ=255)
   cccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCgcg                                                                                                                                              <  2:670102/150‑3 (MQ=255)
    ccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                             >  1:82146/1‑150 (MQ=255)
    ccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                             >  2:65012/1‑150 (MQ=255)
     cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCCCGTGGCACTGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                            <  2:676967/150‑1 (MQ=255)
     cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                            >  1:608001/1‑150 (MQ=255)
      gCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGTAAACCCCGGCATTCATGCCGGTggg                                                                                                                                           <  1:299081/150‑1 (MQ=255)
        aaTCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGa                                                                                                                                         <  2:304355/150‑1 (MQ=255)
         aTCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGAc                                                                                                                                        <  1:1814765/150‑1 (MQ=255)
           ccGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCt                                                                                                                                      <  1:809189/150‑1 (MQ=255)
              gcgcCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATg                                                                                                                                   <  1:1082380/150‑1 (MQ=255)
               cgcCACGCGGGCCCGAGGATCAAGCATGGGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATgg                                                                                                                                  <  2:1603477/150‑1 (MQ=255)
                   aCGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACg                                                                                                                              <  1:587487/150‑1 (MQ=255)
                             gAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCa                                                                                                                    >  1:473927/1‑150 (MQ=255)
                             gAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGACAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                             <  1:767721/125‑1 (MQ=255)
                             gAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGACAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                             >  2:767721/1‑125 (MQ=255)
                                     aGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCa                                                                                                            >  1:1408619/1‑150 (MQ=255)
                                      gCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAc                                                                                                           <  2:82146/150‑1 (MQ=255)
                                      gCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAc                                                                                                           >  2:1825508/1‑150 (MQ=255)
                                         tGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgcg                                                                                                        >  2:1035220/1‑150 (MQ=255)
                                         tGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgcg                                                                                                        <  1:1825508/150‑1 (MQ=255)
                                          gAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGATACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGa                                                                                                       <  2:364359/150‑1 (MQ=255)
                                             tCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACAt                                                                                                    >  1:1266862/1‑150 (MQ=255)
                                               gCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGc                                                                                                                     <  2:1429635/131‑1 (MQ=255)
                                               gCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGc                                                                                                                     >  1:1429635/1‑131 (MQ=255)
                                                cACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGa                                                                                                 <  2:879657/150‑1 (MQ=255)
                                                   cTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGc                                                                                              >  2:852132/1‑150 (MQ=255)
                                                       cGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGACACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATc                                                                                          <  2:1246898/150‑1 (MQ=255)
                                                         aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGTTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAg                                                                                                <  1:1087296/142‑1 (MQ=255)
                                                         aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCAc                                                                                                                               >  1:1131475/1‑111 (MQ=255)
                                                         aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCAc                                                                                                                               <  2:1131475/111‑1 (MQ=255)
                                                         aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAg                                                                                                >  2:1087296/1‑142 (MQ=255)
                                                         aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCgg                                                                                        >  1:1896210/1‑150 (MQ=255)
                                                         aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATAGAGGCCCTCgg                                                                                        >  1:999817/1‑150 (MQ=255)
                                                            gTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGgcgc                                                                                     >  1:581009/1‑150 (MQ=255)
                                                              ccTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCCGCCACGCGACATCGAGGCCATCGGCGCCg                                                                                   <  2:234610/150‑1 (MQ=255)
                                                               cTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCga                                                                                  <  1:14540/150‑1 (MQ=255)
                                                                ttCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCggg                                                                                                                                              >  2:1047252/1‑88 (MQ=255)
                                                                ttCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCggg                                                                                                                                              <  1:1047252/89‑2 (MQ=255)
                                                                ttCGAACTGCTGGCCACCGACGGCAAGGCCCGCCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCgag                                                                                 >  1:320365/1‑150 (MQ=255)
                                                                 tCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCgcca                                                                                <  2:999817/150‑4 (MQ=255)
                                                                     aCTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCAACGTTTCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAgatgat                                                                            >  1:839183/1‑150 (MQ=255)
                                                                       tgctgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATcc                                                                          >  2:1931422/1‑150 (MQ=255)
                                                                       tgctgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACCTAGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATcc                                                                          >  1:74522/1‑150 (MQ=255)
                                                                             ccACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAg                                                                                                >  1:1211052/1‑122 (MQ=255)
                                                                             ccACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAg                                                                                                <  2:1211052/122‑1 (MQ=255)
                                                                              cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACg                                                                                                                              >  1:916440/1‑91 (MQ=255)
                                                                              cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACg                                                                                                                              <  2:916440/91‑1 (MQ=255)
                                                                               aCCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAac                                                                  <  1:1914607/150‑1 (MQ=255)
                                                                                ccGACGGCAAGGCCCGTACGGGTCGCCTAAACTTGTTTCGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAaca                                                                 >  2:992910/1‑150 (MQ=255)
                                                                                 cGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGCACTTCGGCGCCGAGATGATCCTCGGCAacac                                                                <  1:1022170/150‑1 (MQ=255)
                                                                                 cGACGGCAAGGCCCGTCCAGTGCGTATCACCGTCCCACGTGGCACCGTGGAAACCAAGGCATTCATGCAGGTGGGGACCTATAGCACGCTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAacac                                                                <  1:1069600/150‑1 (MQ=255)
                                                                                  gACGGCAAGGCCCGCCGCGGTCGCATCACCTTCCCAGGTGGCAGCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACAcc                                                               >  2:514430/1‑150 (MQ=255)
                                                                                     ggCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTc                                                            <  1:276961/150‑1 (MQ=255)
                                                                                       cAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCa                                                          >  1:313660/1‑150 (MQ=255)
                                                                                       cAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCa                                                          >  1:1241474/1‑150 (MQ=255)
                                                                                            cccGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgt                                                     <  2:313660/150‑1 (MQ=255)
                                                                                             ccGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgtg                                                    <  1:1931422/150‑1 (MQ=255)
                                                                                              cGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACGTCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgtg                                                    >  1:1156156/1‑149 (MQ=255)
                                                                                              cGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACGTCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgtg                                                    <  2:1156156/149‑1 (MQ=255)
                                                                                               gTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAac                                                                  <  1:1032035/134‑1 (MQ=255)
                                                                                               gTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAac                                                                  >  2:1032035/1‑134 (MQ=255)
                                                                                                tCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCCTCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCt                                                 <  2:1408619/150‑1 (MQ=255)
                                                                                                    ggTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTc                                                                        >  1:840621/1‑123 (MQ=255)
                                                                                                    ggTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTc                                                                        <  2:840621/123‑1 (MQ=255)
                                                                                                    ggTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGt                                             <  1:1035220/150‑1 (MQ=255)
                                                                                                    ggTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGt                                             >  1:1668985/1‑150 (MQ=255)
                                                                                                         cATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCggg                                        >  2:1732830/1‑150 (MQ=255)
                                                                                                         cATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCggg                                        >  1:1519187/1‑150 (MQ=255)
                                                                                                          aTCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGAc                                                                                                                                        <  2:574949/53‑1 (MQ=255)
                                                                                                          aTCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGAc                                                                                                                                        >  1:574949/1‑53 (MQ=255)
                                                                                                               cTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGa                                  >  1:1543386/1‑150 (MQ=255)
                                                                                                                  cccACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGt                               <  1:1467028/150‑1 (MQ=255)
                                                                                                                   ccACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTg                              >  1:1720859/1‑150 (MQ=255)
                                                                                                                    cACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGa                             <  2:1831843/150‑1 (MQ=255)
                                                                                                                      cGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATc                           <  2:1796692/150‑1 (MQ=255)
                                                                                                                      cGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATc                           <  2:581009/150‑1 (MQ=255)
                                                                                                                      cGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATc                           <  2:1339864/150‑1 (MQ=255)
                                                                                                                           cACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCaagaa                      <  1:992910/150‑1 (MQ=255)
                                                                                                                            aCCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCaagaag                     <  2:1720859/150‑1 (MQ=255)
                                                                                                                            aCCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCaagaag                     <  2:839183/150‑1 (MQ=255)
                                                                                                                                    cccaccGGCATGAATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACgg             >  2:370526/5‑150 (MQ=255)
                                                                                                                                     aCCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGc            >  2:187882/1‑150 (MQ=255)
                                                                                                                                         cGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGc        >  1:287004/1‑150 (MQ=255)
                                                                                                                                         cGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGc        >  1:1430351/1‑150 (MQ=255)
                                                                                                                                             aTTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCACGa    >  1:540198/1‑150 (MQ=255)
                                                                                                                                              ttCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCACGAt   >  1:1804/1‑150 (MQ=255)
                                                                                                                                               tCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCCAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGAACCGATGTGATCAAGAAGCACAACGGCCTGCACGAtt  <  1:155753/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCACGATT  >  NC_002947/970097‑970389

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