Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,760,260	T→G	9.4%	F93V (TTC→GTC)	pvdM →	dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,760,260	0	T	G	9.4%	158.8 / 10.4	64	F93V (TTC→GTC)	pvdM	dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base G (3/3);  total (32/32)
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.24e-01

CCACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTGCGC  >  NC_002947/4760112‑4760406                                                                                                                                                     |                                                                                                                                                   ccACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGtt                                                                                                                                                   >  2:578189/1‑150 (MQ=255) ccACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGtt                                                                                                                                                   >  2:1062990/1‑150 (MQ=255)  cACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTc                                                                                                                                                  >  2:39881/1‑150 (MQ=255)     ttCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCACTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGAc                                                                                                                                               >  1:651334/1‑150 (MQ=255)     ttCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGaa                                                                                                                                               >  1:596660/1‑149 (MQ=255)      tCGACAAGGACGTCCCCCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGaac                                                                                                                                              <  2:596660/150‑3 (MQ=255)         aCAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTgg                                                                                                                                           <  1:529882/150‑1 (MQ=255)           aaGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGAACTGGcc                                                                                                                                         <  1:39881/150‑1 (MQ=255)               aCGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGCCCTGGCCAAAg                                                                                                                                     <  1:1062990/150‑1 (MQ=255)                       cccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgcgc                                                                                                                             >  1:1274083/1‑150 (MQ=255)                       cccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgcgc                                                                                                                             >  1:901631/1‑150 (MQ=255)                        ccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgcgcg                                                                                                                            >  2:1018301/1‑150 (MQ=255)                         cATAGAGGTCGGGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCCGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCgg                                                                                                                           <  1:641180/150‑1 (MQ=255)                          aTAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCggg                                                                                                                          >  1:672867/1‑150 (MQ=255)                             gagGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGa                                                                                                                       <  2:755542/150‑1 (MQ=255)                                 tCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGCAAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCCGACAAGGATGGCCCGAGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGt                                                                                                                   >  2:500988/1‑150 (MQ=255)                                  cGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTc                                                                                                                  <  2:1199742/150‑1 (MQ=255)                                            gCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGcc                                                                                                        <  2:1235536/150‑1 (MQ=255)                                                aTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                    <  2:1194732/150‑1 (MQ=255)                                                 tGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGAc                                                                                                   >  2:488430/1‑150 (MQ=255)                                                       gCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGAc                                                                                                                                               >  2:482432/1‑100 (MQ=255)                                                       gCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGAc                                                                                                                                               <  1:482432/100‑1 (MQ=255)                                                          agagCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCAAAAGCCACGAAGCTGGGGATGAGATCGGGGGCGAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGc                                                                                          >  1:1120152/1‑150 (MQ=255)                                                                gATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGcc                                                                                    >  1:633370/1‑150 (MQ=255)                                                                gATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGcc                                                                                    >  1:638921/1‑150 (MQ=255)                                                                 aTCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGccc                                                                                   >  2:1076811/1‑150 (MQ=255)                                                                 aTCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGccc                                                                                   <  2:321859/150‑1 (MQ=255)                                                                 aTCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGccc                                                                                   <  1:1100531/150‑1 (MQ=255)                                                                   cATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGa                                                                                 <  1:578189/150‑1 (MQ=255)                                                                   cATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGa                                                                                 <  2:78285/150‑1 (MQ=255)                                                                     tttCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGt                                                                               >  2:175188/1‑150 (MQ=255)                                                                        ggttcgACAGCTATCTAACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCt                                                                            <  2:459122/149‑1 (MQ=255)                                                                         gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTg                                                                           <  2:672867/150‑1 (MQ=255)                                                                         gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGACGACCGTCACTGGCTGCCCTGACCGTCTTCGCTTGGCCCGGGTCCTg                                                                           >  2:103362/1‑150 (MQ=255)                                                                         gttcgACAGCCATCGCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTg                                                                           <  2:519231/150‑1 (MQ=255)                                                                                    aTCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGCTGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACa                                                                >  1:172450/1‑150 (MQ=255)                                                                                       tcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACg                                                             <  2:901631/150‑1 (MQ=255)                                                                                        cACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACgg                                                            >  1:64908/1‑150 (MQ=255)                                                                                         cgccggCCGCGGGGCTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACGGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGc                                                           <  2:172450/145‑1 (MQ=255)                                                                                          cccTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGcc                                                          <  2:751060/150‑1 (MQ=255)                                                                                                  cTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGa                                                                                 <  2:541371/119‑1 (MQ=255)                                                                                                  cTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGa                                                                                 >  1:541371/1‑119 (MQ=255)                                                                                                  cTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACg                                                             >  2:56335/1‑139 (MQ=255)                                                                                                  cTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACg                                                             <  1:56335/139‑1 (MQ=255)                                                                                                   tGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTAAATGTCGTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCaccacc                                                 <  1:997168/150‑1 (MQ=255)                                                                                                   tGGACTTCGGCACCGAAGGCCGCGAAGCTAACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCaccacc                                                 <  1:325571/150‑1 (MQ=255)                                                                                                        ttCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCAccaccccacc                                            >  1:246972/1‑150 (MQ=255)                                                                                                        ttCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCAccaccccacc                                            >  2:692491/1‑150 (MQ=255)                                                                                                        ttCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCAccaccccacc                                            >  2:692479/1‑150 (MQ=255)                                                                                                        ttCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCAccaccccacc                                            >  1:1065479/1‑150 (MQ=255)                                                                                                                cGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGttt                                    >  2:822325/1‑150 (MQ=255)                                                                                                                   aGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTtcgt                                 <  2:1065479/150‑1 (MQ=255)                                                                                                                     gCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGt                                                                                                >  2:313598/1‑85 (MQ=255)                                                                                                                     gCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGt                                                                                                <  1:313598/85‑1 (MQ=255)                                                                                                                       cgcgAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTTGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGAt                             >  1:743690/1‑150 (MQ=255)                                                                                                                       cgcgAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAGGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGAt                             >  1:1148423/1‑150 (MQ=255)                                                                                                                             gCTGACCAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCggg                       <  1:1018301/150‑1 (MQ=255)                                                                                                                                gACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGa                    >  2:199505/1‑150 (MQ=255)                                                                                                                                 aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTtcgtc                                >  1:1223152/1‑137 (MQ=255)                                                                                                                                 aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTtcgtc                                <  2:1223152/137‑1 (MQ=255)                                                                                                                                 aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGAc                   >  2:1001585/1‑150 (MQ=255)                                                                                                                                           cccGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGa         <  2:1120152/150‑1 (MQ=255)                                                                                                                                            tttggCCAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGc                           >  1:640173/4‑131 (MQ=255)                                                                                                                                            tttggCCAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGc                           <  2:640173/128‑1 (MQ=255)                                                                                                                                                 ccAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGATTTTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTgcg   <  2:743690/150‑1 (MQ=255)                                                                                                                                                  cAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTgcgc  <  1:1076811/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTGCGC  >  NC_002947/4760112‑4760406

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