Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,928,017	A→G	7.0%	S719P (TCC→CCC)	cheA ←	chemotaxis histidine kinase CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,928,017	0	A	G	7.0%	198.9 / 4.8	71	S719P (TCC→CCC)	cheA	chemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base G (3/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CCTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCAC  >  NC_002947/4927869‑4928161                                                                                                                                                      |                                                                                                                                                  ccTTCACTGCCATAAACACTCCATTCGGCGGCCCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGa                                                                                                                                                   <  5:264031/149‑1 (MQ=255)  cTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGCGCGCGCCACCGAAATCAAATACCCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATCCGACCGTCGCCGGTAATGGTGGCCCCGGa                                                                                                                                                   >  7:125178/1‑148 (MQ=255)    tGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACAt                                                                                                                                                >  8:94018/1‑149 (MQ=255)           cATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCg‑gg                                                                                                                                        <  8:55941/149‑1 (MQ=255)            aTAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCg‑ggg                                                                                                                                       <  2:300530/149‑1 (MQ=255)             tAAACACTCCATTAGTCGGAGCAACCAGCCAGGGCGCGCCACCGAACTCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGt                                                                                                                                      <  1:264599/149‑1 (MQ=255)             tAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGGCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGGCATGTCCg‑gggc                                                                                                                                      <  5:63616/149‑2 (MQ=255)              aaaCACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGCTTCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGt                                                                                                                                      <  1:117935/149‑1 (MQ=255)                    tCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTg                                                                                                                                >  2:384596/1‑148 (MQ=255)                      cATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTgcag                                                                                                                             <  5:306940/149‑1 (MQ=255)                      cATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTgca                                                                                                                              <  4:40330/148‑1 (MQ=255)                       aTTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTgcag                                                                                                                             >  5:33014/1‑148 (MQ=255)                        ttAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGAACATGCCCG‑GGGTGCCTTgcagca                                                                                                                           >  7:167011/1‑149 (MQ=255)                           ggcggcGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGGCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATc                                                                                                                         <  7:253594/148‑1 (MQ=255)                            gcggcgCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCtt                                                                                                                       >  7:82730/1‑149 (MQ=255)                            gcggcgCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCtt                                                                                                                       >  1:255591/1‑149 (MQ=255)                             cggcgCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTg                                                                                                                      <  7:94018/149‑1 (MQ=255)                                 gcAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGcc                                                                                                                    >  8:156291/1‑147 (MQ=255)                                     cccGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCgg                                                                                                              >  7:148355/1‑149 (MQ=255)                                      ccGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGc                                                                                                             <  1:384596/149‑1 (MQ=255)                                      ccGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGGCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGc                                                                                                             <  8:82730/149‑1 (MQ=255)                                         ccAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTtg                                                                                                          <  7:203388/149‑1 (MQ=255)                                          cAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCGCGGGTGCCTTGCAGCATCTTGCCCAGCGGCTtg                                                                                                          <  1:246165/149‑1 (MQ=255)                                           aGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGAACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTtgat                                                                                                        <  8:167011/149‑1 (MQ=255)                                                gcgcCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCa                                                                                                   >  6:260203/1‑149 (MQ=255)                                                    cACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCActtc                                                                                               >  5:319234/1‑149 (MQ=255)                                                     aCCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCActtct                                                                                              <  5:260203/149‑1 (MQ=255)                                                       cGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACttcttc                                                                                            >  1:179954/1‑149 (MQ=255)                                                           aTCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTgg                                                                                        <  8:148355/149‑1 (MQ=255)                                                            tCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGc                                                                                       <  8:319163/149‑1 (MQ=255)                                                            tCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGc                                                                                       <  8:269953/149‑1 (MQ=255)                                                              aaaTACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGccc                                                                                     >  6:281087/1‑149 (MQ=255)                                                              aaaTACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGcc                                                                                      >  2:255609/1‑148 (MQ=255)                                                              aaaTACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGcc                                                                                      >  8:182763/1‑148 (MQ=255)                                                                 tACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCccacc                                                                                  <  1:249500/149‑1 (MQ=255)                                                                   cgccgcGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCccacca                                                                                 <  7:182763/148‑1 (MQ=255)                                                                   cgccgcGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCccacca                                                                                 <  6:140345/148‑1 (MQ=255)                                                                    gccgcGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGt                                                                               <  1:316122/149‑1 (MQ=255)                                                                    gccgcGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGt                                                                               >  3:339558/1‑149 (MQ=255)                                                                       gcgcgGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTgg                                                                            >  1:201262/1‑149 (MQ=255)                                                                             cgTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTcgacg                                                                      <  2:178627/149‑1 (MQ=255)                                                                                  cGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACaaa                                                                 <  6:319234/149‑1 (MQ=255)                                                                                       gAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCg                                                             <  5:199194/148‑1 (MQ=255)                                                                                        aGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGAt                                                           <  6:33014/149‑1 (MQ=255)                                                                                                 cGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTggg                                                   >  8:178065/1‑149 (MQ=255)                                                                                                  gggACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTggg                                                   <  1:327393/148‑1 (MQ=255)                                                                                                  gggACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTgg                                                    <  7:156291/147‑1 (MQ=255)                                                                                                      cGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGccg                                              <  2:255591/149‑1 (MQ=255)                                                                                                       gTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGccga                                             >  1:44701/1‑149 (MQ=255)                                                                                                          gagaATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTgg                                                    <  3:378229/139‑1 (MQ=255)                                                                                                          gagaATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTgg                                                    >  4:378229/1‑139 (MQ=255)                                                                                                          gagaATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGccgaccg                                          >  6:189176/1‑149 (MQ=255)                                                                                                             aaTCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACa                                       <  5:281087/149‑1 (MQ=255)                                                                                                               tCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTgg                                                    <  4:183557/134‑1 (MQ=255)                                                                                                               tCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTgg                                                    >  3:183557/1‑134 (MQ=255)                                                                                                                 aaCGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGAGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAAt                                   >  8:344503/1‑149 (MQ=255)                                                                                                                  actcgATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGacagacagcag                                   >  8:167236/4‑145 (MQ=255)                                                                                                                  aCGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAAt                                   >  2:87507/1‑148 (MQ=255)                                                                                                                  aCGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATc                                  <  2:179954/149‑1 (MQ=255)                                                                                                                   cgcgATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATca                                 <  7:178065/149‑1 (MQ=255)                                                                                                                   cgcgATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATca                                 >  2:79930/1‑149 (MQ=255)                                                                                                                      caTGCGACCGTCGCCGGTAATGGTGGCCCAGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATcacca                              >  4:131574/2‑149 (MQ=255)                                                                                                                         gCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGt                           >  2:222063/1‑149 (MQ=255)                                                                                                                          cGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTg                          <  2:201262/149‑1 (MQ=255)                                                                                                                                  gCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCg                   <  3:131574/148‑1 (MQ=255)                                                                                                                                   ccGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGtg                 <  2:44701/149‑1 (MQ=255)                                                                                                                                      gTAATGGTGGCCCCGGGCATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGtgctg              >  3:115177/1‑149 (MQ=255)                                                                                                                                      gTAATGGTGGCCCCGGGCATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGtgctg              <  4:115177/149‑1 (MQ=255)                                                                                                                                         caaggcaccccGGGCATGTCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGctc           >  6:382317/8‑149 (MQ=255)                                                                                                                                            gtggCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTc        >  5:6463/1‑149 (MQ=255)                                                                                                                                               accccGGGCATGTCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTg     >  4:97424/2‑149 (MQ=255)                                                                                                                                                  ccGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCAc  >  5:52606/1‑149 (MQ=255)                                                                                                                                                  ccGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCAc  >  8:239885/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCTTGACTGCCATAAACACTCCATTAGGCGGCGCAACCCGCCAGGGCGCGCCACCGAAATCAAATACGCCGCGCGGCGTAACGCTTGAGCATGC‑TCGGGACGTCGAGAATCAACGCGATGCGACCGTCGCCGGTAATGGTGGCCCCGGACATGCCCG‑GGGTGCCTTGCAGCATCTTGCCCAGCGGCTTGATGACCACTTCTTCCTGGCCCACCAGTTGGTCGACGACAAAACCGATGCGCTGGGTGCCGACCGACAGAATCACCACGTGGCCCTCGTGCTGCTCTTCGTGCAC  >  NC_002947/4927869‑4928161

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: ‑