Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,851,958 C→G H334Q (CAC→CAG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,9580CG90.4% 291.8 / 19.3 94H334Q (CAC→CAG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base G (43/42);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAA  >  NC_002947/1851812‑1852104
                                                                                                                                                  |                                                                                                                                                  
gacgTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACt                                                                                                                                                 <  2:1418351/150‑1 (MQ=255)
 acgTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTg                                                                                                                                                <  1:221231/150‑1 (MQ=255)
   gtgcgtgcGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCg                                                                                                                                              <  1:694242/150‑1 (MQ=255)
      cgtgcgGGGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTa                                                                                                                                           <  1:273817/150‑1 (MQ=255)
        tgcggtgcgCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAAt                                                                                                                                         >  1:766574/1‑150 (MQ=255)
                   tGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAgg                                                                                                                              <  1:1603400/150‑1 (MQ=255)
                    gCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAggg                                                                                                                             <  1:1522687/150‑1 (MQ=255)
                    gCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACCGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAggg                                                                                                                             <  2:1121255/150‑1 (MQ=255)
                        gTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCCCACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAg                                                                                                                         <  1:1392695/150‑1 (MQ=255)
                        gTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAg                                                                                                                         <  1:641418/150‑1 (MQ=255)
                         tGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAgg                                                                                                                        >  2:1557136/1‑150 (MQ=255)
                             ggAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGTTAATCGGCCCCGAGGGCAAGGTGCt                                                                                                                    >  2:1014497/1‑150 (MQ=255)
                                    cgcgACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAg                                                                                                                               <  2:802149/132‑1 (MQ=255)
                                    cgcgACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAg                                                                                                                               >  1:802149/1‑132 (MQ=255)
                                        aCTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCTAGGGGAAGGTGCTGGAGGTGCAGa                                                                                                         <  2:17121/150‑1 (MQ=255)
                                            cTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCg                                                                                                     >  1:45190/1‑150 (MQ=255)
                                             tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCg                                                                                                                                              <  1:535891/108‑1 (MQ=255)
                                             tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCg                                                                                                                                              >  2:535891/1‑108 (MQ=255)
                                             tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGt                                                                                                    >  1:122516/1‑150 (MQ=255)
                                             tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGt                                                                                                    >  2:1200013/1‑150 (MQ=255)
                                             tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCAACCTTGGGAGTGTACCGCTCGTTGCAGACTGCGGTAATGCGCCCCCAGGGCAAGGGGCTGGAGGTGCAGATCCGt                                                                                                    >  1:282389/1‑150 (MQ=255)
                                                acCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTAcc                                                                                                 >  1:1266437/1‑150 (MQ=255)
                                                       tGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCa                                                                                          >  2:900237/1‑150 (MQ=255)
                                                        gggCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCAt                                                                                         <  2:1266437/150‑1 (MQ=255)
                                                               gTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCg                                                                                                                                              >  1:374315/1‑90 (MQ=255)
                                                               gTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCg                                                                                                                                              <  2:374315/90‑1 (MQ=255)
                                                               gTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTg                                                                                                                      <  1:1345682/114‑1 (MQ=255)
                                                               gTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTg                                                                                                                      >  2:1345682/1‑114 (MQ=255)
                                                                 gCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAaga                                                                                >  2:951133/1‑150 (MQ=255)
                                                                  cacaCCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAagag                                                                               >  1:1383113/1‑150 (MQ=255)
                                                                   acacCTTGTGGCGGCATATTCCCAAGGAGTTCGCCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAgg                                                                              >  1:726623/1‑150 (MQ=255)
                                                                   acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAgg                                                                              <  2:766574/150‑1 (MQ=255)
                                                                   acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAgg                                                                              >  1:65065/1‑150 (MQ=255)
                                                                   acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAgg                                                                              >  1:1218751/1‑150 (MQ=255)
                                                                    cacCGTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGATCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTCCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGc                                                                             <  1:380011/150‑1 (MQ=255)
                                                                     acCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGcc                                                                            >  1:573810/1‑150 (MQ=255)
                                                                       cTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGa                                                                          >  1:535289/1‑150 (MQ=255)
                                                                         tgtgGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAAc                                                                        <  2:1255531/150‑1 (MQ=255)
                                                                            ggcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTg                                                                     <  2:573810/150‑1 (MQ=255)
                                                                              cggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                                   <  2:26767/150‑1 (MQ=255)
                                                                              cggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                                   <  2:1593237/150‑1 (MQ=255)
                                                                               ggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCg                                                                  >  1:152858/1‑150 (MQ=255)
                                                                                     ttCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATgcg                                                            >  2:1306590/1‑150 (MQ=255)
                                                                                          aaGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAc                                                       >  2:275579/1‑150 (MQ=255)
                                                                                           aGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACt                                                      <  1:231827/150‑1 (MQ=255)
                                                                                            ggAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTg                                                     <  2:65065/150‑1 (MQ=255)
                                                                                            ggAGTGCGACGACGACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGATGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTg                                                     <  2:652122/150‑1 (MQ=255)
                                                                                            gcaGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTg                                                     <  2:54383/148‑1 (MQ=255)
                                                                                                ttCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGgcgc                                                 <  2:152858/150‑1 (MQ=255)
                                                                                                      gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                           <  2:1218751/150‑1 (MQ=255)
                                                                                                      gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCGGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                           <  2:282389/150‑1 (MQ=255)
                                                                                                        cTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCa                                                                                                           >  2:116055/1‑84 (MQ=255)
                                                                                                        cTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCa                                                                                                           <  1:116055/84‑1 (MQ=255)
                                                                                                        cTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAggg                                         <  2:1383113/150‑1 (MQ=255)
                                                                                                          aCATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCa                                       >  1:748730/1‑150 (MQ=255)
                                                                                                          aCATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAACAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCa                                       <  1:1014497/150‑1 (MQ=255)
                                                                                                           cATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGt                                                                                                              >  1:615061/1‑78 (MQ=255)
                                                                                                           cATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGt                                                                                                              <  2:615061/78‑1 (MQ=255)
                                                                                                           cATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCAc                                      <  1:1200013/150‑1 (MQ=255)
                                                                                                           cATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCAc                                      >  2:1337118/1‑150 (MQ=255)
                                                                                                               gCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATcc                                                                                                      <  2:1073215/82‑1 (MQ=255)
                                                                                                               gCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATcc                                                                                                      >  1:1073215/1‑82 (MQ=255)
                                                                                                               gCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGAc                                  >  2:1552268/1‑150 (MQ=255)
                                                                                                                ccAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGCGAAGCGGCCCCGCGGGCAATGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACg                                 >  2:468765/1‑150 (MQ=255)
                                                                                                                  aaCCCCAAGGAGACCAGCTACAGCTCGTGGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTc                               >  2:555805/1‑150 (MQ=255)
                                                                                                                     cccAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGGCGTCAAg                            >  1:914967/1‑150 (MQ=255)
                                                                                                                      ccAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATACGCCCCCACGGCATGCACGAAGAGGCCGAACTTGAGGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                           >  1:768114/1‑150 (MQ=255)
                                                                                                                      ccAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGAGCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                           >  2:52166/1‑150 (MQ=255)
                                                                                                                        aaGGAGAACGGCTACCGTTCGTTGGAGACGGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGccc                         <  2:122516/150‑1 (MQ=255)
                                                                                                                        aaGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGccc                         >  2:630440/1‑150 (MQ=255)
                                                                                                                           gagaACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCGCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGc                      <  1:1557136/150‑1 (MQ=255)
                                                                                                                                cGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                 >  1:651441/1‑150 (MQ=255)
                                                                                                                                 ggCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAAc                <  2:535289/150‑1 (MQ=255)
                                                                                                                                   cTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                 <  2:962322/147‑1 (MQ=255)
                                                                                                                                   cTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                 >  1:962322/1‑147 (MQ=255)
                                                                                                                                   cTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCa              >  1:604925/1‑150 (MQ=255)
                                                                                                                                    tACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc             >  1:1272572/1‑150 (MQ=255)
                                                                                                                                    tACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc             <  2:45190/150‑1 (MQ=255)
                                                                                                                                    tACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc             >  1:1501128/1‑150 (MQ=255)
                                                                                                                                    tACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACtt                                                                      <  2:516812/93‑1 (MQ=255)
                                                                                                                                    tACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACtt                                                                      >  1:516812/1‑93 (MQ=255)
                                                                                                                                     aCCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc             <  2:5814/149‑1 (MQ=255)
                                                                                                                                     aCCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc             >  1:5814/1‑149 (MQ=255)
                                                                                                                                     aCCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACt            <  1:1306590/150‑1 (MQ=255)
                                                                                                                                     aCCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACt            <  1:1337118/150‑1 (MQ=255)
                                                                                                                                      ccGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTa           >  1:944295/1‑150 (MQ=255)
                                                                                                                                       cGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGc                                                                                                            >  2:1362557/1‑52 (MQ=255)
                                                                                                                                       cGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGc                                                                                                            <  1:1362557/52‑1 (MQ=255)
                                                                                                                                        gCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACg         >  2:510251/1‑150 (MQ=255)
                                                                                                                                        gCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACg         >  1:483056/1‑150 (MQ=255)
                                                                                                                                        gCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACg         >  1:598885/1‑150 (MQ=255)
                                                                                                                                           cGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaag      <  2:483056/150‑1 (MQ=255)
                                                                                                                                            gTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaaga     >  2:1001570/1‑150 (MQ=255)
                                                                                                                                               gCAGACTGCGGGAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGaaaa  <  2:978207/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAA  >  NC_002947/1851812‑1852104

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: ‑