Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,154,396	C→A	7.7%	G126W (GGG→TGG)	PP_1909 ←	septum formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,154,396	0	C	A	7.7%	183.7 / 6.4	65	G126W (GGG→TGG)	PP_1909	septum formation protein
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base A (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

TGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGA  >  NC_002947/2154251‑2154542                                                                                                                                                   |                                                                                                                                                   tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGt                                                                                                                                                 >  6:183706/1‑149 (MQ=255) tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGt                                                                                                                                                 >  3:86099/1‑149 (MQ=255)  ggAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                                >  8:13358/1‑149 (MQ=255)    aaCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgc                                                                                                                                              <  5:142360/149‑1 (MQ=255)      cAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcgc                                                                                                                                            >  5:346927/1‑149 (MQ=255)           tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGTggatt                                                                                                                                            >  4:128932/1‑141 (MQ=255)           tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGTggatt                                                                                                                                            <  3:128932/144‑4 (MQ=255)             aCCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                                                                     <  2:301825/149‑1 (MQ=255)               ccccAGCCCCTCAGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTACTGCTGCCACGTAACGCTCGACC‑TGCTCCATGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATACACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                                                                     >  4:205015/1‑147 (MQ=255)                     ccccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCa                                                                                                                             <  7:13358/149‑1 (MQ=255)                     ccccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCa                                                                                                                             <  4:25737/149‑1 (MQ=255)                        cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGcc                                                                                                                          >  1:194989/1‑149 (MQ=255)                           cGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGccc                                                                                                                         <  3:129994/147‑1 (MQ=255)                           cGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                                       <  5:235787/149‑1 (MQ=255)                           cGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                                       >  6:283954/1‑149 (MQ=255)                              cTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTaa                                                                                                                    <  5:183706/149‑1 (MQ=255)                               ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGt                                                                                                                      >  2:191288/1‑146 (MQ=255)                                 gAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGa                                                                                                                  >  7:126675/1‑148 (MQ=255)                                  aaccTCCCTTCGCAATCCCGCCGCTGTTCTGCTGCCACGTAACGCGCTAAC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCTCCTGGCATTGCCCGCTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                                <  8:304120/146‑1 (MQ=255)                                   aaCTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAg                                                                                                               <  4:71118/149‑1 (MQ=255)                                     cTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGc                                                                                                              >  8:124680/1‑148 (MQ=255)                                        ccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                                <  5:208904/143‑1 (MQ=255)                                        ccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                                >  6:208904/1‑143 (MQ=255)                                        ccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCAc                                                                                                          <  8:126675/149‑1 (MQ=255)                                        ccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCAc                                                                                                          >  7:312748/1‑149 (MQ=255)                                          tGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACct                                                                                                        >  2:362033/1‑149 (MQ=255)                                             gcAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACctgct                                                                                                     <  4:182581/149‑1 (MQ=255)                                                   cAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGc                                                                                               >  3:294933/1‑149 (MQ=255)                                                     gCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTg                                                                                             <  5:170881/149‑1 (MQ=255)                                                          tGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                                <  6:313550/93‑1 (MQ=255)                                                          tGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                                >  5:313550/1‑93 (MQ=255)                                                              ctgctgCCACGTTACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAATAGGCTCGCAATCCAGCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCag                                                                                     <  5:245081/148‑1 (MQ=255)                                                              ctgctgCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCag                                                                                     <  8:312748/148‑1 (MQ=255)                                                                    ccACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                                >  2:150370/1‑147 (MQ=255)                                                                    ccACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                                <  1:150370/147‑1 (MQ=255)                                                                    ccACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                               >  7:326829/1‑148 (MQ=255)                                                                       cGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCg                                                                           <  5:283954/149‑1 (MQ=255)                                                                           aCGCTCGACCGAGCCTCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCaa                                                                        <  1:127572/149‑1 (MQ=255)                                                                           aCGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCaa                                                                        <  1:191288/148‑1 (MQ=255)                                                                                 gACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCt                                                                 >  1:186411/1‑149 (MQ=255)                                                                                   cc‑tgctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCg                                                               <  3:205015/149‑1 (MQ=255)                                                                                        ctccctgcttAGTTCACGGAGTGTCACTGTTAATGGCACGCAATCCACGTGGCATTGCCCGTTGTCGCTGTTGAGGAGTGGCAGCCCGGTAAGAAAGCTCACCTGCTTCGCGCTGCAGTCCTGTAGTTGCTCGCAAGCGCGATGGAAc                                                            <  8:326829/148‑1 (MQ=255)                                                                                          ccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTa                                                         >  5:350196/1‑149 (MQ=255)                                                                                          ccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTa                                                         >  7:17476/1‑149 (MQ=255)                                                                                               ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGc                                                      <  5:228998/147‑1 (MQ=255)                                                                                               ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCgg                                                    >  7:36518/1‑149 (MQ=255)                                                                                                      cACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCTCAAGCGCGCTCGAACGTATGCGGCTTGcc                                              <  7:124680/148‑1 (MQ=255)                                                                                                      cACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGccc                                             >  8:229093/1‑149 (MQ=255)                                                                                                       aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCa                                            <  2:186411/149‑1 (MQ=255)                                                                                                       aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGAAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCCCTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTTcc                                              >  4:292983/1‑147 (MQ=255)                                                                                                        cGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAg                                           >  8:26587/1‑149 (MQ=255)                                                                                                         gCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGc                                          <  7:229093/149‑1 (MQ=255)                                                                                                               tCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTg                                     <  6:346927/148‑1 (MQ=255)                                                                                                               tCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGc                                    <  1:362033/149‑1 (MQ=255)                                                                                                                     tAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                                                                            <  2:80142/103‑1 (MQ=255)                                                                                                                     tAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                                                                            >  1:80142/1‑103 (MQ=255)                                                                                                                     tAAAAGGCACGCAATCCACCTGGCAATGCCCGGTGGCGCTGTTCATCAGTGCCAGCCCGGAAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTAGAACGTATGCGGCTAGCACAGCACCTGCTCGcc                               >  4:113490/1‑148 (MQ=255)                                                                                                                        aaGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCa                                         >  6:194733/1‑135 (MQ=255)                                                                                                                        aaGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCa                                         <  5:194733/135‑1 (MQ=255)                                                                                                                                 aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAACGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAGCGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAAcc                  >  6:250987/1‑149 (MQ=255)                                                                                                                                  cgccACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt                 >  8:219891/3‑149 (MQ=255)                                                                                                                                     cACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGt              >  6:180863/1‑149 (MQ=255)                                                                                                                                        cTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGt              <  1:134891/146‑1 (MQ=255)                                                                                                                                              ttGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTa     <  4:86099/149‑1 (MQ=255)                                                                                                                                               tGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTAt    <  7:219891/149‑1 (MQ=255)                                                                                                                                                 cccGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCACCCTGGTCCGAGCCTATGa  >  4:240629/1‑149 (MQ=255)                                                                                                                                                 cccGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGa  <  8:17476/149‑1 (MQ=255)                                                                                                                                                 cccGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAACACCGCAACCTGGTCCGAGCCTATGa  <  6:350196/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACC‑TGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGA  >  NC_002947/2154251‑2154542

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