Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,981,092	T→G	7.8%	intergenic (+31/‑269)	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,092	0	T	G	7.8%	183.4 / 7.4	64	intergenic (+31/‑269)	rpsA/ihfB	30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base T (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 = 6.83e-01

ATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAC  >  NC_002947/1980944‑1981235                                                                                                                                                     |                                                                                                                                                  aTCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCt                                                                                                                                                   >  2:357574/1‑149 (MQ=255)  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                                                                                                                                  <  6:134639/141‑1 (MQ=255)                          gAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCttttttt                                                                                                                                  >  5:134639/1‑141 (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                        >  1:11952/2‑137 (MQ=255)                                                                                                                                      aaaaaaGGGCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCa                        <  2:11952/136‑1 (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)                                                                                                                                                     |                                                                                                                                                  ATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATG‑‑GCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAC  >  NC_002947/1980944‑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: ‑