Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,981,093	C→A	7.9%	intergenic (+32/‑268)	rpsA → / → ihfB	30S ribosomal protein S1/integration host factor subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,981,093	0	C	A	7.9%	197.0 / 5.2	63	intergenic (+32/‑268)	rpsA/ihfB	30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base A (2/3);  total (31/32)
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.76e-01

TCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAC  >  NC_002947/1980945‑1981235                                                                                                                                                     |                                                                                                                                                 tCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGTACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTc                                                                                                                                                  <  1:245431/149‑1 (MQ=255)          gacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGccc                                                                                                                                         <  1:357574/149‑1 (MQ=255)           acgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCt                                                                                                                                        <  5:149299/149‑1 (MQ=255)                 aaGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCtttttac                                                                                                                                  <  5:87788/149‑3 (MQ=255)                   gagCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCGCCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTg                                                                                                                                >  1:125832/1‑149 (MQ=255)                      cgcgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTc                                                                                                                             >  5:234056/1‑149 (MQ=255)                        cgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcg                                                                                                                           >  6:163744/1‑149 (MQ=255)                         gAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCttttttt                                                                                                                                  >  5:134639/1‑141 (MQ=255)                         gAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCttttttt                                                                                                                                  <  6:134639/141‑1 (MQ=255)                             cTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCGCCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGAcc                                                                                                                      <  2:125832/149‑1 (MQ=255)                                tCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCTTTTTTTTGCTCCGCGACCAgg                                                                                                                   <  8:216747/149‑1 (MQ=255)                                  cAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAggtg                                                                                                                 >  3:37068/1‑149 (MQ=255)                                         tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                          >  5:49091/1‑149 (MQ=255)                                         tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                          >  1:338596/1‑149 (MQ=255)                                          gAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggcg                                                                                                         >  4:345550/1‑149 (MQ=255)                                           aaaaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggcg                                                                                                         >  6:259927/1‑148 (MQ=255)                                             aaaCGCTCCGGAAGCTGCTTCCGACACCACCATGTCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGAt                                                                                                      <  3:240405/149‑1 (MQ=255)                                              accGCTCCGGAAGCGGCTGCCGACACCACCATGGCTGCGCTGCGGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATg                                                                                                     <  7:124521/147‑1 (MQ=255)                                               ccGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATg‑‑g                                                                                                  <  7:273750/148‑1 (MQ=255)                                                 gCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCa                                                                                                <  5:163744/149‑1 (MQ=255)                                                 gCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCa                                                                                                >  8:280657/1‑149 (MQ=255)                                                   tCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGAt                                                                                                      <  5:232307/143‑1 (MQ=255)                                                   tCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGAt                                                                                                      >  6:232307/1‑143 (MQ=255)                                                     cGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCActgc                                                                                            >  3:93979/1‑149 (MQ=255)                                                       gAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCActgctg                                                                                          <  5:313984/149‑1 (MQ=255)                                                               tgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGCAAACGc                                                                                  >  7:129727/1‑149 (MQ=255)                                                               tgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGc                                                                                  >  7:134648/1‑149 (MQ=255)                                                                   gACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGc                                                                              >  5:118121/1‑149 (MQ=255)                                                                      accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAggtg                                                                                                                 >  1:92703/1‑113 (MQ=255)                                                                      accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAggtg                                                                                                                 <  2:92703/113‑1 (MQ=255)                                                                       ccaccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCa                                                                          <  2:338596/149‑1 (MQ=255)                                                                         accaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAgg                                                                        <  6:49091/149‑1 (MQ=255)                                                                             tGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTgg                                                                    <  4:37068/149‑1 (MQ=255)                                                                               gctgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGc                                                                   >  1:303725/1‑148 (MQ=255)                                                                                ctgcgctgcTGCGCTAATCAATGGCCAAGCATAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTgg                                                                                         <  2:257330/125‑1 (MQ=255)                                                                                ctgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCCTGAAGCTGGCgt                                                                 <  7:101368/149‑2 (MQ=255)                                                                                ctgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTATTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTgg                                                                                         >  1:257330/1‑125 (MQ=255)                                                                                  gcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAg                                                               >  3:31433/1‑149 (MQ=255)                                                                                   cgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa                                                                <  7:106326/147‑1 (MQ=255)                                                                                               aaccAATGGCCAACCAGAACTGAGTTCTGTTTGATCGGTAAAACGGGCCCCCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGGGGTGCGGCGGATG‑‑GCACTGCTGCAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAg                                                  <  6:234056/146‑1 (MQ=255)                                                                                                aGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCTACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGc                                                 <  4:93979/149‑1 (MQ=255)                                                                                                    aTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCa                                                                                >  3:348713/1‑114 (MQ=255)                                                                                                     tGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGtt                                            >  2:352772/1‑149 (MQ=255)                                                                                                        ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                         <  6:118121/149‑1 (MQ=255)                                                                                                        ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                         >  3:319199/1‑149 (MQ=255)                                                                                                                  cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCa                               <  3:345550/149‑1 (MQ=255)                                                                                                                     aGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAgg                             <  2:303725/148‑1 (MQ=255)                                                                                                                      gttatttttGACCGGAAAAACGGGCGACCCTCGGGGCGCACTTTTTTTTGCTCCGCGACCCGGTGGTGCGTCGGATG‑‑GCACTGCTGGAACACGCCCGCTTCCGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGgc                            <  8:134648/142‑1 (MQ=255)                                                                                                                            tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCa                        <  5:328513/146‑1 (MQ=255)                                                                                                                                  cGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTAt               >  7:167619/1‑149 (MQ=255)                                                                                                                                    gTAAACAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCt             <  5:259927/149‑1 (MQ=255)                                                                                                                                     tAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGAGGACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATc              >  1:107506/1‑149 (MQ=255)                                                                                                                                     aaaaaaGGGCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCa                        <  2:11952/136‑1 (MQ=255)                                                                                                                                     aaaaaaGGGCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCa                        >  1:11952/2‑137 (MQ=255)                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc         >  1:201465/1‑149 (MQ=255)                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATATAAGc         >  7:221859/1‑149 (MQ=255)                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTACGCGACCACGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCCGGTTCACCCGCGAAGCAGGCGCCACTGCCCTATATAAGc         >  6:30657/1‑149 (MQ=255)                                                                                                                                           gggCGACCCGAGGGTCACCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCg                                                           >  3:163010/1‑96 (MQ=255)                                                                                                                                            ggCTACCCTCGGGTTGACTTTTTTTTTGCTCCGCTACCTGGTGGTGCGGCGTTTG‑‑GCACTGCTGTAACACGCCTGTTTCAGGTTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt      <  8:221859/148‑1 (MQ=255)                                                                                                                                            ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg       >  2:266305/1‑147 (MQ=255)                                                                                                                                            ggCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa     <  1:352772/149‑1 (MQ=255)                                                                                                                                               gACCCTCGGGTCGCGCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAc  <  6:98219/149‑1 (MQ=255)                                                                                                                                                  ccTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGc                                <  3:109571/116‑1 (MQ=255)                                                                                                                                                  ccTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGc                                >  4:109571/1‑116 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAC  >  NC_002947/1980945‑1981235

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