Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,533,762	C→G	5.9%	I577M (ATC→ATG)	secA →	protein translocase subunit SecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,533,762	0	C	G	5.9%	212.1 / 4.9	68	I577M (ATC→ATG)	secA	protein translocase subunit SecA
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CTGACCATCGCCACCAACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGC  >  NC_002947/1533613‑1533911                                                                                                                                                      |                                                                                                                                                       cTGACCATCGCCACCAACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATc                                                                                                                                                        <  2:1644633/150‑1 (MQ=255)        tCGCCACCAACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACg‑g                                                                                                                                                >  1:1103123/1‑150 (MQ=255)           ccaccaACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTg                                                                                                                                             >  2:716268/1‑150 (MQ=255)            caccaACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTgg                                                                                                                                            >  2:616255/1‑150 (MQ=255)             accaACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATGGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGc                                                                                                                                           >  2:902792/1‑150 (MQ=255)               caACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCt                                                                                                                                         <  2:553904/150‑1 (MQ=255)               caACAGGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCt                                                                                                                                         <  2:170063/150‑1 (MQ=255)                 aCATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGc                                                                                                                                       <  1:1432494/150‑1 (MQ=255)                    tggctggcCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATg                                                                                                                                    <  1:1449794/150‑1 (MQ=255)                         ggcCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATc                                                                                                                               <  1:98138/150‑1 (MQ=255)                          gcCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCg                                                                                                                              <  1:1174215/150‑1 (MQ=255)                            cGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCt                                                                                                                            >  2:449007/1‑150 (MQ=255)                             gTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCtt                                                                                                                           <  1:194837/150‑1 (MQ=255)                              tGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTc                                                                                                                          <  1:334316/150‑1 (MQ=255)                              tGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTc                                                                                                                          <  1:259813/150‑1 (MQ=255)                               ggTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGGCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCCGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTcc                                                                                                                         <  1:1566997/150‑1 (MQ=255)                                  aCCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAg                                                                                                                      <  1:1045616/150‑1 (MQ=255)                                    cGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAgcg                                                                                                                    >  1:28205/1‑150 (MQ=255)                                    cGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGCAATCGAGACG‑GGTGGCATGCGTGTAATCGCTTCCGAgcg                                                                                                                    >  1:1594627/1‑150 (MQ=255)                                     gACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAgcgc                                                                                                                   <  2:28205/150‑1 (MQ=255)                                           cTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCCAACGATAGCGAGG‑GGGTGCTGGCTTGTAATCGCTTCCGAGCGCCACGaa                                                                                                             >  1:1184588/1‑150 (MQ=255)                                             gCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATc                                                                                                           <  1:1388746/150‑1 (MQ=255)                                                  gCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCGTGTAATCGCTTCCGAGCGCCACGAATCGCGcc                                                                                                      >  1:1208856/1‑150 (MQ=255)                                                      tAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTAt                                                                                                  <  1:1234090/150‑1 (MQ=255)                                                            ggAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCt                                                                                                                                         <  1:433181/105‑1 (MQ=255)                                                            ggAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCt                                                                                                                                         >  2:433181/1‑105 (MQ=255)                                                            ggAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACaa                                                                                            <  2:380943/150‑1 (MQ=255)                                                            ggAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGACCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACaa                                                                                            >  1:1462713/1‑150 (MQ=255)                                                             gAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAAc                                                                                           <  2:1462713/150‑1 (MQ=255)                                                                 ccGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGc                                                                                       <  1:449007/150‑1 (MQ=255)                                                                       tGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGcgcg                                                                                 <  1:616255/150‑1 (MQ=255)                                                                         gCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGAAC‑CGACGAATTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGc                                                                               >  1:515532/1‑150 (MQ=255)                                                                             cccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCt                                                                           >  1:889326/1‑150 (MQ=255)                                                                              ccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATGGAGACG‑GGTGGCCt                                                                                                                                         >  1:1046992/1‑87 (MQ=255)                                                                              ccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATGGAGACG‑GGTGGCCt                                                                                                                                         >  1:1046984/1‑87 (MQ=255)                                                                              ccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATGGAGACG‑GGTGGCCt                                                                                                                                         <  2:1046992/87‑1 (MQ=255)                                                                              ccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATGGAGACG‑GGTGGCCt                                                                                                                                         <  2:1046984/87‑1 (MQ=255)                                                                              ccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTc                                                                          <  1:1213737/150‑1 (MQ=255)                                                                              ccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTc                                                                          >  1:275442/1‑150 (MQ=255)                                                                               cTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTcc                                                                         >  1:1380447/1‑150 (MQ=255)                                                                                tGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCt                                                                                                                                         <  2:800509/85‑1 (MQ=255)                                                                                tGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCt                                                                                                                                         >  1:800509/1‑85 (MQ=255)                                                                                tGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCg                                                                        >  1:1038351/1‑150 (MQ=255)                                                                                       cccGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCa                                                                 >  2:380884/1‑150 (MQ=255)                                                                                        ccGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAg                                                                >  1:81495/1‑150 (MQ=255)                                                                                         cGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAgg                                                               <  1:785297/150‑1 (MQ=255)                                                                                         cGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAgg                                                               <  2:889326/150‑1 (MQ=255)                                                                                               ccGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGAcc                                                         >  2:885918/1‑150 (MQ=255)                                                                                                 gaggcgATGGCCCAGGTCAAGGCCGTTGTGGAAAAGGGTCACCAGCAGGTGGTCGTGACG‑GGTGGGCTTCATGTAGTCGCTTCCGCGGTCCACGAATCGCGCCGTATCGACCAGCACATGCCCGGCCGCTCCGGTCGGCAGGCCGGCCCg                                                       >  2:155947/6‑150 (MQ=255)                                                                                                   gCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCggg                                                     >  1:1068272/1‑150 (MQ=255)                                                                                                   gCAGATCGCCCAGAGCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCggg                                                     <  2:1431174/150‑1 (MQ=255)                                                                                                          gCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGc                                              <  2:1184588/150‑1 (MQ=255)                                                                                                            ccAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCg                                            >  1:354014/1‑150 (MQ=255)                                                                                                                  cAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTa                                      <  1:885918/150‑1 (MQ=255)                                                                                                                   aaGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTAc                                     <  2:591497/150‑1 (MQ=255)                                                                                                                                 gAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAgn                       <  1:380884/150‑2 (MQ=255)                                                                                                                                   aGCGTCACCAGCAAGTGACCGAGACG‑GGTGGCCTGCGTGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACa                     <  2:1208856/150‑1 (MQ=255)                                                                                                                                    gCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGcc                                                                                                      >  2:498506/1‑68 (MQ=255)                                                                                                                                    gCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGcc                                                                                                      <  1:498506/68‑1 (MQ=255)                                                                                                                                    gCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAg                    >  1:898104/1‑150 (MQ=255)                                                                                                                                     cGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGc                   <  2:354014/150‑1 (MQ=255)                                                                                                                                        cACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTg                >  2:995384/1‑150 (MQ=255)                                                                                                                                           cAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGgtcg                                                                    <  2:733687/95‑1 (MQ=255)                                                                                                                                           cAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGgtcg                                                                    >  1:733687/1‑95 (MQ=255)                                                                                                                                           cAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATg             >  1:750284/1‑150 (MQ=255)                                                                                                                                           cAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATg             >  1:906893/1‑150 (MQ=255)                                                                                                                                            aGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATgc            >  1:523877/1‑150 (MQ=255)                                                                                                                                                 gTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCa                                                                 <  2:1570784/92‑1 (MQ=255)                                                                                                                                                 gTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCACCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCa                                                                 >  1:1570784/1‑92 (MQ=255)                                                                                                                                                      cGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGc  <  2:750284/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CTGACCATCGCCACCAACATGGCTGGCCGTGGTACCGACATCCTGCTTGGCGGTAACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACG‑GGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGC  >  NC_002947/1533613‑1533911

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