Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,206,885	C→G	5.1%	H91Q (CAC→CAG)	PP_1950 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,206,885	0	C	G	5.1%	252.7 / 4.3	78	H91Q (CAC→CAG)	PP_1950	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base G (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

ATCAGGCTTACTGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATGTATTGA  >  NC_002947/2206737‑2207033                                                                                                                                                      |                                                                                                                                                     atCAGGCTTCCTGCACGGGCAGTCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACAcc                                                                                                                                                     <  1:1169157/150‑1 (MQ=255)      gCTTACTGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGAc                                                                                                                                                >  2:1668260/1‑150 (MQ=255)        ttACTGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGAccc                                                                                                                                              >  2:1700403/1‑150 (MQ=255)          aCTGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTg                                                                                                                                            >  2:25090/1‑150 (MQ=255)           cTGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGt                                                                                                                                           >  1:206238/1‑150 (MQ=255)            tGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCATCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTc                                                                                                                                          >  1:805136/1‑150 (MQ=255)               aCGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCG‑ACACGCTGCCACTTGATGACGCTGCACACCGTGACCCTGTCATg                                                                                                                                       >  1:1374004/1‑150 (MQ=255)                 gggCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGaa                                                                                                                                     <  2:657116/150‑1 (MQ=255)                 gggCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGaa                                                                                                                                     <  2:1155491/150‑1 (MQ=255)                      ggcaTGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCACTATACTTGTTGCTTAGGTT‑CCGCCGGAGCATGCTGCACACCGTGACCCTGTCATGAAgttgt                                                                                                                                >  2:990571/1‑150 (MQ=255)                       gcaTGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACAcc                                                                                                                                                     <  2:1066569/128‑1 (MQ=255)                       gcaTGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAgttgtt                                                                                                                               >  1:1029936/1‑150 (MQ=255)                       gcaTGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAgttgtt                                                                                                                               <  1:1596844/150‑1 (MQ=255)                       gcaTGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGAACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACAcc                                                                                                                                                     >  1:1066569/1‑128 (MQ=255)                       gcaTGACATGACGCAGAAAGTCTTGAAGGACTGGAAAACATTGACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCAGAACACCGTGCCCCGGTAATGAAgttggt                                                                                                                               >  2:425835/1‑148 (MQ=255)                         aTGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTa                                                                                                                             >  2:189286/1‑150 (MQ=255)                           gACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACt                                                                                                                           <  1:454524/150‑1 (MQ=255)                                gACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGt                                                                                                                      <  2:1244573/150‑1 (MQ=255)                                    cAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATAAAGGTGTTTACTCCCGTGGCa                                                                                                                  >  1:1531282/1‑150 (MQ=255)                                      aaaGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCAtt                                                                                                                <  1:279049/150‑1 (MQ=255)                                        aGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGTCATTGt                                                                                                              <  1:129090/150‑1 (MQ=255)                                           tCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGa                                                                                                           >  1:1117495/1‑150 (MQ=255)                                             ttGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCgtac                                                                                                         <  2:1222103/150‑4 (MQ=255)                                               gAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGt                                                                                                                      <  1:721067/135‑1 (MQ=255)                                               gAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGt                                                                                                                      >  2:721067/1‑135 (MQ=255)                                                     cTGGAAAACATTTGCCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGa                                                                                                 >  2:1257113/1‑150 (MQ=255)                                                      tGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGTGT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGAt                                                                                                >  1:1074955/1‑150 (MQ=255)                                                      tGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGt                                                                                                                                           >  2:1094279/1‑107 (MQ=255)                                                      tGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGt                                                                                                                                           <  1:1094279/107‑1 (MQ=255)                                                      tGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGGCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGAt                                                                                                >  2:1383814/1‑150 (MQ=255)                                                        gAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTa                                                                                              >  2:1583535/1‑150 (MQ=255)                                                        gAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTa                                                                                              >  2:82496/1‑150 (MQ=255)                                                        gAAAACATTTAAATCGACGCTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTa                                                                                              >  2:1912122/1‑150 (MQ=255)                                                              aTTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGaa                                                                                                                                     <  2:1555978/105‑1 (MQ=255)                                                              aTTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGaa                                                                                                                                     >  1:1555978/1‑105 (MQ=255)                                                              aTTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGa                                                                                        <  1:665947/150‑1 (MQ=255)                                                                       gACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATgcgc                                                                               <  1:82496/150‑1 (MQ=255)                                                                            tAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGaaa                                                                                                         >  2:1303674/1‑119 (MQ=255)                                                                            tAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGaaa                                                                                                         <  1:1303674/119‑1 (MQ=255)                                                                                 gCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCtg                                                                     <  1:1810110/150‑1 (MQ=255)                                                                                      tGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGa                                                                <  2:1846857/150‑1 (MQ=255)                                                                                       ggCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGAt                                                               <  1:1000851/150‑1 (MQ=255)                                                                                       ggCCCTCGCGAGCACATTCCCAACATCCTTGTCGCAGAGGAT‑CCGCCGGCTCACGCTGCCCACCGTGCCCCTGTCATGAAGTTGTTTACTCCCGTTGCATTGTCGAAATATCAAGCTTACTTCGATAAATGCGCCGAAATCCTGGTGGAt                                                               <  2:846543/150‑1 (MQ=255)                                                                                         cccTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGa                                                             <  1:1700403/150‑1 (MQ=255)                                                                                           cTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATACATGCGCCGAAATCCTGGTGGATGAcc                                                           >  1:246097/1‑150 (MQ=255)                                                                                            tCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCt                                                          >  1:824698/1‑150 (MQ=255)                                                                                            tCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCt                                                          >  2:143396/1‑150 (MQ=255)                                                                                            tCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCt                                                          >  1:1595123/1‑150 (MQ=255)                                                                                               cgAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTtgt                                                       <  1:25090/150‑1 (MQ=255)                                                                                                    acaTTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTgcggc                                                  >  2:890339/1‑150 (MQ=255)                                                                                                          ccAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGt                                            >  1:729201/1‑150 (MQ=255)                                                                                                           cAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTa                                           >  2:203336/1‑150 (MQ=255)                                                                                                              tatCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTAtgt                                        <  1:143396/150‑1 (MQ=255)                                                                                                               atCCTTTTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTAtgtg                                       <  2:295590/150‑1 (MQ=255)                                                                                                                tCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTgg                                      >  2:1555357/1‑150 (MQ=255)                                                                                                                  cTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGTGCCGAAATCCTGGTGGGTGGCCTTGTGCGGACTCAGTATGTGGAt                                    <  2:805136/150‑1 (MQ=255)                                                                                                                  cTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGAt                                    <  1:1912122/150‑1 (MQ=255)                                                                                                                  cTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGAt                                    <  2:230351/150‑1 (MQ=255)                                                                                                                  cTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGAt                                    <  2:1870539/150‑1 (MQ=255)                                                                                                                   ttGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATg                                   >  1:633049/1‑150 (MQ=255)                                                                                                                   ttGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTGTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATg                                   >  2:643474/1‑150 (MQ=255)                                                                                                                    tGTCGCAGAGGAT‑CCTCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATgg                                  <  2:457597/150‑1 (MQ=255)                                                                                                                      taaacaacttcatgacagggTCACGGTGTGCAGCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATgcg                                                                                >  2:1220765‑M1/21‑103 (MQ=255)                                                                                                                      taaacaacttcatgacagggTCACGGTGTGCAGCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATgcg                                                                                <  1:1220765‑M1/83‑1 (MQ=255)                                                                                                                      tCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCt                                <  2:729201/150‑1 (MQ=255)                                                                                                                        gCAGAGGAT‑CCGCCGGATCACGCTGCACCCCGTGACCCTGTCAAGGAGTCGTTTACGCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCttt                              >  1:712028/1‑150 (MQ=255)                                                                                                                        aacaacttcatgacagggTCACGGTGTGCAGCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCtg                                                                     >  1:366312‑M1/19‑112 (MQ=255)                                                                                                                        aacaacttcatgacagggTCACGGTGTGCAGCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCtg                                                                     <  2:366312‑M1/94‑1 (MQ=255)                                                                                                                         cAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCtttt                             <  1:1583535/150‑1 (MQ=255)                                                                                                                                  ccgccgGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTAGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGtt                     <  2:1029936/150‑1 (MQ=255)                                                                                                                                    gccgGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTTACGTtgc                   >  1:837744/1‑150 (MQ=255)                                                                                                                                     ccgGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTtgct                  >  2:1292404/1‑150 (MQ=255)                                                                                                                                      cgGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTtgctg                 >  1:893970/1‑150 (MQ=255)                                                                                                                                      cgGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCAGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGAGGCGTGAGTATGTGGATGGCTTTTCTGACGTtgctc                 >  1:641531/1‑149 (MQ=255)                                                                                                                                      cgGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCACGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTtgctg                 >  1:1085949/1‑150 (MQ=255)                                                                                                                                        gATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCa               <  1:1292404/150‑1 (MQ=255)                                                                                                                                          TCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAg             <  1:567722‑M1/150‑1 (MQ=255)                                                                                                                                            ACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGtt           <  2:1343143‑M1/150‑1 (MQ=255)                                                                                                                                              GCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTAt         <  1:1444345‑M1/150‑1 (MQ=255)                                                                                                                                               cTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATg        <  1:1257113/150‑1 (MQ=255)                                                                                                                                                tGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATGt       >  2:1643911/1‑150 (MQ=255)                                                                                                                                                   acacCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTtgtg                                                      <  1:617871/100‑1 (MQ=255)                                                                                                                                                   acacCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTtgtg                                                      >  2:617871/1‑100 (MQ=255)                                                                                                                                                   acacCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATGTAtt    >  1:1518564/1‑150 (MQ=255)                                                                                                                                                     acCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATGTATTGa  <  2:633049/150‑1 (MQ=255)                                                                                                                                                     acCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATGTATTGa  >  1:1464323/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ATCAGGCTTACTGCACGGGCAGGCATGACCTGACGCAAAAGGTCTTGAAGGACTGGAAAACATTTACCTCGACGGTAAAGGCTTTTGGCCCTCGCGAGCACATTCCCAATATCCTTGTCGCAGAGGAT‑CCGCCGGATCACGCTGCACACCGTGACCCTGTCATGAAGTTGTTTACTCCCGTGGCATTGTCGAAATATCAAGATTACTTCGATAAATGCGCCGAAATCCTGGTGGATGACCTTGTGCGGCGTGAGTATGTGGATGGCTTTTCTGACGTTGCTGCAAGTTATGTATTGA  >  NC_002947/2206737‑2207033

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