Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,951,546	G→T	6.1%	R243S (CGC→AGC)	fliK ←	flagellar hook‑length control protein FliK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,951,546	0	G	T	6.1%	210.1 / 3.3	66	R243S (CGC→AGC)	fliK	flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base T (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CCCGGTTGACCAGGCCCTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCTTCCTCGCCC  >  NC_002947/4951399‑4951691                                                                                                                                                    |                                                                                                                                                    cccGGTTGACCAGGCCCTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCgg                                                                                                                                                    >  6:201515/1‑149 (MQ=255)  ccGGTTGACCAGGCCCTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCggg                                                                                                                                                   >  7:150053/1‑149 (MQ=255)         aCCAGGCCCTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCgct                                                                                                                                            >  1:237315/1‑149 (MQ=255)           cAGGCCCTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCgctgc                                                                                                                                          >  3:224497/1‑149 (MQ=255)                ccTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGgcg                                                                                                                                     >  6:222063/1‑149 (MQ=255)                ccTCGGCCCAAGCGTTCTGGTTAAACGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTATTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGgcg                                                                                                                                     >  7:239243/1‑149 (MQ=255)                  tCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGgtgt                                                                                                                                                <  3:64511/136‑1 (MQ=255)                  tCGGCCCAAGCGTTCTGGTTCATCGGCAATGCCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGgtgt                                                                                                                                                >  4:64511/1‑136 (MQ=255)                     gCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCtt                                                                                                                                <  5:62940/149‑1 (MQ=255)                          aGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGgtg                                                                                                                            >  6:95332/1‑148 (MQ=255)                                   gTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAgg                                                                                                                   <  6:35699/148‑1 (MQ=255)                                   gTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGc                                                                                                                  >  6:239474/1‑149 (MQ=255)                                     tCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGgt                                                                                                                             <  8:185809/136‑1 (MQ=255)                                     tCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGgt                                                                                                                             >  7:185809/1‑136 (MQ=255)                                       aTCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTc                                                                                                              <  8:150053/149‑1 (MQ=255)                                        tCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCt                                                                                                             >  2:241414/1‑149 (MQ=255)                                          ggCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGGGTCTTTCAGGCCGTCTTc                                                                                                           <  5:201515/149‑1 (MQ=255)                                           gCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTca                                                                                                          <  8:240738/149‑1 (MQ=255)                                            cAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGgt                                                                                                                                                  >  4:222221/1‑108 (MQ=255)                                             aaTGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTcaag                                                                                                        <  8:239243/149‑1 (MQ=255)                                                   tGATGCAGCGGGTTCGCAGTGGCAGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAgg                                                                                                   >  2:344314/1‑148 (MQ=255)                                                    gATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCa                                                                                                 >  2:60807/1‑149 (MQ=255)                                                     aTGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCa                                                                                                 <  5:361263/148‑1 (MQ=255)                                                      tGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCAc                                                                                                >  2:182110/1‑148 (MQ=255)                                                                  gCAGTGGGCTGTACGGCGTTGGCGGTATTGGCGGTGGGCGCCTGTGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGc                                                                                    <  5:142271/148‑1 (MQ=255)                                                                   cAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTc                                                                                                              <  7:203276/121‑1 (MQ=255)                                                                   cAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTc                                                                                                              >  8:203276/1‑121 (MQ=255)                                                                     gTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑Acc                                                                              <  5:239474/149‑1 (MQ=255)                                                                        gCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGTGTGTCGCTGCTGGAGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGa                                                                           <  5:267695/149‑1 (MQ=255)                                                                        gCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGGTCTTGGTGTCTTTCAGGCCGGCTTCAAGCAAGGCACCAAAGGACTTGCC‑‑ACCGGa                                                                           >  8:325736/1‑149 (MQ=255)                                                                            gTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCa                                                                       <  4:224497/149‑1 (MQ=255)                                                                             tACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAg                                                                      <  1:241414/149‑1 (MQ=255)                                                                             tACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAAGCCTTGCC‑‑ACCGGATTCAg                                                                      >  2:153528/1‑149 (MQ=255)                                                                                 gCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCt                                                                  >  1:339551/1‑149 (MQ=255)                                                                                  cGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAg                                                                      >  6:312718/1‑144 (MQ=255)                                                                                  cGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAg                                                                      <  5:312718/144‑1 (MQ=255)                                                                                  cGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCtt                                                                 >  3:190090/1‑149 (MQ=255)                                                                                  cGTTCGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCCAGCAAGGCACCAAAGGCCATGCC‑‑ACCGGATTCAGCGCtt                                                                 >  6:71182/1‑149 (MQ=255)                                                                                                 ccgtggccgCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATcc                                                  >  8:206241/1‑149 (MQ=255)                                                                                                 ccgtggccgCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATc                                                   >  2:142914/1‑148 (MQ=255)                                                                                                   gtggccgCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCTGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTc                                                >  2:283399/1‑149 (MQ=255)                                                                                                   gtggccgCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCTGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTc                                                >  2:288153/1‑149 (MQ=255)                                                                                                     ggccgCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTc                                                <  8:175170/147‑1 (MQ=255)                                                                                                     ggccgCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTc                                                >  7:175170/1‑147 (MQ=255)                                                                                                      gccgccTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGAc                                             <  1:344314/149‑1 (MQ=255)                                                                                                      gccgccTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGAc                                             <  1:182110/149‑1 (MQ=255)                                                                                                       ccgccTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGAc                                             >  1:185035/1‑148 (MQ=255)                                                                                                        cgccTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGcc                                                            <  1:207430/132‑1 (MQ=255)                                                                                                        cgccTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGcc                                                            >  2:207430/1‑132 (MQ=255)                                                                                                        cgccTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGAc                                             >  4:138790/1‑147 (MQ=255)                                                                                                              ggTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTc                                      <  2:237315/148‑1 (MQ=255)                                                                                                                  aGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCaa                                 <  5:222063/149‑1 (MQ=255)                                                                                                                    gCTGGCCAGGCGGTTGGCAAAATCATCGGTGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAgg                               >  1:22006/1‑149 (MQ=255)                                                                                                                      tGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTa                             <  1:142914/149‑1 (MQ=255)                                                                                                                        gCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCACGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTAtt                           >  3:9675/1‑149 (MQ=255)                                                                                                                          cAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTAttaa                         >  5:53376/1‑147 (MQ=255)                                                                                                                          cAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGa                         <  5:95332/149‑1 (MQ=255)                                                                                                                               ggTTGGCAAAATCATCGATGCGGGTGTCGATGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGgcc                    >  2:99315/1‑149 (MQ=255)                                                                                                                                 ttGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGgccgc                  <  4:190090/149‑1 (MQ=255)                                                                                                                                     cAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGc              >  1:286647/1‑149 (MQ=255)                                                                                                                                       aaaTCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGccc            <  7:206241/149‑1 (MQ=255)                                                                                                                                        aaTCATCGATGCTGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCt           <  1:288153/149‑1 (MQ=255)                                                                                                                                        aaTCATCGATGCTGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCt           <  1:283399/149‑1 (MQ=255)                                                                                                                                        aaTCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCt           >  3:173257/1‑149 (MQ=255)                                                                                                                                          tcatcGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCTTc         <  7:21967/149‑1 (MQ=255)                                                                                                                                          tcatcGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCTTc         <  2:22006/149‑1 (MQ=255)                                                                                                                                            atcGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCTTCCt       >  7:63904/1‑149 (MQ=255)                                                                                                                                                   cGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCCTGAGCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCTTCCTCGccc  >  1:307622/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CCCGGTTGACCAGGCCCTCGGCCCAAGCGTTCTGGTTCATCGGCAATGGCTGATGCAGCGGGTTCGCAGTGGCCGGTACGGCGTTGGCGGTCTTGGCCGTGGCCGCCTGGGTCAGGCTGGCCAGGCGGTTGGCAAAATCATCGATGCGGGTGTCGCTGCTGGCGCTCTTGGTGTCTTTCAGGCCGTCTTCAAGCAAGGCACCAAAGGCCTTGCC‑‑ACCGGATTCAGCGCTTTCGCCATCGGTATCCTCGACCAGATTCGCCAAGGTATTGACGGCCGCCTGCCCTTCCTCGCCC  >  NC_002947/4951399‑4951691

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