Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,799,744	A→T	9.9%	T27S (ACA→TCA)	PP_3360 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,799,744	0	A	T	9.9%	195.8 / 10.7	71	T27S (ACA→TCA)	PP_3360	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base T (3/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

ATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGC  >  NC_002947/3799596‑3799892                                                                                                                                                       |                                                                                                                                                      aTCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑Ga                                                                                                                                                       <  8:327671/149‑1 (MQ=255)    gACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACa                                                                                                                                                     <  8:171501/148‑1 (MQ=255)    gACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAg                                                                                                                                                    >  8:49931/1‑149 (MQ=255)         gTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGc                                                                                                                                                >  3:197553/1‑148 (MQ=255)            gcgcCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCAcc                                                                                                                                            >  1:73576/1‑149 (MQ=255)             cgcCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAcca                                                                                                                                           <  6:249636/149‑2 (MQ=255)                 tctTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccag                                                                                                                                          >  1:287137/1‑144 (MQ=255)                   tttGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCTCCACCTTGGTCt                                                                                                                                     >  1:302481/1‑149 (MQ=255)                     tGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccagcaga                                                                                                                                      <  2:287137/146‑7 (MQ=255)                       aCGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCt                                                                                                                                 >  1:174337/1‑149 (MQ=255)                         gTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTgg                                                                                                                               <  2:302481/149‑1 (MQ=255)                           cccAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTc                                                                                                                             >  8:237412/1‑149 (MQ=255)                            ccAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccag                                                                                                                                          <  6:64441/135‑3 (MQ=255)                            ccAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccag                                                                                                                                          >  5:64441/1‑133 (MQ=255)                                    gCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTgcc                                                                                                                    >  5:62826/1‑149 (MQ=255)                                     cACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTgccg                                                                                                                   <  2:174337/149‑1 (MQ=255)                                            gCAACTCAACCCAGAGGCGAGCTACTCATGACCGACCGCATCATCGAGTTTTTCAACGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTTGTCTTCCTGGTCGCCCGTGCCGCCATCa                                                                                                             <  8:234736/148‑1 (MQ=255)                                                aTCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTc                                                                                                        >  5:67022/1‑149 (MQ=255)                                                      ccAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCa                                                                                                             >  7:30602/1‑138 (MQ=255)                                                      ccAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCa                                                                                                             <  8:30602/138‑1 (MQ=255)                                                         gagGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATgcgc                                                                                               <  4:344356/149‑1 (MQ=255)                                                           ggCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCg                                                                                             <  7:237412/149‑1 (MQ=255)                                                               aGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTa                                                                                         <  2:73576/149‑1 (MQ=255)                                                                 gTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACg                                                                                       >  5:150230/1‑149 (MQ=255)                                                                     tCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTATGgcgc                                                                                   >  6:192065/1‑149 (MQ=255)                                                                     tCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGgcgc                                                                                   >  4:264554/1‑149 (MQ=255)                                                                       aTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                                 >  6:214140/1‑149 (MQ=255)                                                                       aTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                                 >  2:24473/1‑149 (MQ=255)                                                                       aTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                                 >  5:116609/1‑149 (MQ=255)                                                                         gAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGt                                                                               >  3:158137/1‑149 (MQ=255)                                                                          aaCGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTc                                                                              <  6:62826/149‑1 (MQ=255)                                                                           aCGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTcc                                                                             <  1:24473/149‑1 (MQ=255)                                                                              aCCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGaa                                                                          <  4:101114/149‑1 (MQ=255)                                                                                cGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAaca                                                                        <  4:197553/149‑1 (MQ=255)                                                                                 gCATCATCGCGTTTTTCACCGAACCCACCGTATTCGTCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACTGCTCTGGTCCGAacac                                                                       <  3:286758/149‑1 (MQ=255)                                                                                       tcGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGt                                                                                          >  2:350976/1‑124 (MQ=255)                                                                                       tcGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGc                                                                 >  6:169282/1‑149 (MQ=255)                                                                                         gcgTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGcc                                                                <  4:301284/148‑1 (MQ=255)                                                                                          cgTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCa                                                              <  5:214140/149‑1 (MQ=255)                                                                                             ttttCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACt                                                            >  2:133712/1‑148 (MQ=255)                                                                                              tttCACCGAAACCACCGTATTCGGCATCTCGATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccag                                                                                                                                          <  7:26205/69‑3 (MQ=25)                                                                                              tttCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccag                                                                                                                                          >  8:26205/1‑67 (MQ=25)                                                                                                tCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGc                                                        >  4:316912/1‑149 (MQ=255)                                                                                                 cACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTATGGCGCTGGTCCGAACACGAAGGCCCACTGAGc                                                        <  5:192065/148‑1 (MQ=255)                                                                                                 cACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACTAAGGCCCACTGAGc                                                        <  1:133712/148‑1 (MQ=255)                                                                                                   ccGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGc                                                        >  7:83037/1‑146 (MQ=255)                                                                                                       aaCCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATaa                                                 >  7:141548/1‑149 (MQ=255)                                                                                                        accaccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAAt                                                <  8:83037/149‑1 (MQ=255)                                                                                                         ccaccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATg                                               <  5:169282/149‑1 (MQ=255)                                                                                                         ccaccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATg                                               <  6:344851/149‑1 (MQ=255)                                                                                                          caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                              >  5:76479/1‑149 (MQ=255)                                                                                                          caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                              >  6:180791/1‑149 (MQ=255)                                                                                                          caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                              >  2:215816/1‑149 (MQ=255)                                                                                                           accGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGc                                             >  6:216796/1‑149 (MQ=255)                                                                                                               tATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑g                                        <  5:216796/149‑1 (MQ=255)                                                                                                               tATTCGGCATCACCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑g                                        <  6:116609/149‑1 (MQ=255)                                                                                                                   cGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTg                                     <  6:76479/148‑1 (MQ=255)                                                                                                                    ggCATCTCCATTCTCAATCTGCTGGTGGCG‑CTAGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGc                                    >  1:42664/1‑149 (MQ=255)                                                                                                                           ccATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCgg                             >  7:330016/1‑148 (MQ=255)                                                                                                                            cATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTcc                               <  3:264554/145‑1 (MQ=255)                                                                                                                            cATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCg                              <  3:316912/146‑1 (MQ=255)                                                                                                                            cATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCa                           >  5:23495/1‑149 (MQ=255)                                                                                                                            cATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAG‑CAGATAATGGCCAA‑GGTGCTTTCCg                              <  6:67022/145‑1 (MQ=255)                                                                                                                                     tCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACt                  >  7:146820/1‑149 (MQ=255)                                                                                                                                         ctgGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCt              <  1:215816/149‑1 (MQ=255)                                                                                                                                           ggtggCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCtg             >  2:276808/1‑148 (MQ=255)                                                                                                                                           ggtggCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCtg             <  1:276808/148‑1 (MQ=255)                                                                                                                                           gggggCG‑CT‑GACAGGCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCtgt            <  8:74033/146‑1 (MQ=255)                                                                                                                                            gtggCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTc                                >  2:361/1‑128 (MQ=255)                                                                                                                                            gtggCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTc                                <  1:361/128‑1 (MQ=255)                                                                                                                                            ggggCG‑CG‑GACAGGCGCCACCTTGGTCTTCCTGGTCGGCCGTGCCGCCATCACTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGCAGGCCCACCGAGCCAGATAATGGCCAACGGTGCTTTCCGGCACCAGCAACTTTCt              <  1:197629/145‑1 (MQ=255)                                                                                                                                               gcg‑cT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCtgttgt         <  8:141548/148‑1 (MQ=255)                                                                                                                                                       aCAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGc  <  8:146820/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGC  >  NC_002947/3799596‑3799892

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