Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,799,748	Δ1 bp	10.0%	coding (83/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,748	0	T	.	10.0%	199.9 / 16.9	70	coding (83/1125 nt)	PP_3360	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base . (3/4);  total (34/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 = 1.00e+00

GTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGC  >  NC_002947/3799604‑3799892                                                                                                                                                   |                                                                                                                                                  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                                                                                                                                          >  5:64441/1‑133 (MQ=255)                    ccAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACT‑GTCAG‑CGCCAccag                                                                                                                                          <  6:64441/135‑3 (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                                                                                                             <  8:30602/138‑1 (MQ=255)                                              ccAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCa                                                                                                             >  7:30602/1‑138 (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                                              >  6:180791/1‑149 (MQ=255)                                                                                                  caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                              >  5:76479/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                                <  1:361/128‑1 (MQ=255)                                                                                                                                    gtggCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTc                                >  2:361/1‑128 (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)                                                                                                                                                   |                                                                                                                                                  GTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CT‑GACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAA‑GGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGC  >  NC_002947/3799604‑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: ‑