Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,176,459	T→A	5.7%	V587E (GTG→GAG)	PP_2787 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,176,459	0	T	A	5.7%	221.2 / 2.5	70	V587E (GTG→GAG)	PP_2787	transporter
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (2/2);  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 = 5.71e-01

ACCACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGAT  >  NC_002947/3176317‑3176607                                                                                                                                               |                                                                                                                                                     aCCACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGc                                                                                                                                                 <  6:332707/148‑1 (MQ=255)  ccACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAg                                                                                                                                               <  5:248845/149‑1 (MQ=255)   cACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAgg                                                                                                                                              >  1:34782/1‑149 (MQ=255)   cACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGCACCCGCAg                                                                                                                                               >  8:336702/1‑148 (MQ=255)   cACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAAAGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAgg                                                                                                                                              >  1:286633/1‑149 (MQ=255)        aaaCCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGt                                                                                                                                          <  2:34782/148‑1 (MQ=255)        aaaCCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGt                                                                                                                                          <  8:321062/148‑1 (MQ=255)         aaCCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTgg                                                                                                                                        <  7:325351/149‑1 (MQ=255)         aaCCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTgg                                                                                                                                        <  2:187020/149‑1 (MQ=255)         aaCCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGCACCCGCAggaaagca                                                                                                                                        <  7:336702/149‑7 (MQ=255)          aCCATGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGa                                                                                                                                       <  4:332420/149‑1 (MQ=255)                     ccTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACgg                                                                                                                            >  6:310942/1‑149 (MQ=255)                     ccTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACg                                                                                                                             <  1:84233/148‑1 (MQ=255)                                tGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCgatgat                                                                                                                 >  4:57722/1‑149 (MQ=255)                                  gACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCgatgat                                                                                                                 <  1:63105/147‑1 (MQ=255)                                   aCCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGAc                                                                                                              >  3:181333/1‑149 (MQ=255)                                     cAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACACTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACgg                                                                                                            >  5:284875/1‑149 (MQ=255)                                        ttCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGCACCCGCAggaaa                                                                                                                                           >  2:299889/1‑112 (MQ=255)                                        ttCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGCACCCGCAggaaa                                                                                                                                           <  1:299889/115‑4 (MQ=255)                                           atcaAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCgatga                                                                                                                  >  3:8856/1‑137 (MQ=255)                                           atcaAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCgatga                                                                                                                  <  4:8856/137‑1 (MQ=255)                                           atcaAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAAcc                                                                                                      >  2:203431/1‑149 (MQ=255)                                             caAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCg                                                                                                     >  5:234699/1‑148 (MQ=255)                                                   ccAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCt                                                                                              >  6:126279/1‑149 (MQ=255)                                                   ccAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCt                                                                                              <  7:222449/149‑1 (MQ=255)                                                     aGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTtc                                                                                            <  1:203431/149‑1 (MQ=255)                                                       ccGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTtcga                                                                                          >  4:141753/1‑149 (MQ=255)                                                              tGAGTTACACCTATTCAGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTtcgatgtcgat                                                                                   >  5:78809/1‑149 (MQ=255)                                                                    acacCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAAc                                                                             >  6:203439/1‑149 (MQ=255)                                                                      acCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGTAACTc                                                                           >  7:106775/1‑149 (MQ=255)                                                                        cTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGga                                                                                                   >  5:228905/1‑123 (MQ=255)                                                                        cTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGga                                                                                                   <  6:228905/123‑1 (MQ=255)                                                                        cTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGTAACTCGc                                                                         <  8:106775/149‑1 (MQ=255)                                                                        cTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGc                                                                         <  5:203439/149‑1 (MQ=255)                                                                        cTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGc                                                                         <  3:141753/149‑1 (MQ=255)                                                                           ttCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGa                                                                                                                                       <  4:136929/84‑1 (MQ=255)                                                                           ttCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGa                                                                                                                                       >  3:136929/1‑84 (MQ=255)                                                                           ttCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTtcgatgt                                                                                       >  1:143663/1‑132 (MQ=255)                                                                           ttCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTtcgatgt                                                                                       <  2:143663/132‑1 (MQ=255)                                                                              ggTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATc                                                                    <  4:126565/148‑1 (MQ=255)                                                                                tGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGaa                                                                 >  4:170808/1‑149 (MQ=255)                                                                                    gACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGaacaac                                                             <  6:78809/149‑1 (MQ=255)                                                                                       cTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAAcaca                                                          <  5:310942/149‑1 (MQ=255)                                                                                         aCATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGCAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACAcat                                                         >  4:62534/1‑148 (MQ=255)                                                                                            tCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCatg                                                     <  3:170808/149‑1 (MQ=255)                                                                                              aaGCGCATGAACCTCGGGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCatgat                                                   >  6:344332/1‑149 (MQ=255)                                                                                               aGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCatgatg                                                  >  3:4265/1‑149 (MQ=255)                                                                                                  gcATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTc                                               <  2:286633/149‑1 (MQ=255)                                                                                                     tGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCgag                                            <  6:234699/149‑1 (MQ=255)                                                                                                              tGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTggc                                   <  1:148983/149‑1 (MQ=255)                                                                                                              tGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTggc                                   <  4:181333/149‑1 (MQ=255)                                                                                                                 tCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTggcgg                                 >  1:352024/1‑148 (MQ=255)                                                                                                                         aTGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACtt                        <  5:344332/149‑1 (MQ=255)                                                                                                                          tGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTg                       <  3:57722/149‑1 (MQ=255)                                                                                                                             aCCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAACCCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCatgat                                                   >  8:151873/1‑118 (MQ=255)                                                                                                                               ccTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATgg                                                                                 <  5:271402/86‑1 (MQ=255)                                                                                                                               ccTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATgg                                                                                 >  6:271402/1‑86 (MQ=255)                                                                                                                               ccTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGaa                  <  4:4265/149‑1 (MQ=255)                                                                                                                                 tGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCa                     <  8:17414/144‑1 (MQ=255)                                                                                                                                 tGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCa                     >  7:17414/1‑144 (MQ=255)                                                                                                                                  gCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGt               >  8:16089/1‑149 (MQ=255)                                                                                                                                    tttCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGc                      <  4:88237/140‑1 (MQ=255)                                                                                                                                    tttCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGc                      >  3:88237/1‑140 (MQ=255)                                                                                                                                    tttCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGc             >  3:59787/1‑149 (MQ=255)                                                                                                                                        cTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCa                     >  7:134731/1‑137 (MQ=255)                                                                                                                                        cTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCa                     <  8:134731/137‑1 (MQ=255)                                                                                                                                        cTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCaa         >  6:158989/1‑149 (MQ=255)                                                                                                                                        cTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCaa         >  3:338907/1‑149 (MQ=255)                                                                                                                                         tGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAAt        >  2:14786/1‑149 (MQ=255)                                                                                                                                            gggTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGcc     <  5:126279/149‑1 (MQ=255)                                                                                                                                               tGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGAt  <  7:16089/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     ACCACTGAAACCGTGCAGTTCCTCGACCGTCTGGACCAGTTCATCAAGGCCCAGCCGAACGTGAGTTACACCTATTCGGTGGCGACCTACATCAAGCGCATGAACCTCGTGCTCAACGACATGAACCCTGCTTTCCTGCGGGTGCCGCAGGCGGTGGAGCAGGTGACGGTGGTCGATGATGACGGCAAACCGGAGACCTTCGATGTCGATGGCAACTCGCTCATCGAACAACACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGAT  >  NC_002947/3176317‑3176607

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