Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,060,559	C→A	7.4%	intergenic (+139/+42)	PP_5311 → / ← PP_5312	signal transduction protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,060,559	0	C	A	7.4%	214.4 / 5.0	67	intergenic (+139/+42)	PP_5311/PP_5312	signal transduction protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base A (3/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCTGAATCCC  >  NC_002947/6060421‑6060699                                                                                                                                           |                                                                                                                                             gcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCca                                                                                                                                   <  1:1028747/150‑1 (MQ=255)  cattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGc                                                                                                                                        <  2:297894/144‑1 (MQ=255)  cattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGc                                                                                                                                        >  1:297894/1‑144 (MQ=255)    ttaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                <  1:241016/150‑1 (MQ=255)       acacATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTa                                                                                                                             >  2:963040/1‑150 (MQ=255)           aTAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGc                                                                                                                         <  2:289508/150‑1 (MQ=255)            tAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCg                                                                                                                        >  2:134674/1‑150 (MQ=255)               ccTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtca                                                                                                                     >  2:867291/1‑150 (MQ=255)               ccTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtca                                                                                                                     <  2:17471/150‑1 (MQ=255)                cTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtcat                                                                                                                    >  1:940724/1‑150 (MQ=255)                          aaCTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtc                                                                                                                      >  2:1184255/1‑138 (MQ=255)                          aaCTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtc                                                                                                                      <  1:1184255/138‑1 (MQ=255)                                 ttGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTc                                                                                                   <  1:618232/150‑1 (MQ=255)                                  tGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCa                                                                                                  <  1:867291/150‑1 (MQ=255)                                       aaTTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCtt                                                                                             <  1:805577/150‑1 (MQ=255)                                           tCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTttctt                                                                                         <  2:1148681/150‑1 (MQ=255)                                               aCAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCGCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCttt                                                                                     >  1:826186/1‑150 (MQ=255)                                                            ttttGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTAtt                                                                        >  2:617446/1‑150 (MQ=255)                                                               tGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGACGGTTGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtca                                                                     <  1:801383/150‑1 (MQ=255)                                                                    tgtgGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtc                                                                                                                      <  1:604011/96‑1 (MQ=255)                                                                    tgtgGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtc                                                                                                                      >  2:604011/1‑96 (MQ=255)                                                                      tgGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAAccc                                                              <  2:1158303/150‑1 (MQ=255)                                                                      tgGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAAccc                                                              >  1:281247/1‑150 (MQ=255)                                                                       gggAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCttt                                                                                     <  1:145912/126‑1 (MQ=255)                                                                       gggAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCttt                                                                                     >  2:145912/1‑126 (MQ=255)                                                                       gggAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAAcccc                                                             >  1:248330/1‑150 (MQ=255)                                                                        ggAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                <  1:14196/82‑1 (MQ=255)                                                                        ggAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                >  2:14196/1‑82 (MQ=255)                                                                             gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAAc                                                       >  2:361667/1‑150 (MQ=255)                                                                             gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGAGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAAc                                                       <  1:617446/150‑1 (MQ=255)                                                                               gCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGCCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACca                                                     >  1:450318/1‑150 (MQ=255)                                                                                  TGCCCGCGAGCACCGGCGAAGCCGGTGCCCTGCACCGCGTCGCTTGTTTCGCGGGCCCGCCCGTTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCCGGTATTTCATCAACCCCTGGAACCACAt                                                  >  2:815563‑M1/1‑150 (MQ=255)                                                                                     ccgcgaaacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCAc                                               >  1:1121891‑M2/55‑150 (MQ=255)                                                                                     ccgcgaaacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaAGCCCGCTCCCACAGTTAAAGCGtcatc                                                                                                                   <  1:710624‑M2/28‑1 (MQ=255)                                                                                     ccgcgaaacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaAGCCCGCTCCCACAGTTAAAGCGtcatc                                                                                                                   >  2:710624‑M2/55‑82 (MQ=255)                                                                                     CCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCt                                                            <  2:230127‑M1/137‑1 (MQ=255)                                                                                     CCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCt                                                            >  1:230127‑M1/1‑137 (MQ=255)                                                                                      cgcgaaacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaCGCGCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTGTTCTTTGGTTTCAGGTATTTCATCAACTCCTGGAACCACATCAcc                                              >  1:442800‑M2/54‑150 (MQ=255)                                                                                      cgcgaaacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaAGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCAcc                                              >  2:216952‑M2/54‑150 (MQ=255)                                                                                       gcgaaacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaAGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCa                                             >  1:348406‑M2/53‑150 (MQ=255)                                                                                             acCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCg                                       <  2:940724/150‑1 (MQ=255)                                                                                                gcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGaa                                                        <  1:119737‑M2/87‑1 (MQ=255)                                                                                                gcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGaa                                                        >  2:119737‑M2/44‑130 (MQ=255)                                                                                                   acgcggtgcatggcaccggcttcgccggtgttcgcgggcaAGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGc                                 <  1:216952‑M2/110‑1 (MQ=255)                                                                                                      gCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCACCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCt                              >  1:223558/1‑150 (MQ=255)                                                                                                         ggTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCttt                                                                                     <  1:904135/92‑1 (MQ=255)                                                                                                         ggTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCttt                                                                                     >  2:904135/1‑92 (MQ=255)                                                                                                         ggTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCGCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTgt                           <  2:826186/150‑2 (MQ=255)                                                                                                           tGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATc                         >  1:385314/1‑150 (MQ=255)                                                                                                            atggcaccggcttcgccggtgttcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCt                        >  1:830610‑M2/32‑150 (MQ=255)                                                                                                             ccATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCtt                                                                                             <  1:947616/80‑1 (MQ=255)                                                                                                             ccATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCtt                                                                                             >  2:947616/1‑80 (MQ=255)                                                                                                                tGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAt                    <  1:134674/150‑1 (MQ=255)                                                                                                                caccggcttcgccggtgttcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAt                    <  2:830610‑M2/123‑1 (MQ=255)                                                                                                                 accggcttcgccggtgttcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAtg                   <  2:1121891‑M2/124‑1 (MQ=255)                                                                                                                  cACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAtgt                  <  2:450318/150‑1 (MQ=255)                                                                                                                  cACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAtgt                  <  1:361667/150‑1 (MQ=255)                                                                                                                   aCCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAtgtg                 >  1:1011607/1‑150 (MQ=255)                                                                                                                    ccGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGc                <  2:385314/150‑1 (MQ=255)                                                                                                                    ccGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGc                <  2:223558/150‑1 (MQ=255)                                                                                                                      gcgTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCtt              <  1:914058/150‑1 (MQ=255)                                                                                                                       cgTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTg             >  1:1090868/1‑150 (MQ=255)                                                                                                                         tCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTc           >  2:628256/1‑150 (MQ=255)                                                                                                                          cGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTcc          >  1:109965/1‑150 (MQ=255)                                                                                                                               gTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGTCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATc                                                 <  2:411313/106‑1 (MQ=255)                                                                                                                               gTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATc                                                 >  1:411313/1‑106 (MQ=255)                                                                                                                                  gcgcgggcaCGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCTGAATccc  >  2:830532‑M2/10‑150 (MQ=255)                                                                                                                                           |                                                                                                                                             GCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCTGAATCCC  >  NC_002947/6060421‑6060699

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