Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,812,781	G→A	12.3%	intergenic (+32/‑78)	rplT → / → pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,812,781	0	G	A	12.3%	157.9 / 10.9	64	intergenic (+32/‑78)	rplT/pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (4/4);  total (30/34)
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.82e-01

CTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCC  >  NC_002947/2812634‑2812928                                                                                                                                                    |                                                                                                                                                     cTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCt                                                                                                                                                    >  2:35338/1‑150 (MQ=255)    gCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGcc                                                                                                                                                 >  1:803953/1‑150 (MQ=255)      cTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCgg                                                                                                                                               >  1:389056/1‑150 (MQ=255)      cTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCgg                                                                                                                                               >  1:1238574/1‑150 (MQ=255)        gAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCgg                                                                                                                                               <  2:1307421/148‑1 (MQ=255)        gAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCgg                                                                                                                                               >  1:1307421/1‑148 (MQ=255)              ggCTTCGATCGAAGTCGACCGTAAGGTTCTGGCTGATCGGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGaa                                                                                                                                       <  1:1891220/150‑1 (MQ=255)               gCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                                      <  2:1386913/150‑1 (MQ=255)                cTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGTGCTGCCGGGCATTGAACg                                                                                                                                     <  1:913107/150‑1 (MQ=255)                      tcgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                                      >  2:1733837/1‑143 (MQ=255)                      tcgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                                      <  1:1733837/143‑1 (MQ=255)                           aTCGACCGTCAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGCCGGGAAGACTAAAGACGTTCTGGCTTAAGTACCCGCGACACTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAgg                                                                                                                          <  1:611419/150‑1 (MQ=255)                            tCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAggg                                                                                                                         <  2:1811455/150‑1 (MQ=255)                                 cGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGTGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAaga                                                                                                                    >  2:642893/1‑150 (MQ=255)                                 cGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAaga                                                                                                                    <  2:389056/150‑1 (MQ=255)                                 cGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAaga                                                                                                                    <  2:1914539/150‑1 (MQ=255)                                 cGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgg                                                                                                                                               >  1:1759039/1‑123 (MQ=255)                                 cGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgg                                                                                                                                               <  2:1759039/123‑1 (MQ=255)                                    aaGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGcc                                                                                                                 >  1:1098543/1‑150 (MQ=255)                                      ggTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCtt                                                                                                               <  1:1246101/150‑1 (MQ=255)                                       gTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTATGGCTTAAGTACCCGACCCAGTGACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTc                                                                                                              <  2:731420/150‑1 (MQ=255)                                         tCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCaa                                                                                                            <  2:1317487/150‑1 (MQ=255)                                            ggCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGc                                                                                                         >  2:634634/1‑150 (MQ=255)                                                gATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGACTTCCTTGTTCGCTACAATCTCTCGGCGGTGCTAGCGCGGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCtt                                                                                                     <  1:118454/150‑1 (MQ=255)                                                   cTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTccc                                                                                                  <  2:1238574/150‑1 (MQ=255)                                                     ggCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCt                                                                                                >  2:1793768/1‑150 (MQ=255)                                                            aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAg                                                                                                                           >  2:162937/1‑116 (MQ=255)                                                            aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAg                                                                                                                           <  1:162937/116‑1 (MQ=255)                                                               gAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTAt                                                                                      <  2:803953/150‑1 (MQ=255)                                                                 aaaaGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                                    >  1:252977/1‑150 (MQ=255)                                                                     gcggcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAg                                                                                >  1:617044/1‑150 (MQ=255)                                                                      cggcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGNAGt                                                                               <  2:1187019/150‑1 (MQ=255)                                                                           tttGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTa                                                                          <  2:617044/150‑1 (MQ=255)                                                                                tgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGa                                                                                 <  1:1137236/138‑1 (MQ=255)                                                                                tgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGa                                                                                 >  2:1137236/1‑138 (MQ=255)                                                                                   gATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGACAGTCATCGATAGGGGAAGAGCCTTCACAGCTCGTCCCCAATTTTCGTGTCTGGAAGCTGCACATGGcaa                                                                  >  1:205270/1‑150 (MQ=255)                                                                                          gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGa                                                          >  2:854638/1‑150 (MQ=255)                                                                                          gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGa                                                          >  1:333132/1‑150 (MQ=255)                                                                                               aGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATgcgc                                                     >  1:884529/1‑150 (MQ=255)                                                                                                  tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                                      >  2:325978/1‑67 (MQ=255)                                                                                                  tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                                      <  1:325978/67‑1 (MQ=255)                                                                                                  tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgggtga                                                                                                                                           <  1:65030/62‑4 (MQ=18)                                                                                                  tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgggtga                                                                                                                                           >  2:65030/1‑59 (MQ=18)                                                                                                   aaaGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGt                                                 <  2:1079817/150‑1 (MQ=255)                                                                                                     aGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGtct                                               >  1:1808714/1‑150 (MQ=255)                                                                                                       ccGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTTAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTcc                                             <  2:919651/150‑1 (MQ=255)                                                                                                        cGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTccc                                            <  1:854638/150‑1 (MQ=255)                                                                                                               gCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCt                                     >  2:1508427/1‑150 (MQ=255)                                                                                                                  tAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTaga                                  <  2:1767596/150‑1 (MQ=255)                                                                                                                   aaGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTagag                                 <  2:1294993/150‑1 (MQ=255)                                                                                                                     gTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGc                               <  1:634634/150‑1 (MQ=255)                                                                                                                        cccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGt                            <  1:1793768/150‑1 (MQ=255)                                                                                                                         ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTcccc                                                                                                 >  2:969716/1‑81 (MQ=255)                                                                                                                         ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTcccc                                                                                                 <  1:969716/81‑1 (MQ=255)                                                                                                                           gcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGa                         >  1:1248263/1‑150 (MQ=255)                                                                                                                                aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg                    >  2:599187/1‑150 (MQ=255)                                                                                                                                aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg                    >  2:1377628/1‑150 (MQ=255)                                                                                                                                    aCCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTCCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCCAGGTGGCGCATG‑GTGCTAATGCCTAGAGGCCGTGGAACGCGCTGa                >  1:1866924/1‑150 (MQ=255)                                                                                                                                     cccGGCGGTGCTCGCGCTGCCGGGCATTGAACGTCATCGATAGGGTAAAAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCGAAGCCCTAGAGGCCGTGGAACGCGCTGaa               <  1:108456/150‑1 (MQ=255)                                                                                                                                     cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa               <  2:1574206/150‑1 (MQ=255)                                                                                                                                     cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa               <  2:1574040/150‑1 (MQ=255)                                                                                                                                     cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa               >  1:1574206/1‑150 (MQ=255)                                                                                                                                     cccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa               >  1:1574040/1‑150 (MQ=255)                                                                                                                                                 aGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGcc   <  1:1377628/150‑1 (MQ=255)                                                                                                                                                  gcgcTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGccc  <  1:642893/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAA‑CCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCC  >  NC_002947/2812634‑2812928

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