Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,243	C→G	6.3%	intergenic (‑137/+39)	PP_4606 ← / ← PP_4607	ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,229,243	0	C	G	6.3%	200.4 / 4.7	64	intergenic (‑137/+39)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base G (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGC  >  NC_002947/5229095‑5229385                                                                                                                                                     |                                                                                                                                                 ggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                                                 >  2:1193592/1‑150 (MQ=255)  gTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTCCCTGGAAACATTTTCGCTGTGGCTGCCTCTTCGCGGGTAAACCCGCCCCTa                                                                                                                                                >  2:919943/1‑150 (MQ=255)        ttCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGAGTTTCCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                         <  1:919943/150‑1 (MQ=255)             ggACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTa                                                                                                                                    <  2:140487/150‑1 (MQ=255)              gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAc                                                                                                                                   >  1:32470/1‑150 (MQ=255)               aCGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAcc                                                                                                                                  >  1:1630520/1‑150 (MQ=255)                  cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGt                                                                                                                               <  1:1193592/150‑1 (MQ=255)                       cTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGTCGGCCTCGTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAaga                                                                                                                          >  2:100273/1‑150 (MQ=255)                       cTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAaga                                                                                                                          >  2:755714/1‑150 (MQ=255)                         gtgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGGTGCCACGTGCGGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTAAAGAgg                                                                                                                        >  1:1604956/1‑150 (MQ=255)                                  tGATTATGGCGAGTCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCagag                                                                                                               <  1:337593/150‑1 (MQ=255)                                   gATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              >  1:942040/1‑150 (MQ=255)                                   gATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                              >  1:606355/1‑150 (MQ=255)                                       attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAAc                                                                                                          >  1:301601/1‑150 (MQ=255)                                          gCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTAc                                                                                                       <  2:606355/150‑1 (MQ=255)                                              atcaatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑aa                                                                                                  >  2:1098073/1‑150 (MQ=255)                                                   tcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCCCCTGTGCCTGGAAATCTCCTTGCTGTTGCTGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTGCCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAgg                                                                                             <  1:851511/150‑1 (MQ=255)                                                   tcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAgg                                                                                             >  1:742461/1‑150 (MQ=255)                                                             aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcaca                                                                                                                <  2:376109/122‑1 (MQ=255)                                                             aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcaca                                                                                                                >  1:376109/1‑122 (MQ=255)                                                               tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCGACGGG‑GGGTTTACCCGTGAAGAGGCTGGAACAGGCAACTACCG‑AACCAGGGTCAAACGcccc                                                                                 <  2:679973/150‑1 (MQ=255)                                                                      acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACgg                                                                          >  2:471484/1‑150 (MQ=255)                                                                          aTGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTg                                                                      >  1:1329863/1‑150 (MQ=255)                                                                           tGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGc                                                                     <  2:52606/150‑1 (MQ=255)                                                                                    cTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaaaa                                                            <  2:1329863/150‑1 (MQ=255)                                                                                     tGTTACCTCACCTGTTCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaa                                                               <  1:164818/146‑1 (MQ=255)                                                                                     tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaa                                                               >  2:164818/1‑146 (MQ=255)                                                                                     tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                                           <  2:742461/150‑1 (MQ=255)                                                                                      gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgg                                                                                                                                             <  1:1215846/69‑1 (MQ=255)                                                                                      gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgg                                                                                                                                             >  2:1215846/1‑69 (MQ=255)                                                                                      gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                          <  1:755714/150‑1 (MQ=255)                                                                                      gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                          >  1:866479/1‑150 (MQ=255)                                                                                       ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAccc                                                         >  1:709329/1‑150 (MQ=255)                                                                                          ccTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                                 <  2:52277/139‑1 (MQ=255)                                                                                          ccTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTg                                                                 >  1:52277/1‑139 (MQ=255)                                                                                            tCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCa                                                                                               <  2:601941/107‑1 (MQ=255)                                                                                            tCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCa                                                                                               >  1:601941/1‑107 (MQ=255)                                                                                              aCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                                  >  1:8883/1‑150 (MQ=255)                                                                                                 tgtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACg                                                                           >  2:944969/1‑122 (MQ=255)                                                                                                 tgtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACg                                                                           <  1:944969/122‑1 (MQ=255)                                                                                                      cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                          <  2:709329/150‑1 (MQ=255)                                                                                                       tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAACAGGCGGGCCCAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                                         >  1:667638/1‑150 (MQ=255)                                                                                                                 tttGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcgc                               <  2:301601/150‑1 (MQ=255)                                                                                                                        tgGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCagcagc                        >  1:781054/1‑150 (MQ=255)                                                                                                                           cggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGc                     <  2:866479/150‑1 (MQ=255)                                                                                                                             gcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTc                   <  2:942040/150‑1 (MQ=255)                                                                                                                             gcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTc                   >  1:132428/1‑150 (MQ=255)                                                                                                                              ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  <  2:1630520/150‑1 (MQ=255)                                                                                                                              ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  >  2:687779/1‑150 (MQ=255)                                                                                                                              ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  <  2:32470/150‑1 (MQ=255)                                                                                                                              ccTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg                  <  2:1417663/150‑1 (MQ=255)                                                                                                                               ctctTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCgg                 >  2:294708/1‑150 (MQ=255)                                                                                                                               ctctTCACGGGTAAACCC‑CCCGTA‑GGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                                  >  2:1069864/1‑115 (MQ=38)                                                                                                                               ctctTCACGGGTAAACCC‑CCCGTA‑GGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCt                                                  <  1:1069864/115‑1 (MQ=38)                                                                                                                                tctTCGCGAGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                >  1:711328/1‑150 (MQ=255)                                                                                                                                  ttCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCa              >  2:312504/1‑150 (MQ=255)                                                                                                                                         gTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGaa       >  1:1455062/1‑150 (MQ=38)                                                                                                                                          tAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAg      >  1:643807/1‑150 (MQ=255)                                                                                                                                              cccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc  <  2:8883/150‑1 (MQ=255)                                                                                                                                              cccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAAC‑ACCGCCTCCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc  <  1:100273/150‑1 (MQ=255)                                                                                                                                                 gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGc                                                                                    <  2:562563/65‑1 (MQ=255)                                                                                                                                                 gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGc                                                                                    >  1:562563/1‑65 (MQ=255)                                                                                                                                                 gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                >  1:740455/1‑133 (MQ=255)                                                                                                                                                 gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                <  2:740455/133‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCG‑AACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGC  >  NC_002947/5229095‑5229385

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