Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,789,637	G→C	6.1%	*63S (TGA→TCA)	PP_2441 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,789,637	0	G	C	6.1%	208.0 / 4.5	66	*63S (TGA→TCA)	PP_2441	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/25);  new base C (2/2);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

AAAGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTC  >  NC_002947/2789491‑2789777                                                                                                                                                   |                                                                                                                                             aaaGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGACATCTACGCCGACGGTGGCGCGGTGTGCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACg                                                                                                                                           >  1:1390272/1‑150 (MQ=255)   aGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGc                                                                                                                                         >  1:1019880/1‑150 (MQ=255)      ccGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCa                                                                                                                                      >  1:732054/1‑150 (MQ=255)      atctgcagcaAGATCAGGACCTGCACGGGGGCTACCTGGTAAGCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCa                                                                                                                                      <  1:13282/146‑1 (MQ=255)       cGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAg                                                                                                                                     >  1:586912/1‑150 (MQ=255)           cagcaAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGACCTACTGACCGGGAGCTGGCCGAGCGCAAGACCTGACGCAATTTCGACGGCCTCGAGCCGTGAGGGCTCAgtggt                                                                                                                                 >  2:13282/1‑146 (MQ=255)            agcaAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTc                                                                                                                                >  2:846904/1‑150 (MQ=255)             gcaAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCa                                                                                                                                      >  1:664779/1‑143 (MQ=255)             gcaAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCa                                                                                                                                      <  2:664779/143‑1 (MQ=255)               aaGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTCACGGCTcggg                                                                                                                                    <  1:676416/143‑4 (MQ=255)               aaGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTCACGGCTcggg                                                                                                                                    >  2:676416/1‑140 (MQ=255)                agagCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGc                                                                                                                            >  2:955659/1‑150 (MQ=255)                           tGCACTGGGGATAGCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGgcg                                                                                                                 >  2:1060195/1‑150 (MQ=255)                                tGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCtt                                                                                                            >  2:855173/1‑150 (MQ=255)                                tGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGACCTGGCCGAGCGCAAGAGCTGCCGCAATTTCACCGACCCCGAGCCGTGACGGCTCAGCCTAACAGCTCAGCCGGGCGCGCtt                                                                                                            >  2:387565/1‑150 (MQ=255)                                 ggggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTc                                                                                                           >  1:983582/1‑150 (MQ=255)                                   ggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGc                                                                                                                                    <  2:508380/123‑1 (MQ=255)                                   ggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGc                                                                                                                                    >  1:508380/1‑123 (MQ=255)                                   ggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCa                                                                                                         <  2:1019880/150‑1 (MQ=255)                                          tGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTc                                                                                                                                <  1:48867/120‑1 (MQ=255)                                          tGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTc                                                                                                                                >  2:48867/1‑120 (MQ=255)                                            gAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCaa                                                                                                >  2:1233664/1‑150 (MQ=255)                                             aaaTCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAAt                                                                                               <  1:855173/150‑1 (MQ=255)                                                tttACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCTTACTGACCGGGAGCTGGCCGAGCGCAATAGCTGCCGCAATTTCCCCGTCCCCGAGCCGTGACTGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCg                                                                                            <  2:570667/148‑1 (MQ=255)                                                  tACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCggg                                                                                          <  2:983582/150‑1 (MQ=255)                                                     gCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggctt                                                                                       >  2:602892/1‑150 (MQ=255)                                                     gCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggctt                                                                                       <  2:732054/150‑1 (MQ=255)                                                                       ttCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGt                                                                     >  2:775113/1‑150 (MQ=255)                                                                       ttCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGt                                                                     >  2:788207/1‑150 (MQ=255)                                                                          gtcgACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGtctt                                                                  >  1:711601/1‑150 (MQ=255)                                                                           tcgACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCGATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCTGGCTTCGAGGGCTTCGATGCGGtcttc                                                                 <  1:537695/150‑1 (MQ=255)                                                                            cgACCAGCGCCCTACTGACCGGTAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGtcttct                                                                <  1:782993/150‑1 (MQ=255)                                                                              aCCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCATCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGa                                                              >  1:50474/1‑150 (MQ=255)                                                                                 aGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGacac                                                           >  1:316945/1‑150 (MQ=255)                                                                                 aGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGacac                                                           >  1:1381844/1‑150 (MQ=255)                                                                                       cTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTgcc                                                     >  2:56075/1‑150 (MQ=255)                                                                                        tACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTgccg                                                    <  1:955659/150‑1 (MQ=255)                                                                                          cTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTgccgcc                                                  >  1:142487/1‑150 (MQ=255)                                                                                            gACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGCGCTTCGATGAGGAATTCTGACACACTGCCGCCAc                                                <  2:59281/150‑1 (MQ=255)                                                                                               cGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCggg                                                                                          >  2:1375611/1‑105 (MQ=255)                                                                                               cGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCggg                                                                                          <  1:1375611/105‑1 (MQ=255)                                                                                                  gAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTgc                                          >  2:864946/1‑150 (MQ=255)                                                                                                          cGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGcgccg                                  <  1:788207/150‑1 (MQ=255)                                                                                                             gcgcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCt                               >  1:1036801/1‑150 (MQ=255)                                                                                                             gcgcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCt                               <  1:1060195/150‑1 (MQ=255)                                                                                                               gcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTc                             <  2:1278104/150‑1 (MQ=255)                                                                                                               gcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTc                             <  1:1233664/150‑1 (MQ=255)                                                                                                               gcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTCGTTTCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTc                             >  1:397416/1‑150 (MQ=255)                                                                                                                    agCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCctgct                        >  2:413832/1‑150 (MQ=255)                                                                                                                        gCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCg                    <  1:413832/150‑1 (MQ=255)                                                                                                                         ccGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGt                   >  2:8242/1‑150 (MQ=255)                                                                                                                         ccGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGt                   <  1:864946/150‑1 (MQ=255)                                                                                                                          cGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTg                  >  1:353425/1‑150 (MQ=255)                                                                                                                           gCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggc                                                                                         <  2:504552/78‑1 (MQ=255)                                                                                                                           gCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggc                                                                                         >  1:504552/1‑78 (MQ=255)                                                                                                                           gCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGc                 >  2:1281991/1‑150 (MQ=255)                                                                                                                              aTTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCg                                                                       >  1:28900/1‑93 (MQ=255)                                                                                                                              aTTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCg                                                                       <  2:28900/93‑1 (MQ=255)                                                                                                                                ttCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCga            <  1:56075/150‑1 (MQ=255)                                                                                                                                 tCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCgag           <  1:1196332/150‑1 (MQ=255)                                                                                                                                  ccccgaccccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCgaga          >  2:1272867/1‑150 (MQ=255)                                                                                                                                    ccgaccccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAAt        >  1:1167647/1‑150 (MQ=255)                                                                                                                                         cccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCAAGGGTGCGCGCTGCGCCGTCCTTCTGCGGTGCAGCGAGAATCGCCt   >  2:944167/1‑150 (MQ=255)                                                                                                                                          ctgaGCCGTCACGGCTCGGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTACCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCt   >  1:1261960/3‑149 (MQ=255)                                                                                                                                          ctgaGCCGTCACGGCTCGGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTACCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCt   <  2:1261960/147‑1 (MQ=255)                                                                                                                                          ccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTc  <  1:1281991/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AAAGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTC  >  NC_002947/2789491‑2789777

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