Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,447,094 Δ1 bp 10.1% coding (440/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,0940T.10.1% 198.7 / 17.6 69coding (440/1011 nt)nagZbeta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base T (36/26);  new base . (5/2);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GGCTGAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCC  >  NC_002947/2446945‑2447236
                                                                                                                                                     |                                                                                                                                               
ggCTGAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTTGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCt                                                                                                                                                 <  2:589745/150‑1 (MQ=255)
   tGAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGc                                                                                                                                              <  1:40099/150‑1 (MQ=255)
       cAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgc                                                                                                                                         <  1:1138674/150‑1 (MQ=255)
       cAGTGCGGCTGGCTGATGGCGACCGAGGTGACGGTGTTTTGCCTGGACCTCAGCTTCGGTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgc                                                                                                                                         <  2:564855/150‑1 (MQ=255)
                tggctgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATc                                                                                                                                >  2:528989/1‑150 (MQ=255)
                tggctgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATc                                                                                                                                >  2:603655/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:434308/150‑1 (MQ=255)
                                   tGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGcggcgg                                                                                                             <  1:1053847/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                                                                                                                                 <  2:1894851/116‑4 (MQ=255)
                                                 cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                 <  2:1899420/116‑4 (MQ=255)
                                                 cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                 >  1:1899420/1‑113 (MQ=255)
                                                 cTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttgca                                                                                                                                 >  1:1894851/1‑113 (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:2041138/1‑150 (MQ=255)
                                                                                       ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                         >  2:104460/1‑150 (MQ=255)
                                                                                       ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                         >  1:2326967/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                                                                                             >  2:217391/1‑64 (MQ=255)
                                                                                                                                         ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTgg                                                                                             <  1:217391/64‑1 (MQ=255)
                                                                                                                                          ggCAACTCGCCTGGCCGGCGCCTTCATCAGTGGCATGCATGCGGTGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCa      <  2:391705/149‑1 (MQ=255)
                                                                                                                                           gCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAg     >  2:467435/1‑150 (MQ=255)
                                                                                                                                            gaaCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGt    >  1:1597095/2‑150 (MQ=255)
                                                                                                                                              aCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTcc  >  1:2450974/1‑150 (MQ=255)
                                                                                                                                              aCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTcc  >  2:1464382/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GGCTGAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCC  >  NC_002947/2446945‑2447236

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: