Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,043,929	G→A	5.9%	S13N (AGC→AAC)	PP_3565 →	LysE family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,043,929	0	G	A	5.9%	192.2 / 2.4	68	S13N (AGC→AAC)	PP_3565	LysE family transporter
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

CATGAAACCACGTCCAGAAAAGACCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTCGGCCTCAG  >  NC_002947/4043785‑4044076                                                                                                                                                 |                                                                                                                                                    cATGAAACCACGTCCAGAAAAGACCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCg                                                                                                                                                >  1:422374/1‑150 (MQ=255)   tGAAACCACGTCCAGAAAAGACCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGc                                                                                                                                              >  2:25836/1‑150 (MQ=255)                aGAAAAGACCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTAAGTGCGTGCCAGCATCGGCCTGTCGCTGAccc                                                                                                                                 >  1:1171469/1‑150 (MQ=255)                   gaactcccTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCgg                                                                                                                              <  1:492382/145‑1 (MQ=255)                    aaGACCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCggg                                                                                                                             <  1:227157/150‑1 (MQ=255)                      caCCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGc                                                                                                                           <  1:645095/149‑1 (MQ=255)                       aCCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGcc                                                                                                                          >  2:1840170/1‑150 (MQ=255)                        cccTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGcct                                                                                                                         <  2:682075/150‑2 (MQ=255)                         ccTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCGAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGACAGCATCGGCCTGTCGCTAACGTCggggggtg                                                                                                                        >  2:492382/1‑146 (MQ=255)                          cTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCaa                                                                                                                       <  2:1762357/150‑1 (MQ=255)                           ttttATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCt                                                                                                                                      <  1:1122809/134‑1 (MQ=255)                           ttttATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCt                                                                                                                                      >  2:1122809/1‑134 (MQ=255)                                   ccAGATAATACTCTTTCAGCGTGGGCGTTGCCTCCTCTACTGTCGCCTTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGtt                                                                                                              <  2:2134306/150‑1 (MQ=255)                                     aGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACGGTCGCCCCCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTTGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGACCCAACACCCTGTTgg                                                                                                            >  1:1066699/1‑150 (MQ=255)                                         aataCTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGt                                                                                                         <  1:1189958/149‑1 (MQ=255)                                         aataCTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTAACCCCGGGGCCCAACAGCCTGTTGGCGt                                                                                                         >  2:1189958/1‑149 (MQ=255)                                            aCTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGAc                                                                                                     >  1:1203142/1‑150 (MQ=255)                                            aCTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGAc                                                                                                     >  2:1446405/1‑150 (MQ=255)                                              tctTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAACATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGAccc                                                                                                   >  1:73868/1‑150 (MQ=255)                                              tctTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAACATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGAccc                                                                                                   >  1:73776/1‑150 (MQ=255)                                                   cAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAACGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGgc                                                                                              <  2:1686136/150‑1 (MQ=255)                                                        tGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGct                                                                                         <  2:771271/150‑1 (MQ=255)                                                              ttGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATgg                                                                                   <  2:2244567/150‑1 (MQ=255)                                                                 cgcCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATg                                                                                    <  1:1118026/146‑1 (MQ=255)                                                                 cgcCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATg                                                                                    >  2:1118026/1‑146 (MQ=255)                                                                    cATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACg                                                                             <  2:1649478/150‑1 (MQ=255)                                                                     aTAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGc                                                                                                                   <  1:1487476/111‑1 (MQ=255)                                                                     aTAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGc                                                                                                                   >  2:1487476/1‑111 (MQ=255)                                                                             tCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCt                                                                    <  1:25836/150‑1 (MQ=255)                                                                               gCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGt                                                                  >  1:481507/1‑150 (MQ=255)                                                                                cctgctcgcGACCAAGGCAGTTTTCAACATGGTACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACCGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGtt                                                                 <  1:443917/146‑1 (MQ=255)                                                                                  ctctctCGACCAAGGCAGTTTTCGACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCa                                                               <  1:968862/150‑1 (MQ=255)                                                                                     tctcGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCGTTTCGCTGACCCAGGTGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCa                                                            >  1:755838/1‑150 (MQ=255)                                                                                     tctcGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGgc                                                                                              <  2:874074/116‑1 (MQ=255)                                                                                     tctcGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGgc                                                                                              >  1:874074/1‑116 (MQ=255)                                                                                      ctcGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATTGCGCGCTCTATGGCGCACGGCGAAc                                                                       >  2:271358/1‑138 (MQ=255)                                                                                      ctcGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAAc                                                                       <  1:271358/138‑1 (MQ=255)                                                                                          aCCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGGCCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACGATCGTg                                                       >  1:377148/1‑150 (MQ=255)                                                                                            cAAGGCAGTGTTTTACATGGCACTCGCCACTTGGCTCATCTACCTGCTGGCGAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTggg                                                     >  2:1256571/1‑150 (MQ=255)                                                                                              aGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAACATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCg                                                   <  2:73868/150‑1 (MQ=255)                                                                                              aGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAACATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCg                                                   <  2:73776/150‑1 (MQ=255)                                                                                               gggaGTTTTCAACCTGGCTCTCTTCGCCTGGCTTATCTACCTGTTGGCCGGCATCGGCCTGTCGCTGTCCCGGGGGCCCAACAGCCCGTTGGGGTGGACCCATGGAGCGCGCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCgg                                                  <  2:377148/147‑1 (MQ=255)                                                                                                gCAGTTTTCAACATGGCACTGTACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCAAAATCGTGGGCggc                                                 <  2:1203142/150‑2 (MQ=255)                                                                                                gCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCggg                                                 <  1:1840170/150‑1 (MQ=255)                                                                                                gCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCg                                                   >  2:160184/1‑148 (MQ=255)                                                                                                gCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCg                                                   <  1:160184/148‑1 (MQ=255)                                                                                                     tttCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCtt                                            >  1:2070050/1‑150 (MQ=255)                                                                                                       tCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTgg                                          <  2:422374/150‑1 (MQ=255)                                                                                                         aaCATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCGGGGGt                                        >  1:821307/1‑150 (MQ=255)                                                                                                                 aCTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCggg                                                 >  1:1241297/1‑133 (MQ=255)                                                                                                                 aCTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCggg                                                 <  2:1241297/133‑1 (MQ=255)                                                                                                                   tGGACACCTGGCGCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGccc                              >  1:1827068/1‑150 (MQ=255)                                                                                                                    ggACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATCGCGCCCGGCGTGCGCGGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCt                             <  2:1834529/150‑1 (MQ=255)                                                                                                                          cTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCa                                                               <  2:267247/110‑1 (MQ=255)                                                                                                                          cTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCa                                                               >  1:267247/1‑110 (MQ=255)                                                                                                                          cTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTggg                                         >  2:595704/1‑132 (MQ=255)                                                                                                                          cTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTggg                                         <  1:595704/132‑1 (MQ=255)                                                                                                                               tCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGa                  >  1:1348421/1‑150 (MQ=255)                                                                                                                                cATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTTGGTTTCAGCGCCCTCATCGCCTTGAc                 <  2:481507/150‑1 (MQ=255)                                                                                                                                 aTCTACCTGCTGACCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCCGCGCCCTCATCGCCTTGAcc                >  2:155491/1‑150 (MQ=255)                                                                                                                                      ccTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGtt           >  2:2238501/1‑150 (MQ=255)                                                                                                                                       ctgctgGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTc          <  2:1348421/150‑1 (MQ=255)                                                                                                                                        tgctgGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTCg         <  2:2070050/150‑1 (MQ=255)                                                                                                                                        tgctgGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTCg         >  1:246447/1‑150 (MQ=255)                                                                                                                                        tgctgGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTCg         >  2:1587904/1‑150 (MQ=255)                                                                                                                                             gCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGtg                                              >  2:1776376/1‑107 (MQ=255)                                                                                                                                             gCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTc                                              <  1:1776376/108‑1 (MQ=255)                                                                                                                                               cAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTCGGCCTCAg  >  1:953541/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CATGAAACCACGTCCAGAAAAGACCCTTTTATGACCAGATAATACTCTTTCAGCGTGAGCGTTGCGCCATAGACTGTCGCCCTCTCTCGACCAAGGCAGTTTTCAACATGGCACTGGACACCTGGCTCATCTACCTGCTGGCCAGCATCGGCCTGTCGCTGACCCCGGGGCCCAACAGCCTGTTGGCGTTGACCCATGGCGCGCTCTATGGCGCACGGCGAACGCTGTTCACCATCGTGGGCGGGGTCTTGGGTTTCAGCGCCCTCATCGCCTTGACCATGTTCGGCCTCAG  >  NC_002947/4043785‑4044076

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