Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,447,098:1 +A 10.0% coding (444/1011 nt) nagZ → beta‑N‑acetylglucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,447,0981.A10.0% 245.9 / 23.9 70coding (444/1011 nt)nagZbeta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base . (36/27);  new base A (5/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

CAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGCC  >  NC_002947/2446952‑2447247
                                                                                                                                                   |                                                                                                                                                     
cAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgc                                                                                                                                                    <  1:1138674/150‑1 (MQ=255)
cAGTGCGGCTGGCTGATGGCGACCGAGGTGACGGTGTTTTGCCTGGACCTCAGCTTCGGTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgc                                                                                                                                                    <  2:564855/150‑1 (MQ=255)
         tggctgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATc                                                                                                                                           >  2:603655/1‑150 (MQ=255)
         tggctgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATc                                                                                                                                           >  2:528989/1‑150 (MQ=255)
             tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  2:2120443/1‑150 (MQ=255)
                 ggCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCAt                                                                                                                                   <  1:2089053/150‑1 (MQ=255)
                  gCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATg                                                                                                                                  >  2:2141884/1‑150 (MQ=255)
                     aCCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTGCGGTCCGCGACCTGGACCACCGTCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGCCAtcact                                                                                                                               >  2:694354/1‑146 (MQ=255)
                     aCCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGTGCGTTCATCCGTGGCATGAAt                                                                                                                               >  2:440044/1‑150 (MQ=255)
                     aCCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAAt                                                                                                                               <  2:1448707/150‑1 (MQ=255)
                      ccGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATg                                                                                                                              >  2:1576560/1‑150 (MQ=255)
                           gTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCGCAGCGGGGGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTAATCCGTAGCATTAATgcggcg                                                                                                                         >  1:1332281/1‑150 (MQ=255)
                            tGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGcggcgg                                                                                                                        <  1:1053847/150‑1 (MQ=255)
                            tGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGcggcgg                                                                                                                        <  1:434308/150‑1 (MQ=255)
                                 gCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATg                                                                                                                   <  1:1792286/150‑1 (MQ=255)
                                        gCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCAAGTCCGGCCTGCGGCGTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCt                                                                                                            >  2:1380633/1‑150 (MQ=255)
                                          cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtgt                                                                                                          >  2:594071/1‑150 (MQ=255)
                                          cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtgt                                                                                                          >  2:737448/1‑150 (MQ=255)
                                          cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                            >  1:1894851/1‑113 (MQ=255)
                                          cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                            >  1:1899420/1‑113 (MQ=255)
                                          cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                            <  2:1899420/116‑4 (MQ=255)
                                          cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                            <  2:1894851/116‑4 (MQ=255)
                                           tGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtgtg                                                                                                         >  2:290730/1‑150 (MQ=255)
                                            ggACCTCAGCTTCGTCCCGGTGCTCGACCTGGACCACCAGCGTAGACCCCTGGTGGGCAGCCGCGCGTTCGAGGGTGCCCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTgg                                                                                                        <  1:1453957/150‑1 (MQ=255)
                                                  cAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCa                                                                                                  >  2:1989677/1‑150 (MQ=255)
                                                     cTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCAttt                                                                                               >  1:1457606/1‑150 (MQ=255)
                                                      ttCGCTGCGGGGCGCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTc                                                                                              >  2:1172732/1‑150 (MQ=255)
                                                          cTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCgg                                                                                          <  1:737448/150‑1 (MQ=255)
                                                          cTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGAAGCCGCGCGTTCGAGGTTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCgg                                                                                          <  2:2099467/150‑1 (MQ=255)
                                                           tCCGGTGCTCGACCTGGACCACCGGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGACCGCCTGTGGCAAGCATTTCACggg                                                                                         >  2:1191040/1‑150 (MQ=255)
                                                            ccGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCgggg                                                                                        <  1:2141884/150‑1 (MQ=255)
                                                                   tCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCt                                                                                 >  2:2379944/1‑150 (MQ=255)
                                                                       ccTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGc                                                                             >  1:1246941/1‑150 (MQ=255)
                                                                       ccTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGAGCATGCGTGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGc                                                                             <  1:2120443/150‑1 (MQ=255)
                                                                        cTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGcc                                                                            <  1:1576560/150‑1 (MQ=255)
                                                                          ggACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGa                                                                          >  2:2289176/1‑150 (MQ=255)
                                                                          ggACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGa                                                                          >  1:879119/1‑150 (MQ=255)
                                                                           gACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTCGACGGCAGGTTCATCCGTGGTGGGAAGGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGaa                                                                         <  1:290730/150‑1 (MQ=255)
                                                                                ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                                    >  2:104460/1‑150 (MQ=255)
                                                                                ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                                    >  1:2326967/1‑150 (MQ=255)
                                                                                ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                                    >  2:2041138/1‑150 (MQ=255)
                                                                                 caGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGAt                                                                   >  1:1212562/1‑150 (MQ=255)
                                                                                 caGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGAt                                                                   >  1:211505/1‑150 (MQ=255)
                                                                                  aGCGTAGCGCCGTGGTGGGCAGCCGGGACTGCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGAtt                                                                  <  1:603655/150‑1 (MQ=255)
                                                                                           ccGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTgg                                                         >  1:75422/1‑150 (MQ=255)
                                                                                              tggtggGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTGGCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCa                                                      >  1:391705/1‑150 (MQ=255)
                                                                                                       gCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGTGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtg                                             <  1:440044/150‑1 (MQ=255)
                                                                                                        ccGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtga                                            <  2:1457606/150‑1 (MQ=255)
                                                                                                         cgcgcgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtgat                                           <  1:2041138/150‑1 (MQ=255)
                                                                                                             cgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAAc                                       <  1:2257466/150‑1 (MQ=255)
                                                                                                                tCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCa                                    >  2:319607/1‑150 (MQ=255)
                                                                                                                  gAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGc                                  >  2:714198/1‑150 (MQ=255)
                                                                                                                     ggTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTg                               >  1:1854247/1‑150 (MQ=255)
                                                                                                                        gACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaa                            <  1:1989677/150‑1 (MQ=255)
                                                                                                                          cccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaaca                          >  2:2129714/1‑150 (MQ=255)
                                                                                                                           ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGTGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                         >  2:462096/1‑150 (MQ=255)
                                                                                                                           ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                         >  1:2094405/1‑150 (MQ=255)
                                                                                                                           ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                         >  2:914185/1‑150 (MQ=255)
                                                                                                                               tgctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgc                     <  1:1172732/150‑1 (MQ=255)
                                                                                                                                 ctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgcgc                   <  1:2379944/150‑1 (MQ=255)
                                                                                                                                 ctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgcgc                   <  1:2329682/150‑1 (MQ=255)
                                                                                                                                  ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTgg                                                                                                        <  1:217391/64‑1 (MQ=255)
                                                                                                                                  ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTgg                                                                                                        >  2:217391/1‑64 (MQ=255)
                                                                                                                                   ggCAACTCGCCTGGCCGGCGCCTTCATCAGTGGCATGCATGCGGTGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCa                 <  2:391705/149‑1 (MQ=255)
                                                                                                                                    gCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAg                >  2:467435/1‑150 (MQ=255)
                                                                                                                                     gaaCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGt               >  1:1597095/2‑150 (MQ=255)
                                                                                                                                       aCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTcc             >  2:1464382/1‑150 (MQ=255)
                                                                                                                                       aCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTcc             >  1:2450974/1‑150 (MQ=255)
                                                                                                                                               ggCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGt     <  2:2326967/150‑1 (MQ=255)
                                                                                                                                                gCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTg    <  1:528989/150‑1 (MQ=255)
                                                                                                                                                 ccGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGc   <  2:1854247/150‑1 (MQ=255)
                                                                                                                                                  cGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGcc  >  1:613676/1‑150 (MQ=255)
                                                                                                                                                  cGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGcc  <  2:75422/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGCC  >  NC_002947/2446952‑2447247

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: