Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	892,893	G→C	6.7%	S274R (AGC→AGG)	pta ←	phosphate acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	892,893	0	G	C	6.7%	247.4 / 6.9	75	S274R (AGC→AGG)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base C (2/3);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAG  >  NC_002947/892749‑893036                                                                                                                                                 |                                                                                                                                                ggTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGa                                                                                                                                             >  4:119411/1‑149 (MQ=255)        tACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAg                                                                                                                                      >  4:177563/1‑149 (MQ=255)          cGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGATGCCATTGATCGCTGCCAGGCTGACAGCGAGga                                                                                                                                    <  1:36748/149‑1 (MQ=255)          cGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGga                                                                                                                                    >  7:255187/1‑149 (MQ=255)           gAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgat                                                                                                                                   >  2:314933/1‑149 (MQ=255)             gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGatcaatg                                                                                                                                 <  8:296538/149‑7 (MQ=255)             gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgatga                                                                                                                                 >  2:231515/1‑149 (MQ=255)                ggTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATgt                                                                                                                              <  8:255187/149‑1 (MQ=255)                       aCCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTcgcg                                                                                                                       >  6:307146/1‑149 (MQ=255)                         cGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGat                                                                                                                                      >  5:143169/1‑131 (MQ=255)                         cGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGat                                                                                                                                      <  6:143169/132‑2 (MQ=255)                          gAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGa                                                                                                                                       <  8:28042/130‑1 (MQ=255)                           aaaGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCg                                                                                                                   >  5:34123/1‑149 (MQ=255)                            aaGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGc                                                                                                                  >  3:105946/1‑149 (MQ=255)                              ggATCGGCAGGCCGGCCTGCATGGCACCGCGGCACAGGCCAAGGATCCGCACGTCAGGCTTGCTGTAGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGccc                                                                                                                >  5:105355/1‑149 (MQ=255)                              ggATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGccc                                                                                                                >  1:48335/1‑149 (MQ=255)                                aTCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCgg                                                                                                              >  5:258316/1‑149 (MQ=255)                                 tCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                             <  1:231515/149‑1 (MQ=255)                                      aGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTac                                                                                                        <  3:177563/149‑1 (MQ=255)                                       ggccggccTTCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGTCTCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTacc                                                                                                       <  7:205043/149‑1 (MQ=255)                                       ggccggccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTacc                                                                                                       >  2:101400/1‑149 (MQ=255)                                             ccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAgg                                                                                                 <  1:101400/149‑1 (MQ=255)                                             ccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAgg                                                                                                 <  6:34123/149‑1 (MQ=255)                                                       aCCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAgg                                                                                       >  6:174683/1‑149 (MQ=255)                                                        ccGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                             >  5:43399/1‑126 (MQ=255)                                                        ccGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                             <  6:43399/126‑1 (MQ=255)                                                          gcgGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                    <  3:273094/149‑1 (MQ=255)                                                             gCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTa                                                                                 <  3:119411/149‑1 (MQ=255)                                                                 aGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCgg                                                                             <  5:307146/149‑1 (MQ=255)                                                                         gATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCgtgt                                                                     >  2:71224/1‑149 (MQ=255)                                                                           tACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                   <  1:314933/149‑1 (MQ=255)                                                                           tACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                   <  2:48335/149‑1 (MQ=255)                                                                              gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                    <  7:299936/129‑1 (MQ=255)                                                                              gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                    >  8:299936/1‑129 (MQ=255)                                                                              gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGc                                                                 >  2:186301/1‑148 (MQ=255)                                                                                aCGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCAc                                                              >  5:292034/1‑149 (MQ=255)                                                                                aCGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCAc                                                              >  5:90299/1‑149 (MQ=255)                                                                                     aGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtgc                                                         <  8:264808/149‑1 (MQ=255)                                                                                         ttGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtgcgtg                                                      >  1:24580/1‑148 (MQ=255)                                                                                           gCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGca                                                    <  6:258316/148‑1 (MQ=255)                                                                                              gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                              >  4:193023/1‑119 (MQ=255)                                                                                              gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                              <  3:193023/119‑1 (MQ=255)                                                                                              gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAga                                                <  4:105946/149‑1 (MQ=255)                                                                                                 gCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatga                                             >  5:128385/1‑149 (MQ=255)                                                                                                  cTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTgg                                                                  <  1:146092/127‑1 (MQ=255)                                                                                                  cTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTgg                                                                  >  2:146092/1‑127 (MQ=255)                                                                                                    ggTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatga                                          >  6:320454/1‑149 (MQ=255)                                                                                                     gTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatgat                                         >  2:12062/1‑149 (MQ=255)                                                                                                      tCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATc                                        <  6:292034/149‑1 (MQ=255)                                                                                                         gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCtt                                      >  6:2527/1‑148 (MQ=255)                                                                                                         gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCtt                                      <  5:2527/148‑1 (MQ=255)                                                                                                         gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTg                                     <  1:186301/149‑1 (MQ=255)                                                                                                             caAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                   <  4:113280/115‑1 (MQ=255)                                                                                                             caAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                   >  3:113280/1‑115 (MQ=255)                                                                                                                      cAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCgct                        >  7:188637/1‑149 (MQ=255)                                                                                                                          gggACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                    <  6:90299/149‑1 (MQ=255)                                                                                                                          gggACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                    <  8:188637/149‑1 (MQ=255)                                                                                                                                 cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCTCCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCg                   <  5:171601/143‑1 (MQ=255)                                                                                                                                 cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgt                                                           >  7:144297/1‑103 (MQ=255)                                                                                                                                 cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgt                                                           <  8:144297/103‑1 (MQ=255)                                                                                                                                 cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCg                   >  6:171601/1‑143 (MQ=255)                                                                                                                                 cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCg             >  6:102104/1‑149 (MQ=255)                                                                                                                                  aTTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGc            >  8:200773/1‑149 (MQ=255)                                                                                                                                   ttGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCCCAGATGATGATCTTGCTCATGCGCCGCTGCTCGTACTCGcc           >  6:305346/1‑149 (MQ=255)                                                                                                                                   ttGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGcc           <  6:128385/149‑1 (MQ=255)                                                                                                                                     gATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGc                                                                               >  7:157670/1‑79 (MQ=255)                                                                                                                                     gATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGc                                                                               <  8:157670/79‑1 (MQ=255)                                                                                                                                       tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                    >  5:63978/1‑136 (MQ=255)                                                                                                                                       tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                    <  6:63978/136‑1 (MQ=255)                                                                                                                                       tCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCa       >  7:123876/1‑149 (MQ=255)                                                                                                                                            gCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAg  >  4:173139/1‑149 (MQ=255)                                                                                                                                            gCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAg  <  2:24580/149‑1 (MQ=255)                                                                                                                                             ccAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                              <  4:28536/72‑1 (MQ=255)                                                                                                                                             ccAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                              >  3:28536/1‑72 (MQ=255)                                                                                                                                             ccAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAg  <  5:174683/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                GGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAG  >  NC_002947/892749‑893036

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