Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,045,272	G→A	10.7%	intergenic (+42/+64)	PP_0905 → / ← PP_0906	hypothetical protein/RND family multidrug transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,045,272	0	G	A	10.7%	162.1 / 5.5	61	intergenic (+42/+64)	PP_0905/PP_0906	hypothetical protein/RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (4/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTG  >  NC_002947/1045123‑1045420                                                                                                                                                      |                                                                                                                                                     gcgTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgg                                                                                                                                                      <  2:94298/150‑1 (MQ=255)   gTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgggg                                                                                                                                                    >  1:632224/1‑150 (MQ=255)        cgcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCttt                                                                                                                                               <  1:1051229/150‑1 (MQ=255)         gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                              <  2:424617/150‑1 (MQ=255)         gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGGCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTCGCTGCGAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                              <  1:1002174/150‑1 (MQ=255)          cgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCCCGCAATCGCGAGGCCTCTgg                                                                                                                                             <  1:387914/150‑1 (MQ=255)            ggATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtt                                                                                                                                           >  2:220688/1‑150 (MQ=255)                        tcagcagcagGCCCTGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGGTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGc                                                                                                                               <  1:269015/149‑1 (MQ=255)                         cagcagcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGcc                                                                                                                              >  2:604159/1‑150 (MQ=255)                              gcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGcc                                                                                                                         <  2:528975/150‑1 (MQ=255)                                          cgagcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgcg                                                                                                             <  2:632224/150‑1 (MQ=255)                                               cAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtc                                                                                                                                                <  1:1032273/110‑2 (MQ=39)                                               cAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtc                                                                                                                                                >  2:1032273/1‑109 (MQ=38)                                                 gACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGc                                                                                                                                                   >  2:287148/1‑105 (MQ=255)                                                 gACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGc                                                                                                                                                   <  1:287148/105‑1 (MQ=255)                                                     aGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAGCCCTGCGTTAGCTGGAAAAGACAGAAGCCCCGCCACTGCGAGCCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTAGCGGGCAgcccgcc                                                                                                  >  1:875650/1‑150 (MQ=255)                                                      ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgccc                                                                                                 >  1:1149794/1‑150 (MQ=255)                                                      ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgc                                                                                                       >  2:511842/1‑144 (MQ=255)                                                      ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgc                                                                                                       <  1:511842/144‑1 (MQ=255)                                                          cTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAg                                                                                             <  1:490503/150‑1 (MQ=255)                                                               gCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttttt                                                                                                                                        <  2:462037/102‑1 (MQ=255)                                                               gCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttttt                                                                                                                                        >  1:462037/1‑102 (MQ=255)                                                                cTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGt                                                                                       >  2:761253/1‑150 (MQ=255)                                                                 tGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGtt                                                                                      <  1:740934/150‑1 (MQ=255)                                                                 tGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGtt                                                                                      <  2:562906/150‑1 (MQ=255)                                                                  gTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCctct                                                                                                                  >  1:135708/1‑121 (MQ=255)                                                                  gTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCctct                                                                                                                  <  2:135708/121‑1 (MQ=255)                                                                  gTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTc                                                                                     >  2:1165238/1‑150 (MQ=255)                                                                   ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGtt                                                                                      >  1:1207192/1‑148 (MQ=255)                                                                   ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGtt                                                                                      <  2:1207192/148‑1 (MQ=255)                                                                      cTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGc                                                                                 >  1:475342/1‑150 (MQ=255)                                                                               cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAgg                                                                        >  1:371333/1‑150 (MQ=255)                                                                               cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAACCAAAGGCCCCGCAATCGCGAGGCCTCTggc                                                                                                                                            >  2:1183132/1‑81 (MQ=17)                                                                               cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAACCAAAGGCCCCGCAATCGCGAGGCCTCTggc                                                                                                                                            <  1:1183132/82‑2 (MQ=17)                                                                               cTGCGGCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAgg                                                                        >  2:239868/1‑150 (MQ=255)                                                                                tGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCTGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAgga                                                                       <  1:220688/150‑2 (MQ=255)                                                                                            tGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTc                                                           <  1:239868/150‑1 (MQ=255)                                                                                                 tCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCtt                                                      <  2:475342/150‑1 (MQ=255)                                                                                                   acaGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTGGCCGTCCTCTTGGCGGGCAGCCCGCCCCAAGTTCGGTTCAAGCGATTGCAGGTGCCGTGGCGCTCGCCTtgt                                                    >  2:1156934/1‑150 (MQ=255)                                                                                                    caGGGATGAACCCTGCGTTAGCTGGGAAAGACAGAAGCCTCGCGAGTGCGGGGCCGTTGGTTTTTTAAAGCTGTCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgtg                                                   >  1:687912/1‑150 (MQ=255)                                                                                                    caGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCggg                                                                                                           <  2:1127955/94‑1 (MQ=255)                                                                                                    caGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCggg                                                                                                           >  1:1127955/1‑94 (MQ=255)                                                                                                    caGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgtg                                                   >  2:122252/1‑150 (MQ=255)                                                                                                     aGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGCCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgc                                                  <  2:875650/150‑1 (MQ=255)                                                                                                       ggATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgcgc                                                <  1:1165238/150‑1 (MQ=255)                                                                                                         aTGAACCCTGCGTTAGCTGGAAAATCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGccc                                              >  1:1198703/1‑149 (MQ=255)                                                                                                         aTGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgtg                                                   >  2:627413/1‑145 (MQ=255)                                                                                                         aTGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgtg                                                   <  1:627413/145‑1 (MQ=255)                                                                                                               ccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAggg                                        <  2:1149794/150‑1 (MQ=255)                                                                                                                cTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                                       <  1:122252/150‑1 (MQ=255)                                                                                                                 tGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgtg                                      >  1:304046/1‑150 (MQ=255)                                                                                                                  gCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgtgt                                     <  2:304046/150‑1 (MQ=255)                                                                                                                     ttAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGa                                  <  2:1198703/150‑1 (MQ=255)                                                                                                                          tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgc                                                                                                              <  1:1078942/69‑1 (MQ=255)                                                                                                                          tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgc                                                                                                              >  2:1078942/1‑69 (MQ=255)                                                                                                                                   caaagGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAAcaca                    >  1:1081747/4‑150 (MQ=255)                                                                                                                                   cAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCtt                                                                  <  1:561720/104‑1 (MQ=255)                                                                                                                                   cAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCtt                                                                  >  2:561720/1‑104 (MQ=255)                                                                                                                                    agagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaa                   >  1:761149/1‑150 (MQ=255)                                                                                                                                        gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgt                                                    <  2:260576/113‑1 (MQ=255)                                                                                                                                        gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgt                                                    >  1:260576/1‑113 (MQ=255)                                                                                                                                               gATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAaac        >  1:550761/1‑150 (MQ=255)                                                                                                                                                    cGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGt   >  1:72298/1‑150 (MQ=255)                                                                                                                                                     gggtCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTg  <  1:604159/150‑1 (MQ=255)                                                                                                                                                     ggggCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTg  >  2:891995/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTG  >  NC_002947/1045123‑1045420

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