Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,452,047 C→A 5.3% A148A (GCC→GCA) ‡ nagZ → beta‑N‑acetylglucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,452,0470CA 5.3% 1146.9 / 35.4 336A148A (GCC→GCA) ‡nagZbeta‑N‑acetylglucosaminidase
Reads supporting (aligned to +/- strand):  ref base C (169/149);  new base A (8/10);  total (177/159)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 339 total aligned reads displayed.

GAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGCC  >  NC_002947/2451898‑2452196
                                                                                                                                                     |                                                                                                                                                     
   cAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgc                                                                                                                                                    <  2:3290233/150‑1 (MQ=255)
   cAGTACGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGgc                                                                                                                                                    <  2:195204/150‑1 (MQ=255)
      tGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGcgcg                                                                                                                                                 >  2:6407485/1‑150 (MQ=255)
        cGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGtt                                                                                                                                               >  1:1007022/1‑150 (MQ=255)
         ggctggctGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCGGCCAGGCGCGttg                                                                                                                                              <  2:921765/150‑2 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  2:1367292/1‑150 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       <  2:5717759/150‑1 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  2:4234918/1‑150 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  2:5878551/1‑150 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  1:2849568/1‑150 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  1:2559684/1‑150 (MQ=255)
                tgATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTg                                                                                                                                       >  1:2372669/1‑150 (MQ=255)
                 gATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTgg                                                                                                                                      >  2:5931440/1‑150 (MQ=255)
                   tGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCAGCGCGTTCATCCGTGGCa                                                                                                                                    >  1:418839/1‑150 (MQ=255)
                        aCCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAAt                                                                                                                               >  2:3948387/1‑150 (MQ=255)
                          cGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGt                                                                                                                                                <  1:6028155/131‑1 (MQ=255)
                            aGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATgcgg                                                                                                                           >  2:2109537/1‑150 (MQ=255)
                                  tGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTa                                                                                                                     >  1:445820/1‑150 (MQ=255)
                                      ggTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGgc                                                                                                                 <  2:556314/150‑1 (MQ=255)
                                            ccTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtg                                                                                                           >  1:5479307/1‑150 (MQ=255)
                                            ccTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtg                                                                                                           <  1:1821713/150‑1 (MQ=255)
                                              tGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCtgtg                                                                                                         <  1:6741957/150‑1 (MQ=255)
                                                 aCCTCAGCTTCGCTCCGGTGCCCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCa                                                                                                      >  2:4630304/1‑150 (MQ=255)
                                                  ccTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCaa                                                                                                     <  1:2662536/150‑1 (MQ=255)
                                                  ccTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCaa                                                                                                     <  1:4448299/150‑1 (MQ=255)
                                                  ccTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCaa                                                                                                     <  2:5296881/150‑1 (MQ=255)
                                                  ccTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCaa                                                                                                     <  1:3780232/150‑1 (MQ=255)
                                                   cTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGt                                                                                                                                        >  1:3996735/1‑114 (MQ=255)
                                                   cTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGt                                                                                                                                        <  2:3996735/114‑1 (MQ=255)
                                                   cTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAg                                                                                                    >  1:4980088/1‑150 (MQ=255)
                                                    tCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGAACGCGCCTGGCCGGCGCGTTCaa                                                                                                                                            >  1:234726/1‑107 (MQ=255)
                                                    tCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGAACGCGCCTGGCCGGCGCGTTCaa                                                                                                                                            <  2:234726/108‑2 (MQ=255)
                                                       gCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCAtt                                                                                                >  1:1352425/1‑150 (MQ=255)
                                                        cTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCAttt                                                                                               <  1:5713263/150‑1 (MQ=255)
                                                        cTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCAttt                                                                                               >  2:1111379/1‑150 (MQ=255)
                                                        cTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCAttt                                                                                               >  2:4093255/1‑150 (MQ=255)
                                                        cTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCAttt                                                                                               <  1:833552/150‑1 (MQ=255)
                                                            gCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCg                                                                                           >  1:4384823/1‑150 (MQ=255)
                                                              tCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCggg                                                                                         <  1:4629140/150‑1 (MQ=255)
                                                                        gACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTgg                                                                               >  1:5111318/1‑150 (MQ=255)
                                                                        gACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTgg                                                                               >  2:4358207/1‑150 (MQ=255)
                                                                         aCCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTggg                                                                              <  2:5448796/150‑1 (MQ=255)
                                                                         aCCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTggg                                                                              >  1:1525294/1‑150 (MQ=255)
                                                                           cTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGcc                                                                            >  1:5915416/1‑150 (MQ=255)
                                                                           cTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGcc                                                                            >  2:1644445/1‑150 (MQ=255)
                                                                           cTGGACAACCAGCGTAGCGCCGAGGTGGGAAGACGCGCGTTAGAGGGTGACCCGCTGATTGCAAAGCGCCTGTCCGGCGCGTTCATCCGTGGCATGAAAGCGGCGGGTACGGCCGCCGGTGGCAAGCATTTCCCGGGGCATGGCTGGGcc                                                                            >  1:714374/1‑150 (MQ=255)
                                                                              gACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGaa                                                                         <  2:5577333/150‑1 (MQ=255)
                                                                              gACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGaa                                                                         <  2:6202108/150‑1 (MQ=255)
                                                                                ccaccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGc                                                                       <  1:3527424/150‑1 (MQ=255)
                                                                                ccaccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGc                                                                       <  2:6724071/150‑1 (MQ=255)
                                                                                 caccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGcc                                                                      >  2:5027409/1‑150 (MQ=255)
                                                                                 caccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGcc                                                                      >  2:6029375/1‑150 (MQ=255)
                                                                                   ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                                    >  2:3790676/1‑150 (MQ=255)
                                                                                   ccaGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGa                                                                    >  2:3197390/1‑150 (MQ=255)
                                                                                     aGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGAtt                                                                  <  2:2512084/150‑1 (MQ=255)
                                                                                     aGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGAta                                                                  >  1:476431/1‑149 (MQ=255)
                                                                                      gCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGACGCCGATTc                                                                 >  2:4452022/1‑150 (MQ=255)
                                                                                      gCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTc                                                                 <  2:4579201/150‑1 (MQ=255)
                                                                                        gTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGc                                                               <  1:4092967/150‑1 (MQ=255)
                                                                                           gcgcCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACg                                                            <  1:3683242/150‑1 (MQ=255)
                                                                                              ccGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTgg                                                         <  2:1550691/150‑1 (MQ=255)
                                                                                                       gCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCa                                                >  2:3692304/1‑150 (MQ=255)
                                                                                                        cAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCAc                                               <  2:2224291/150‑1 (MQ=255)
                                                                                                        cAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCAc                                               <  1:1439181/150‑1 (MQ=255)
                                                                                                         aGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACt                                              <  1:6816385/150‑1 (MQ=255)
                                                                                                         aGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACt                                              <  1:3095086/150‑1 (MQ=255)
                                                                                                             gcgcgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACtgatg                                          <  1:6351166/150‑1 (MQ=255)
                                                                                                                cgTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAAc                                       <  1:1861919/150‑1 (MQ=255)
                                                                                                                   tCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCa                                    >  2:3364644/1‑150 (MQ=255)
                                                                                                                      aGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGcc                                 <  1:1140945/150‑1 (MQ=255)
                                                                                                                        ggTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTg                               >  2:357837/1‑150 (MQ=255)
                                                                                                                        ggTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTg                               <  2:2651327/150‑1 (MQ=255)
                                                                                                                        ggTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTg                               >  2:1428324/1‑150 (MQ=255)
                                                                                                                          tGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGa                             <  2:775638/150‑1 (MQ=255)
                                                                                                                           gACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGc                                                               <  2:3461966/115‑1 (MQ=255)
                                                                                                                           gACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGc                                                               >  1:3461966/1‑115 (MQ=255)
                                                                                                                           gACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaa                            <  1:3386790/150‑1 (MQ=255)
                                                                                                                              ccGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGaacaa                         >  1:3824607/1‑150 (MQ=255)
                                                                                                                                gctgctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACt                       <  1:1644445/150‑1 (MQ=255)
                                                                                                                                  tgctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgc                     <  1:2079496/150‑1 (MQ=255)
                                                                                                                                    ctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgcgc                   <  2:3496599/150‑1 (MQ=255)
                                                                                                                                    ctTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTgcgc                   >  1:3027891/1‑150 (MQ=255)
                                                                                                                                     ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGcc                  >  2:2924173/1‑150 (MQ=255)
                                                                                                                                     ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGcc                  >  1:4151041/1‑150 (MQ=255)
                                                                                                                                     ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGcc                  >  2:1205005/1‑150 (MQ=255)
                                                                                                                                     ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGcc                  >  1:92248/1‑150 (MQ=255)
                                                                                                                                     ttGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCCGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGcc                  >  2:4893702/1‑150 (MQ=255)
                                                                                                                                      tGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCa                 >  1:4312735/1‑150 (MQ=255)
                                                                                                                                      tGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCa                 <  1:974848/150‑1 (MQ=255)
                                                                                                                                      atgaaCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGccc                 >  1:226936/4‑149 (MQ=255)
                                                                                                                                       gCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAg                >  1:3295579/1‑150 (MQ=255)
                                                                                                                                         aaCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCCGCCTGGaac                           <  2:1798526/137‑1 (MQ=255)
                                                                                                                                          aCGCGCCGGCCAGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTcc             >  1:5070438/1‑150 (MQ=255)
                                                                                                                                           cgcgCCGGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCg            <  1:6596542/150‑1 (MQ=255)
                                                                                                                                                cTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTg       <  2:2849568/150‑1 (MQ=255)
                                                                                                                                                  ggCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGt     <  1:3222105/150‑1 (MQ=255)
                                                                                                                                                  ggCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGt     <  2:4812641/150‑1 (MQ=255)
                                                                                                                                                   gCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCTGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTg    <  1:217225/150‑1 (MQ=255)
                                                                                                                                                   gCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGTCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTg    <  1:5627678/150‑1 (MQ=255)
                                                                                                                                                     cGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGcc  <  2:6544077/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCAGTGCGGCTGGCTGATGGCGACCGAGGTGCTGGCGGTTGGCCTGGACCTCAGCTTCGCTCCGGTGCTCGACCTGGACCACCAGCGTAGCGCCGTGGTGGGCAGCCGCGCGTTCGAGGGTGACCCGCTGCTTGCAACGCGCCTGGCCGGCGCGTTCATCCGTGGCATGAATGCGGCGGGTATGGCCGCCTGTGGCAAGCATTTCCCGGGGCATGGCTGGGCCGAAGCCGATTCGCACGTGGCCATCCCCACTGATGAACGCAGCCTGGAACAACTGCGCCAGTCCGACCTGGTGCC  >  NC_002947/2451898‑2452196

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