Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,799,740:1	+A	9.9%	coding (75/1125 nt)	PP_3360 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,799,740	1	.	A	9.9%	262.1 / 24.4	71	coding (75/1125 nt)	PP_3360	membrane protein
Reads supporting (aligned to +/- strand):  ref base . (32/32);  new base A (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 = 5.07e-01

CGATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGT  >  NC_002947/3799594‑3799885                                                                                                                                                    |                                                                                                                                                    cGATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGgcg‑c                                                                                                                                                    <  4:70704/148‑1 (MQ=255) cGATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGgcg‑c                                                                                                                                                    <  4:212736/148‑1 (MQ=255)   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                                                                          >  5:116609/1‑149 (MQ=255)                                                                         aTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                          >  2:24473/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                                       >  6:180791/1‑149 (MQ=255)                                                                                                            caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                       >  2:215816/1‑149 (MQ=255)                                                                                                            caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                       >  5:76479/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)                                                                                                                                                    |                                                                                                                                                    CGATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGT  >  NC_002947/3799594‑3799885

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