Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,606,687	G→A	7.3%	D501D (GAC→GAT)	cyaB ←	cyclolysin secretion/processing ATP‑binding protein CyaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,606,687	0	G	A	7.3%	178.0 / 3.1	65	D501D (GAC→GAT)	cyaB	cyclolysin secretion/processing ATP‑binding protein CyaB
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCAATACCCGCCCCCGCTCAGGTGAGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATACCCCTGTCTGCTGGAAGCTG  >  NC_002947/5606545‑5606834                                                                                                                                               |                                                                                                                                                     gCAATACCCGCCCCCGCTCAGGTGAGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGcc                                                                                                                                               >  2:925643/1‑150 (MQ=255) gCAATACCCGCCCCCGCTCAGGTGAGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGcc                                                                                                                                               >  1:1861127/1‑150 (MQ=255)   aaTACCCGCCCCCGCTCAGGTGAGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGa                                                                                                                                             >  2:1883111/1‑150 (MQ=255)                 cTCAGGTGAGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAacacac                                                                                                                               <  2:268477/150‑1 (MQ=255)                         aGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGa                                                                                                                       <  1:348186/150‑1 (MQ=255)                              aGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCg                                                                                                                 >  1:1795733/1‑150 (MQ=255)                                ccGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGAGCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGGTGACATAGAGACCCCCCTCCCCGCCACCCGGCCGATAGCGAAAACACACACTATCGAGC‑TCGAt                                                                                                               <  2:19813/150‑1 (MQ=255)                                ccGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGAt                                                                                                               >  2:1139463/1‑150 (MQ=255)                                ccGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGAt                                                                                                               <  2:3639/150‑1 (MQ=255)                                 cGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATc                                                                                                              >  2:526378/1‑150 (MQ=255)                                    tGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGc                                                                                                           >  1:217802/1‑150 (MQ=255)                                     gAATCAAGCGTGTCAGCGTGCTCTTGCCTGTCCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGcc                                                                                                          <  1:1827265/150‑1 (MQ=255)                                          aaGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACa                                                                                                     >  2:1460056/1‑150 (MQ=255)                                                 tCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGccc                                                                                                         >  1:1500901/1‑139 (MQ=255)                                                 tCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGccc                                                                                                         <  2:1500901/139‑1 (MQ=255)                                                                     cccGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAgt                                                                                                  <  2:1298729/126‑1 (MQ=255)                                                                     cccGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAgt                                                                                                  >  1:1298729/1‑126 (MQ=255)                                                                     cccGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCccttgcct                                                                          >  2:1290529/1‑150 (MQ=255)                                                                     cccGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCATCCGGCCGATAGCGAAAaca                                                                                                                                  >  1:1637543/1‑95 (MQ=255)                                                                     cccGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCATCCGGCCGATAGCGAAAaca                                                                                                                                  <  2:1637543/95‑1 (MQ=255)                                                                       cGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCccttgccttg                                                                        <  1:1428551/150‑1 (MQ=255)                                                                        gAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGATC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTgg                                                                       <  2:1838431/150‑1 (MQ=255)                                                                           cGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCTACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTg                                                                    >  1:1557843/1‑150 (MQ=255)                                                                              ccTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGAcc                                                                 >  2:440781/1‑150 (MQ=255)                                                                                 aCAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTc                                                               <  2:424483/149‑1 (MQ=255)                                                                                 aCAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTc                                                               >  1:424483/1‑149 (MQ=255)                                                                                 aCAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCg                                                              >  2:596167/1‑150 (MQ=255)                                                                                  cAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCgg                                                             >  2:896667/1‑150 (MQ=255)                                                                                   aaCGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGAc                                                                  >  2:1793203/1‑144 (MQ=255)                                                                                   aaCGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGAc                                                                  <  1:1793203/144‑1 (MQ=255)                                                                                         tATCACTTCTCCGGCGGCGATCTCCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGgt                                                      <  1:395850/150‑1 (MQ=255)                                                                                              cTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTga                                                 <  1:566172/150‑1 (MQ=255)                                                                                              cTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTga                                                 >  2:1625428/1‑150 (MQ=255)                                                                                                     ggcggcGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTc                                          <  1:1663224/150‑1 (MQ=255)                                                                                                          agATCTTCAGGCTGACAGCCGGCACCACTTCGCCCCCGTCCGGCCGATAGCGTAGCCA‑GTACGCTCTTGCTTACAGCTACCCGACAAGGGGTGGCACCACCCTGCCCCTCAATTGGTTGGCCTCGTGACGGGTGTTGAGAATGTCGCCCa                                     <  1:925643/149‑1 (MQ=255)                                                                                                           gATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAg                                    >  2:1621096/1‑150 (MQ=255)                                                                                                             tCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑ACGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGc                                  <  1:526378/150‑1 (MQ=255)                                                                                                              cTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCg                                 >  2:123125/1‑150 (MQ=255)                                                                                                               ttCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGt                                <  2:1557843/150‑1 (MQ=255)                                                                                                               ttCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGt                                >  1:1629244/1‑150 (MQ=255)                                                                                                               ttCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGt                                <  2:1861127/150‑1 (MQ=255)                                                                                                                tCAGGCTGACATCGCGCACCACTTCGCCCCCGCCAGGCCCACACCGACCACACCCACGACTGAAT‑GCTGGCTTCCGGTGAAGGGTTGGTGGCACGGTTCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGtt                               <  2:217802/150‑1 (MQ=255)                                                                                                                  aGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGa                             >  1:1431806/1‑150 (MQ=255)                                                                                                                    gCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTGGGTACGGGTGTGGAGAATGTCGCCCAGGTGTTGAAt                           >  2:556732/1‑150 (MQ=255)                                                                                                                     cTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATc                          <  1:1621096/150‑1 (MQ=255)                                                                                                                       gACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCATGCGTTGAATCGa                        <  2:1629244/150‑1 (MQ=255)                                                                                                                        aCATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTggg                                                                                             <  1:1056273/80‑1 (MQ=255)                                                                                                                        aCATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTggg                                                                                             >  2:1056273/1‑80 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGtt                                                                     >  2:1353363/1‑102 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGtt                                                                     <  1:1353363/102‑1 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATg                                            >  1:1446689/1‑127 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATg                                            <  2:1446689/127‑1 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATAc                     >  2:753857/1‑150 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCATCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACgg                                                        >  2:1083186/1‑115 (MQ=255)                                                                                                                          aTCGCGCACCACTTCGCCCCCATCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACgg                                                        <  1:1083186/115‑1 (MQ=255)                                                                                                                            cgcgCACCACTTCGCCCCCGTCCGGCAGCTTGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATAccc                   <  1:1460056/150‑1 (MQ=255)                                                                                                                               gcACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAgtgt                                                                                                >  1:287619/1‑70 (MQ=255)                                                                                                                               gcACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAgtgt                                                                                                <  2:287619/70‑1 (MQ=255)                                                                                                                               gcACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATACCCCTg                >  2:1167992/1‑150 (MQ=255)                                                                                                                                caccacTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATACCCCTGt               <  1:1139463/150‑1 (MQ=255)                                                                                                                                   cacTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTg                                                  <  1:558114/112‑1 (MQ=255)                                                                                                                                   cacTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTg                                                  >  2:558114/1‑112 (MQ=255)                                                                                                                                     cTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATACCCCTGTctgct          <  2:1795733/150‑1 (MQ=255)                                                                                                                                     cTTCGCCCCCGTCCGGCCGATAACGAAAACACACACGATCGAAC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATTTCGCCTAGGCGTTGAATCGATACCCCTGTctgct          <  2:256235/150‑1 (MQ=255)                                                                                                                                             ccGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATACCCCTGTCTGCTGGAAGCTg  >  2:1279748/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                     GCAATACCCGCCCCCGCTCAGGTGAGTACAGCCGCTGAATCAAGCGTGTCAGCGTGCTCTTGCCTGACCCCGAGCGACCTACAACGCCTATCACTTCTCCGGCGGCGATCTTCAGGCTGACATCGCGCACCACTTCGCCCCCGTCCGGCCGATAGCGAAAACACACACGATCGAGC‑TCGATCTGCCCGACAAGTGTGGGCAGCACCGTGCCCTTGCCTTGGTTGACCTCGGTACGGGTGTTGAGAATGTCGCCCAGGCGTTGAATCGATACCCCTGTCTGCTGGAAGCTG  >  NC_002947/5606545‑5606834

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