Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,918,174	C→G	7.7%	L191L (CTC→CTG)	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,918,174	0	C	G	7.7%	203.4 / 4.3	64	L191L (CTC→CTG)	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base G (2/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAA  >  NC_002947/4918037‑4918318                                                                                                                                          |                                                                                                                                                 tACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCggcaggca                                                                                                                                       <  1:203837/149‑1 (MQ=255)   ccGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTGAGCggcacgaag                                                                                                                                      <  7:241483/148‑6 (MQ=255)   ccGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGgcaggcag                                                                                                                                      <  8:77888/148‑1 (MQ=255)           tGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACt                                                                                                                              >  1:127067/1‑148 (MQ=255)             gTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAg                                                                                                                           >  1:94787/1‑149 (MQ=255)              tCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGCGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACt                                                                                                                              <  5:194572/145‑1 (MQ=255)              tCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGAGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACt                                                                                                                              >  6:194572/1‑145 (MQ=255)                aGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCAt                                                                                                                        >  5:245129/1‑149 (MQ=255)                 gCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCAtt                                                                                                                       >  8:144630/1‑149 (MQ=255)                  ccTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTg                                                                                                                      <  4:117280/149‑1 (MQ=255)                   cTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGc                                                                                                                          >  1:253055/1‑144 (MQ=255)                   cTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGc                                                                                                                          <  2:253055/144‑1 (MQ=255)                     gCTCAACAGCAACCAGGCCGGGGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGccc                                                                                                                   <  1:199565/149‑1 (MQ=255)                       tCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCaa                                                                                                                 >  5:45512/1‑149 (MQ=255)                            aGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCggca                                                                                                                                           <  4:218740/118‑1 (MQ=255)                            aGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgct                                                                                                            <  6:107314/149‑1 (MQ=255)                             gCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgctg                                                                                                           >  5:171669/1‑149 (MQ=255)                              cAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGTCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCTTTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCATCAgctgc                                                                                                          <  5:132637/149‑1 (MQ=255)                                 ccAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGa                                                                                                       <  7:149516/149‑1 (MQ=255)                                  cAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGATTCGTGCAGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGAc                                                                                                      >  4:170099/1‑149 (MQ=255)                                     gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCa                                                                                                   <  8:163006/149‑1 (MQ=255)                                      ccGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAg                                                                                                  <  3:170099/149‑1 (MQ=255)                                                acTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCa                                                                                                                            >  6:133431/1‑113 (MQ=255)                                                acTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCa                                                                                                                            <  5:133431/113‑1 (MQ=255)                                                acTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGc                                                                                         >  6:24133/1‑148 (MQ=255)                                                  cgAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGcc                                                                                        >  4:197950/2‑147 (MQ=255)                                                    aGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCaa                                                                                      <  7:144630/147‑1 (MQ=255)                                                         cGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGc                                                                               >  3:206631/1‑149 (MQ=255)                                                           aCAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTa                                                                             <  6:45512/149‑1 (MQ=255)                                                                 cgcgcccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGc                                                                       <  2:94787/149‑1 (MQ=255)                                                                         cTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCaa                                                               <  6:93565/149‑1 (MQ=255)                                                                            gTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGc                                                            >  3:154516/1‑149 (MQ=255)                                                                            gTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGc                                                            >  3:6071/1‑149 (MQ=255)                                                                               aTTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGTCAAGCCTCGCTACCCGGCCTGCTCAACGCCCt                                                         >  2:73434/1‑149 (MQ=255)                                                                               aTTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGccc                                                          <  4:206631/148‑1 (MQ=255)                                                                                      gCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAg                                                  <  8:20658/149‑1 (MQ=255)                                                                                            tCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCaaaa                                                                                             <  1:136677/100‑1 (MQ=255)                                                                                            tCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCaaaa                                                                                             >  2:136677/1‑100 (MQ=255)                                                                                            tCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACa                                                   >  7:238432/1‑142 (MQ=255)                                                                                            tCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACa                                                   <  8:238432/142‑1 (MQ=255)                                                                                              aGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCcg                                          >  8:43203/1‑149 (MQ=255)                                                                                              aGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCcg                                          <  6:171669/149‑1 (MQ=255)                                                                                                   ccTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGAcc                                     >  1:112562/1‑149 (MQ=255)                                                                                                    cTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGAccc                                    <  3:197950/149‑1 (MQ=255)                                                                                                       gTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGTCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAg                                  <  1:73434/148‑1 (MQ=255)                                                                                                            gggTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTc                             >  7:244529/1‑148 (MQ=255)                                                                                                               tGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATg                         <  2:112562/149‑1 (MQ=255)                                                                                                                  tCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATgg                        <  7:134064/147‑1 (MQ=255)                                                                                                                  tCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATgg                        >  8:134064/1‑147 (MQ=255)                                                                                                                     gTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAcc                   >  3:257157/1‑149 (MQ=255)                                                                                                                        aTTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgc                >  8:29521/1‑149 (MQ=255)                                                                                                                            cgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgcg             <  4:6071/148‑1 (MQ=255)                                                                                                                              ctgcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACACAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGTAGTCGATGGCGACCTgc                <  5:208318/138‑1 (MQ=255)                                                                                                                              ctgcctGCCGCTGAGCGGAAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGCCCCAGCAAAACCGCCACGCCTCGCTACCCGGCCTGCTCACCGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgagctc          >  6:208318/6‑147 (MQ=255)                                                                                                                               ttCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAg         >  3:161968/1‑149 (MQ=255)                                                                                                                                gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGc                                                                       >  3:19818/4‑86 (MQ=255)                                                                                                                                gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGc                                                                       <  4:19818/83‑1 (MQ=255)                                                                                                                                   tGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCt                  <  3:43498/136‑1 (MQ=255)                                                                                                                                   tGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCt                  >  4:43498/1‑136 (MQ=255)                                                                                                                                      cGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGaa  >  5:30758/1‑149 (MQ=255)                                                                                                                                       gCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGc                                <  7:68690/118‑1 (MQ=255)                                                                                                                                       gCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGc                                >  8:68690/1‑118 (MQ=255)                                                                                                                                       gCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGcc           <  1:267878/139‑1 (MQ=255)                                                                                                                                       gCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGcc           >  2:267878/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                                 TACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGATCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCCGAA  >  NC_002947/4918037‑4918318

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