Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	755,882	T→C	5.7%	S159P (TCG→CCG)	PP_0646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	755,882	0	T	C	5.7%	213.7 / 2.0	71	S159P (TCG→CCG)	PP_0646	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/39);  new base C (2/2);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

ACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCATGGCGGCCGCCGAGC  >  NC_002947/755734‑756026                                                                                                                                                     |                                                                                                                                                 acCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCGCCGGAGCATGTCGTTCCCGCAGCTTCTGTTGGAGCTGGCTAGCCATGATTc                                                                                                                                                 >  2:799845/1‑150 (MQ=255) acCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTCGAGCTGGCTAGCCATGATTc                                                                                                                                                 >  2:370719/1‑150 (MQ=255)  cccGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCg                                                                                                                                                >  2:345801/1‑150 (MQ=255)   ccgccgGCCCTGGGTGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTAGCAAATCCTCTCACGTGATTGTCGCGGGATTCTATGTGGTGCCCCGGCGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGa                                                                                                                                               <  1:924769/150‑1 (MQ=255)    cgccgGCCCTGGGAGTTCTAGCGCTTTGGCAGGGGCTTCGACCGCCATGGTTGGCGGGTTGCCTGGCCAGATGCAGCCCGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGAc                                                                                                                                              <  2:168059/150‑1 (MQ=255)     gccgGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGAcc                                                                                                                                             <  1:1681734/150‑1 (MQ=255)              gggAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCAtt                                                                                                                                    <  2:726431/150‑1 (MQ=255)               ggAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCAttt                                                                                                                                   <  1:841602/150‑1 (MQ=255)                gAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTc                                                                                                                                  <  1:345801/150‑1 (MQ=255)                   ctctAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTa                                                                                                                               <  1:1236741/150‑1 (MQ=255)                      tAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCAGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCtt                                                                                                                                           <  2:1239496/135‑1 (MQ=255)                      tAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCAGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCtt                                                                                                                                           >  1:1239496/1‑135 (MQ=255)                          gcTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCaa                                                                                                                        <  2:816320/150‑1 (MQ=255)                              tGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGcc                                                                                                                    <  2:1065533/150‑1 (MQ=255)                                gCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCtt                                                                                                                  <  1:1233138/150‑1 (MQ=255)                                  aGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGACTCGACCTTGCCCATTTCGTACCTGCAACGCCTTaa                                                                                                                <  1:1117646/150‑1 (MQ=255)                                           gcccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGTTTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACttt                                                                                                       <  1:210670/150‑1 (MQ=255)                                            cccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTc                                                                                                      >  1:937710/1‑150 (MQ=255)                                             ccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTAc                                                                                                                              <  2:734667/125‑1 (MQ=255)                                             ccgccATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTAc                                                                                                                              >  1:734667/1‑125 (MQ=255)                                                 cATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCa                                                                                                 >  2:84121/1‑150 (MQ=255)                                                    ccatGGCGGGCTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAggg                                                                                              >  1:1920647/4‑150 (MQ=255)                                                     gTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAgggg                                                                                             <  2:937710/150‑1 (MQ=255)                                                      ttGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGt                                                                                            >  2:524685/1‑150 (MQ=255)                                                            ggTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAgg                                                                                      <  1:1029572/150‑1 (MQ=255)                                                                gCCTGGGCAGATGCAGCCGGGATTCTATGTGGCGCACCGTAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAgg                                                                                  <  1:593276/150‑1 (MQ=255)                                                                 ccTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAggg                                                                                 >  1:1853561/1‑150 (MQ=255)                                                                     ggCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCACCTTCTGTTGGACCAAAATCTGCATGATTCGACCGTGCCCATTTCGTACCTTCAACGCCTTAACCCGACGTTCGAGCAGGGGTTCAAGGCAGGGGaaa                                                                             <  1:1304162/150‑1 (MQ=255)                                                                      gCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAg                                                                                                >  1:183890/1‑130 (MQ=255)                                                                      gCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCAGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAg                                                                                                <  2:183890/130‑1 (MQ=255)                                                                          aTGCAGCCGTGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTc                                                                        <  1:84121/150‑1 (MQ=255)                                                                          aTGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTc                                                                        <  1:1615885/150‑1 (MQ=255)                                                                          aTGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGCTCTCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACATACACCGCTATCTCCCGACTTTCGAGCAGGCGTTCAAGGCCGGGGAGATCTTc                                                                        <  2:761182/150‑1 (MQ=255)                                                                                ccGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATc                                                                  <  1:703896/150‑1 (MQ=255)                                                                                 cGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCg                                                                 >  1:955305/1‑150 (MQ=255)                                                                                  gggATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCgg                                                                >  2:155074/1‑150 (MQ=255)                                                                                           tgtgGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGa                                                       <  2:955305/150‑1 (MQ=255)                                                                                            gtggtgCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCt                                                                                                                            >  1:1677686/1‑80 (MQ=255)                                                                                            gtggtgACCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCt                                                                                                                            <  2:1677686/80‑1 (MQ=255)                                                                                                tgCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGATCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGCCCTCTGCCCACTGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACg                                                  >  1:1574663/1‑150 (MQ=255)                                                                                                   cccGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGtt                                               >  2:608537/1‑150 (MQ=255)                                                                                                     cGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAggg                                                                                 >  1:1436637/1‑114 (MQ=255)                                                                                                     cGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAggg                                                                                 <  2:1436637/114‑1 (MQ=255)                                                                                                     cGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGACCAGGTGTTCAAGTCAGTGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTAt                                             <  1:530943/150‑1 (MQ=255)                                                                                                       gAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           <  1:524685/150‑1 (MQ=255)                                                                                                        aGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGcc                                          <  1:1126954/150‑1 (MQ=255)                                                                                                                gTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCa                                  <  2:1208528/150‑1 (MQ=255)                                                                                                                 ttCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           <  1:1838243/140‑1 (MQ=255)                                                                                                                 ttCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           >  2:1838243/1‑140 (MQ=255)                                                                                                                 ttCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCag                                 >  2:608232/1‑150 (MQ=255)                                                                                                                   cccGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCagga                               <  2:1853561/150‑1 (MQ=255)                                                                                                                    ccGCAGATTCTGTTGGATCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggag                              >  1:383380/1‑150 (MQ=255)                                                                                                                    ccGCAGATTCTGTTGGATCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggag                              >  1:936755/1‑150 (MQ=255)                                                                                                                     cGCAGCTTCTGTTGGAGCTGGCTAGCCATGCTGCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACGTTAGAGCAGGGCTTTAAGGCAGGGGAAATCTTCGTCATCGGCGAACCGCATAACGGTTATGCCGGCACCCcggcgg                             >  1:41325/1‑150 (MQ=255)                                                                                                                     cGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggagg                             <  1:512490/150‑1 (MQ=255)                                                                                                                     cGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGACGTTGGCGATCCGGATAACGGTTATGCCTGCACCCaggagg                             >  1:509937/1‑150 (MQ=255)                                                                                                                      gCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCaggagga                            <  1:608232/150‑1 (MQ=255)                                                                                                                           ttCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGcc                       >  2:1748189/1‑150 (MQ=255)                                                                                                                              tGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAAc                                                   <  1:832501/119‑1 (MQ=255)                                                                                                                              tGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAAc                                                   >  2:832501/1‑119 (MQ=255)                                                                                                                               gTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACt                   >  2:456026/1‑150 (MQ=255)                                                                                                                                ttGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTc                  >  2:1290846/1‑150 (MQ=255)                                                                                                                                 tGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCa                 <  1:456026/150‑1 (MQ=255)                                                                                                                                    aGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCATgg              <  1:370719/150‑1 (MQ=255)                                                                                                                                              ccATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCATGGCGGCCGCCGa    >  1:912119/1‑150 (MQ=255)                                                                                                                                               cATGATCCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           <  1:861493/110‑1 (MQ=255)                                                                                                                                               cATGATCCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           >  2:861493/1‑110 (MQ=255)                                                                                                                                               cATGATCCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           >  2:861546/1‑110 (MQ=255)                                                                                                                                               cATGATCCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGc                                           <  1:861546/110‑1 (MQ=255)                                                                                                                                                aTGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCATGGCGGCCGCCGAgc  <  1:608537/150‑1 (MQ=255)                                                                                                                                                aTGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCATGGCGGCCGCCGAgc  <  1:1368045/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 ACCCGCCGGCCCTGGGAGCTCTAGCGCTTTGGCAGGGGCTTCGCCCGCCATGGTTGGCGGGTTGCCTGGGCAGATGCAGCCGGGATTCTATGTGGTGCCCCGGAGCATGTCGTTCCCGCAGATTCTGTTGGAGCTGGCTAGCCATGATTCGACCTTGCCCATTTCGTACCTGCAACGCCTTAACCCGACTTTCGAGCAGGGGTTCAAGGCAGGGGAAATCTTCGTCATCGGCGATCCGGATAACGGTTATGCCTGCACCCAGGAGGAGGGCCAACTCATGGCGGCCGCCGAGC  >  NC_002947/755734‑756026

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