Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,282,920 G→A 5.3% intergenic (+83/+55) trmA → / ← pcaG tRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,282,9200GA 5.3% 350.0 / 5.6 114intergenic (+83/+55)trmA/pcaGtRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (57/49);  new base A (3/3);  total (60/54)
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.89e-01

Only 100 of 117 total aligned reads displayed.

CATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGC  >  NC_002947/5282773‑5283064
                                                                                                                                                   |                                                                                                                                                
cATTGTTTGACCTGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCTGCGCTTATCGTTGTATTGCCGGGGCTGCTTTTCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCTGATCGtt                                                                                                                                                <  1:4496161/150‑1 (MQ=255)
cATTGTTTGACCAGGTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTAGTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                                                <  1:2138680/150‑1 (MQ=255)
 aTTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTg                                                                                                                                               <  2:2079646/150‑1 (MQ=255)
  ttGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACAATGATCGGCGTCGGTGTAAGCAGCTCGgtgt                                                                                                                                              >  1:1374452/1‑150 (MQ=255)
        gACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGCTCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGc                                                                                                                                        <  2:1985723/150‑1 (MQ=255)
         aCCAGTTCCAGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCg                                                                                                                                       <  2:2373829/150‑1 (MQ=255)
            aggtCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTTGGTGTACGCAGATCGTTGTAGGAGCGGGc                                                                                                                                    <  1:2907521/147‑1 (MQ=255)
            aGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGc                                                                                                                                    <  2:2542570/150‑1 (MQ=255)
             ggtCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCt                                                                                                                                   <  2:883535/148‑1 (MQ=255)
             gTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCt                                                                                                                                   <  1:3951945/150‑1 (MQ=255)
              ttCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCtt                                                                                                                                  <  1:3235635/150‑1 (MQ=17)
                cccGTAGACCCACCATAGGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTAGTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGcaagc                                                                                                                                <  2:1214859/150‑5 (MQ=255)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCTTg                                                                                                                                 <  1:1292303/146‑1 (MQ=17)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCTTg                                                                                                                                 >  2:1292303/1‑146 (MQ=17)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGt                                                                                                                                              >  1:1164280/1‑133 (MQ=255)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGt                                                                                                                                              <  2:1164280/133‑1 (MQ=255)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCggg                                                                                                                                     <  1:5134012/142‑1 (MQ=255)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCggg                                                                                                                                     >  2:5134012/1‑142 (MQ=255)
                   gTATACCCACCATATGGAGAGTGGGGTTCTGCGGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCcgc                                                                                                                            <  2:1930884/150‑1 (MQ=255)
                      taCCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCtt                                                                                                                                  <  1:542603/142‑1 (MQ=17)
                      taCCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCtt                                                                                                                                  >  2:542603/1‑142 (MQ=17)
                         ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATAGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATg                                                                                                                       >  1:2950716/1‑150 (MQ=255)
                         ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCCCTCCCACAATTGTCGGCGTCGGTGTAAGCAGATCGTTGTTGGTGCGGGCGTGCCCCCGAATg                                                                                                                       >  1:1505552/1‑150 (MQ=38)
                         ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCCCACCGACAATGATAGGCGGCGGTGTACCCAGATAGGTGTAGGAGCGGGCGTGCCCGCGaaag                                                                                                                       >  1:219917/1‑148 (MQ=38)
                          caccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATgg                                                                                                                      >  2:577654/1‑150 (MQ=255)
                          caccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATgg                                                                                                                      >  2:2599786/1‑150 (MQ=255)
                          caccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCCCGAATgg                                                                                                                      >  2:577628/1‑150 (MQ=255)
                           accaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGCCATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATggg                                                                                                                     <  2:5271093/150‑1 (MQ=255)
                            ccaTATGGAGCGGGGGGTACGGCGGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATCCGGGGGCACGCCCGCTCCGACTTTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGTGCGGGCTTGCCCGCGAATGGGc                                                                                                                    <  2:5006423/150‑1 (MQ=255)
                            ccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGCTCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGc                                                                                                                    <  2:4949607/150‑1 (MQ=255)
                                atGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCa                                                                                                                >  1:4022706/1‑150 (MQ=255)
                                      agTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCgg                                                                                                          >  2:1489357/1‑150 (MQ=255)
                                        tGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGcc                                                                                                        >  2:5390540/1‑150 (MQ=255)
                                        tGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGcc                                                                                                        >  2:3985270/1‑150 (MQ=255)
                                        tGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGcc                                                                                                        <  2:420257/150‑1 (MQ=255)
                                        tGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGcc                                                                                                        >  2:1988335/1‑150 (MQ=255)
                                        tGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCCCACCGACAATGATCGGCGTCGGTGTACGCACATCGTTGTAGGAGCGGGCTTGCCCCCCAAATGGCCGCAGAGCGGcc                                                                                                        >  1:2795611/1‑150 (MQ=38)
                                          gggTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGACGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGcccc                                                                                                      >  2:4510716/1‑150 (MQ=255)
                                               ctgctgGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCt                                                                                                 >  1:4931438/1‑150 (MQ=255)
                                                  ctgGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCa                                                                                              >  2:4458308/1‑150 (MQ=255)
                                                   tgGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGTTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCAc                                                                                             <  2:3307248/150‑1 (MQ=255)
                                                   tgGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGTTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCAc                                                                                             <  2:3307218/150‑1 (MQ=255)
                                                                ctCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGCCATTGATCGGCGTCGGTGTACGCAGATCTTTGTAGGTGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCgaa                                                                                <  1:3102355/149‑1 (MQ=255)
                                                                  cGTTTTATTGCCGGGGCTGCTTTGCAGCCAATGCGCGGGCAAGCCCGCTCCGACATTGAACGGCGTCGGTTTACGAAGATCGTTGTAGGTGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCgaaga                                                                              <  2:1815201/150‑1 (MQ=38)
                                                                  cGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCgaaga                                                                              >  2:3625329/1‑150 (MQ=255)
                                                                    ttGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGaagaag                                                                            >  1:4047802/1‑150 (MQ=255)
                                                                    ttGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCCCACCGACCGTGATCGGCGGCGGTGTACCCCCCTCGGGGTAGGAGCGGGCGTGCCCGCGCAAGGGCCGCAGAGCCGCCCCCAGCGTCACCACAGAAGTCGcagaag                                                                            >  1:1176076/1‑150 (MQ=25)
                                                                        aTTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGTCgg                                                                        >  1:1322498/1‑150 (MQ=255)
                                                                         ttGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGCCATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGt                                                                       >  1:3165305/1‑150 (MQ=255)
                                                                         ttGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGt                                                                       >  1:3165754/1‑150 (MQ=255)
                                                                            cgggggCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGtctc                                                                    <  2:3238617/148‑1 (MQ=255)
                                                                            ccGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGCCAGTGATCGGCGTCGGTGTACGCCGATGGTTGTAGGCGCGGGCGTGCCCGCGAATGGGCCGCAGAGCGGCCCCCAGCGTCACCACAGAAGGCGGAGAAGAAGGtctc                                                                    >  1:4798128/1‑150 (MQ=37)
                                                                            ccGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGtctc                                                                    <  2:3778961/150‑1 (MQ=255)
                                                                            ccGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCACCGACATTTATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGtctc                                                                    >  2:158038/1‑150 (MQ=255)
                                                                                    gctTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTccc                                                            <  2:782844/150‑1 (MQ=255)
                                                                                       ttGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCt                                                         <  2:390662/150‑1 (MQ=255)
                                                                                             ccctgtCGCGGGCACGCCCGCTCCGACATGGCTCGGCGTGGGTGTCCGCTGATGGTTGTTGTAGGGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGc                                                   <  1:4316927/145‑1 (MQ=255)
                                                                                               cATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCgg                                                 <  2:4791957/150‑1 (MQ=255)
                                                                                                aTTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGa                                                <  2:2826014/150‑1 (MQ=255)
                                                                                                aTTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGa                                                <  2:4000524/150‑1 (MQ=255)
                                                                                                aTTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGa                                                >  1:2609513/1‑150 (MQ=255)
                                                                                                  tCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTGGGTGTACGCAGATCGTTGTTGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGata                                              <  2:680935/150‑1 (MQ=255)
                                                                                                    gcgGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATc                                            >  2:922936/1‑150 (MQ=255)
                                                                                                    gcgGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATc                                            >  1:5423565/1‑150 (MQ=255)
                                                                                                      gggCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGACGAAGACGGTCTCCCCTTCCCCCTGAAGGCGGATATCaa                                          >  1:4452008/1‑150 (MQ=14)
                                                                                                        gCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAg                                        <  1:3545937/150‑1 (MQ=255)
                                                                                                        gCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAg                                        >  1:971771/1‑150 (MQ=255)
                                                                                                          aaGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCCAAGCg                                      >  2:1661773/1‑150 (MQ=255)
                                                                                                          aaGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCg                                      >  2:4752293/1‑150 (MQ=255)
                                                                                                          aaGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCg                                      >  2:4414855/1‑150 (MQ=255)
                                                                                                          aaGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCg                                      >  2:1661759/1‑150 (MQ=255)
                                                                                                           aGCCCGCTCCGCCATTGATCGGCGTCGGTGTACGCAGAGCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCgg                                     <  2:158435/150‑1 (MQ=255)
                                                                                                           aGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCgg                                     >  1:4684330/1‑150 (MQ=255)
                                                                                                           aGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCgg                                     >  1:1945045/1‑150 (MQ=255)
                                                                                                              ccGCTCCCACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATAACAAAGCGGTAc                                  <  2:797777/150‑1 (MQ=255)
                                                                                                                gCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGaa                                                                                      >  2:4587241/1‑96 (MQ=255)
                                                                                                                gCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGaa                                                                                      <  1:4587241/96‑1 (MQ=255)
                                                                                                                    cGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTc                            >  2:114104/1‑150 (MQ=255)
                                                                                                                    cGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTc                            >  2:3704969/1‑150 (MQ=255)
                                                                                                                      aCATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGc                                                   <  1:1871542/125‑1 (MQ=255)
                                                                                                                      aCATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTACCCCCTGAATGc                                                   >  2:1871542/1‑125 (MQ=255)
                                                                                                                       cttggctCGGGGTGGGTTTACGCTGTTCGTGTTTGGCGGGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTc                         <  1:982228/144‑1 (MQ=255)
                                                                                                                        aTTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTcc                        >  1:1055253/1‑150 (MQ=255)
                                                                                                                             tCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATcc                   >  1:3216473/1‑150 (MQ=255)
                                                                                                                              cGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGc                                                                                                         >  1:3236473/1‑63 (MQ=255)
                                                                                                                              cGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGc                                                                                                         <  2:3236473/63‑1 (MQ=255)
                                                                                                                                gCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAATAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACt                >  1:4445360/1‑150 (MQ=255)
                                                                                                                                  gTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTc              >  1:3052590/1‑150 (MQ=255)
                                                                                                                                  gTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTc              >  1:3052130/1‑150 (MQ=255)
                                                                                                                                      gtgtACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCgg                                                 <  1:5287310/111‑1 (MQ=255)
                                                                                                                                      gtgtACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCgg                                                 >  2:5287310/1‑111 (MQ=255)
                                                                                                                                            gCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTg    >  2:2473783/1‑150 (MQ=255)
                                                                                                                                             ctgATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTgg   <  2:219917/148‑1 (MQ=255)
                                                                                                                                             cctttggTTGTAGGCGGGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTgg   <  1:577628/144‑1 (MQ=255)
                                                                                                                                             cccaTCGTTTTCGGCGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTgg   <  1:577654/147‑1 (MQ=255)
                                                                                                                                             cAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTgg   <  2:2950716/150‑1 (MQ=255)
                                                                                                                                             cAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTgg   <  2:1945045/150‑1 (MQ=255)
                                                                                                                                              aGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGc  >  2:4378679/1‑150 (MQ=255)
                                                                                                                                                   gTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATcc                   <  2:4959137/128‑1 (MQ=255)
                                                                                                                                                   gTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATcc                   >  1:4959137/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGC  >  NC_002947/5282773‑5283064

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