Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,962,593	C→A	7.6%	L286M (CTG→ATG)	PP_1760 →	ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,962,593	0	C	A	7.6%	179.0 / 2.1	57	L286M (CTG→ATG)	PP_1760	ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base A (3/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGAAACCCTCGATGCCATCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCTTAATCTGCA  >  NC_002947/1962446‑1962738                                                                                                                                                         |                                                                                                                                                   gCGGAAACCCTCGATGCCATCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGcc                                                                                                                                                    <  7:157593/148‑1 (MQ=255)    gAAACCCTCGATGCCATCCAGGAAATCCGCGCCGGCAACCGCAAGGGCTTCTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTTCGCGGTGGACTCGCAATGGAAAATCGATGCCAGTGGCCGTGCCAGTGGCct‑gct                                                                                                                                               <  3:231266/149‑1 (MQ=255)     aaaCCCTCGATGCCATCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTGCACTCCGCCTGGGCCTCCCCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGcc                                                                                                                                                    <  1:111549/149‑1 (MQ=255)              aTGCCATCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAACAGC‑AGGCCACTGGCACGGCCACTGGCATCGCTTTTCCattgc                                                                                                                                     <  4:148210/149‑5 (MQ=25)                  cATCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCa                                                                                                                                 <  3:292919/149‑1 (MQ=255)                   aTCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTc                                                                                                                                     <  1:18369/144‑1 (MQ=255)                   aTCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTc                                                                                                                                     >  2:18369/1‑144 (MQ=255)                       aGGCAAACCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTCGACTCGCAAGGGAAACGCGATCCCAGTGGCCGTCCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGAc                                                                                                                            <  8:227641/149‑1 (MQ=255)                       aGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCa                                                                                                                                 >  5:219875/1‑144 (MQ=255)                       aGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCa                                                                                                                                 <  6:219875/144‑1 (MQ=255)                       aGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGAc                                                                                                                            <  3:277736/149‑1 (MQ=255)                          aaaTCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCgg                                                                                                                                   >  3:216434/1‑139 (MQ=255)                          aaaTCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCgg                                                                                                                                   <  4:216434/139‑1 (MQ=255)                                gcgcCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGt                                                                                                                                             <  8:305070/123‑1 (MQ=255)                                gcgcCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGt                                                                                                                                             >  7:305070/1‑123 (MQ=255)                                        aaCCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAACAGC‑AGGCCACTGGCACGGCCACTGGCATCGCTTTTCCa                                                                                                                                         >  6:290517/1‑119 (MQ=21)                                        aaCCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAACAGC‑AGGCCACTGGCACGGCCACTGGCATCGCTTTTCCa                                                                                                                                         <  5:290517/119‑1 (MQ=21)                                            gCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGCCTCGCAAGGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCAG‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATg                                                                                                       >  5:279148/1‑149 (MQ=255)                                               aGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTg                                                                                                    >  4:285589/1‑149 (MQ=255)                                                   cTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCg                                                                                                >  1:77035/1‑149 (MQ=255)                                                   cTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCg                                                                                                >  3:151090/1‑149 (MQ=255)                                                       ttCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCCCGAGGTGCGAGATTACACCGTGGACTCGCAATGGAAAACCGATGCCAGTCGCCGTGACAGTGGCAT‑GCTGTTCCAGATCGGCATCGCCATCTTCCGCGCGGCGGCCATGCTGACCGTGCt                                                                                            >  7:225461/1‑149 (MQ=255)                                                        tCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTg                                                                                           >  1:62135/1‑149 (MQ=255)                                                         ccTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTcc                                                                                                                     >  4:10116/1‑122 (MQ=255)                                                         ccTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTcc                                                                                                                     <  3:10116/122‑1 (MQ=255)                                                                   cgccTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGa                                                                                   <  4:151090/149‑1 (MQ=255)                                                                    gccTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGAc                                                                                  >  5:188508/1‑149 (MQ=255)                                                                     ccTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTc                                                                                      >  5:35749/1‑144 (MQ=255)                                                                     ccTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTc                                                                                      <  6:35749/144‑1 (MQ=255)                                                                      cTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTg                                                                                              >  1:130477/1‑135 (MQ=255)                                                                      cTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTg                                                                                              <  2:130477/135‑1 (MQ=255)                                                                                gCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCg                                                                        >  6:75506/1‑149 (MQ=255)                                                                                       cgcgAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCGAGTGGCCGTGCCAGTGGCCT‑GGTGTTCCAGTTCGGCATCGAGATCTTCCGCGCGGCGGGCATGCTGACCGTGCTGTTTTCGGACCTGTCGATCGGCCAGAt                                                                 <  1:319200/149‑1 (MQ=255)                                                                                       cgcgAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGAACGGCCAGAt                                                                 >  3:68113/1‑149 (MQ=255)                                                                                               gcgcGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCgg                                                         >  4:263487/1‑149 (MQ=255)                                                                                               gcgcGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCgg                                                         >  4:236015/1‑149 (MQ=255)                                                                                                cgcgATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCgg                                                         <  6:188508/148‑1 (MQ=255)                                                                                                     ttACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCa                                                   <  5:75506/149‑1 (MQ=255)                                                                                                          ccGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTAc                                              <  2:77035/149‑1 (MQ=255)                                                                                                                aCTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTgg                                        <  3:236015/149‑1 (MQ=255)                                                                                                                  tCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGt                                       >  3:298834/1‑148 (MQ=255)                                                                                                                   cGCCATGGAAACGCGATGCCAGGGGCCGTGCCAGGGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTc                                     <  6:279148/149‑1 (MQ=255)                                                                                                                    gCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCa                                    >  4:8698/1‑149 (MQ=255)                                                                                                                    gCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCa                                    <  2:62135/149‑1 (MQ=255)                                                                                                                     cAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCAt                                   <  1:212305/149‑1 (MQ=255)                                                                                                                        tGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTAc                                              >  6:143547/1‑135 (MQ=255)                                                                                                                        tGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTAc                                              <  5:143547/135‑1 (MQ=255)                                                                                                                          gAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATc                               <  3:285589/148‑1 (MQ=255)                                                                                                                                  aTGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTg                      <  4:298834/149‑1 (MQ=255)                                                                                                                                    gCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGt                       >  7:205562/1‑146 (MQ=255)                                                                                                                                    gCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGt                       <  8:205562/146‑1 (MQ=255)                                                                                                                                     ccAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGa                    <  1:91275/148‑1 (MQ=255)                                                                                                                                     ccAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAg                   <  3:263487/149‑1 (MQ=255)                                                                                                                                       aGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGa                                 >  5:323621/1‑133 (MQ=255)                                                                                                                                       aGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGa                                 <  6:323621/133‑1 (MQ=255)                                                                                                                                             cGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCt           <  3:8698/149‑1 (MQ=255)                                                                                                                                             cGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCt           >  2:12210/1‑149 (MQ=255)                                                                                                                                                   aGTGGCCT‑GCTGTTCCAGTTAGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCTTAATCt     <  4:238270/149‑1 (MQ=255)                                                                                                                                                     tGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCTTAATCTg    >  4:61476/1‑148 (MQ=255)                                                                                                                                                       gCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCTTAATCTGCa  <  4:68113/148‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                   GCGGAAACCCTCGATGCCATCCAGGAAATCCGCGCCGGCAACCGCCAGGGCTACTTCCTTG‑‑‑GCCGCCTGGGCTT‑‑GCGTGCCCGCGAGGTGCGCGATTACGCCGTGGACTCGCAATGGAAAAGCGATGCCAGTGGCCGTGCCAGTGGCCT‑GCTGTTCCAGTTCGGCATCGACATCTTCCGCGCGGCGGCCATGCTGACCGTGCTGTTTTCCGACCTGTCGATCGGCCAGATGCTGGCGGTGTTCAGCTACCTGTGGTTCATGATCGGTCCGGTGGAGCAACTGCTTAATCTGCA  >  NC_002947/1962446‑1962738

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