Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,877,159	C→A	8.7%	A83D (GCC→GAC)	PP_3425 →	RND family transporter MFP subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,877,159	0	C	A	8.7%	205.2 / 6.9	66	A83D (GCC→GAC)	PP_3425	RND family transporter MFP subunit
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base A (3/3);  total (32/34)
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

GCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTG  >  NC_002947/3877011‑3877301                                                                                                                                                      |                                                                                                                                               gCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCCCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGgcgc                                                                                                                                                >  6:251691/1‑149 (MQ=255)      gcCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCGGTGAAGGCGACCTgg                                                                                                                                           <  1:241354/149‑1 (MQ=255)      gcCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAgg                                                                                                                                           <  4:62333/149‑1 (MQ=255)          aaGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCGGTGAAGGCGACCTGGTcga                                                                                                                                       <  7:91374/149‑3 (MQ=255)           aGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTc                                                                                                                                         >  6:230709/1‑146 (MQ=255)           aGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCATTCACCGAAGGCGCCCAGGTc                                                                                                                                         >  6:339045/1‑146 (MQ=255)              tAAGTGTGGGCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGTCCGCCTCGAAGCCCCGTAAACCGTCGAAGTGCGCCCGCGTGTCTCCG‑GCCAGATCGACCTGGTCGCCTTCACCGAAGGCGCCGAGGTCAAGaaa                                                                                                                                   <  7:189875/149‑1 (MQ=255)              tAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGaa                                                                                                                                    <  2:195187/148‑1 (MQ=255)              tAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGaa                                                                                                                                    >  1:195187/1‑148 (MQ=255)                aGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAgg                                                                                                                                 <  4:157814/149‑1 (MQ=255)                 gtgtGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGc                                                                                                                                <  1:308996/149‑1 (MQ=255)                        cAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGCGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACctg                                                                                                                          <  1:102296/149‑1 (MQ=255)                               aTCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTcc                                                                                                                   <  8:191352/148‑1 (MQ=255)                                tCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAg                                                                                                                 <  5:230709/149‑1 (MQ=255)                                tCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAg                                                                                                                 <  2:69965/149‑1 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAg                                                                                                                                  <  2:88508/129‑1 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAg                                                                                                                                  >  1:88508/1‑129 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGt                                                                                                                      >  1:124615/1‑141 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGt                                                                                                                      <  2:124615/141‑1 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTc                                                                                                                    >  6:284086/1‑143 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTc                                                                                                                    >  8:327018/1‑143 (MQ=255)                                   aacaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTCCCCCCGCGGGTCTCCC‑GCCAGATCGACCACGTCGCCGTCCACGAAGGCGCCCAGCTCAAGCAAGGCGACCTCCTCATCCCGAt                                                                                                               >  4:215587/1‑148 (MQ=255)                                        ccGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGAccc                                                                                                         <  6:220175/149‑1 (MQ=255)                                        ccGATCACCGAATGGGACGAGTTCACCGGCCACCTCGAAGCCCCGGAAACCGTCGAAGTGCGCACGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGAccc                                                                                                         >  7:242585/1‑149 (MQ=255)                                         cGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTcc                                                                                                                   >  4:195400/1‑138 (MQ=255)                                         cGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTcc                                                                                                                   <  3:195400/138‑1 (MQ=255)                                              aCCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGccc                                                                                                   >  2:69581/1‑149 (MQ=255)                                              aCCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGccc                                                                                                   <  7:327018/149‑1 (MQ=255)                                                cGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGt                                                                                                 <  1:343860/149‑1 (MQ=255)                                                      ggACGAGTTCACCTGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAgg                                                                                           <  2:198737/149‑1 (MQ=255)                                                         cGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTg                                                                                        <  5:284086/149‑1 (MQ=255)                                                          gcgTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGCCCCGCGCCCGTTCCAGGCTg                                                                                        >  7:12743/3‑148 (MQ=255)                                                           aGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCCCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGaa                                                                                      <  5:251691/149‑1 (MQ=255)                                                                     ccgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTACG‑GCCAGATCGAACAGGTCGACTTCACCGAAGGCGCCCAGGTCAATAAAGGCGACCTGCTGTTCCAGATCGACACGCGCCCGTTACAGGCTGAATTAcgccgc                                                                             >  8:235218/1‑148 (MQ=255)                                                                            gAAGCCCCGGAAACCGTCGAAGTGCTGCCGCGGGTCTCCG‑GCCAGATCGCGCAGGGCGCCTTGTCGGAAGGCGCGCAGGTCATTATTGGCGACATGGTGTTCCAGATCGACCGGCGCGAGTTCGAGGCTGTAGTCCGCCGCCTCGAAGc                                                                     <  4:29703/149‑1 (MQ=255)                                                                            gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCATGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAg                                                                      <  5:339045/148‑1 (MQ=255)                                                                                   cGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTg                                                              >  3:354060/1‑149 (MQ=255)                                                                                     gAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGc                                                             >  7:240010/1‑148 (MQ=255)                                                                                          cGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTg                                                              <  5:183237/142‑1 (MQ=255)                                                                                          cGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTg                                                              >  6:183237/1‑142 (MQ=255)                                                                                          cGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACagg                                                       >  6:217768/1‑149 (MQ=255)                                                                                          cGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACagg                                                       >  3:299621/1‑149 (MQ=255)                                                                                                   gcgcCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGa                                                                        >  3:177056/1‑123 (MQ=255)                                                                                                   gcgcCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGa                                                                        <  4:177056/123‑1 (MQ=255)                                                                                                    cgcCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCAc                                             >  8:213762/1‑149 (MQ=255)                                                                                                       ccGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaaggc                                                >  1:243231/1‑143 (MQ=255)                                                                                                       ccGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaaggc                                                <  2:243231/143‑1 (MQ=255)                                                                                                       ccGAGGGTCTCCG‑GCCAGATCGACCACGTCGCCTTCACCGAAGGCGCCCAGTTCTAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCATTTCCAGGCTGAATTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCg                                           <  1:132383/148‑1 (MQ=255)                                                                                                             gTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaa                                                   <  1:191929/134‑1 (MQ=255)                                                                                                             gTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaa                                                   >  2:191929/1‑134 (MQ=255)                                                                                                                     gCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgcgc                              >  1:90748/1‑148 (MQ=255)                                                                                                                      ctttcttgACCTGGGCGCCTTCGGTGAAGGCGACCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGcc                             >  3:33606/8‑148 (MQ=38)                                                                                                                      ctttcttgACCTGGGCGCCTTCGGTGAAGGCGACCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGcc                             >  3:56749/8‑148 (MQ=38)                                                                                                                        aGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAAc                          <  7:213762/149‑1 (MQ=255)                                                                                                                          aTCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACg                         >  2:64447/1‑148 (MQ=255)                                                                                                                              aCCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGcc                    >  6:121847/1‑149 (MQ=255)                                                                                                                                  ggTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCccgcc                >  3:348786/1‑149 (MQ=255)                                                                                                                                    gcGCCTTCGGTGAAGGCGACCTGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCcgcg                                                                                                      >  2:132477/2‑61 (MQ=21)                                                                                                                                    gcGCCTTCGGTGAAGGCGACCTGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCcgcg                                                                                                      <  1:132477/60‑1 (MQ=21)                                                                                                                                       ccTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGc           <  8:242585/149‑1 (MQ=255)                                                                                                                                       catttaCCTAATGCTCCCATGTTAATAAATGCGACCTGCTTTTTCATATCGTCCCTCTCCCGTTTCATGGTGTTTTCCGCCGCCTCTAATCCCAATTTCAATATTCCAATGCCACCTTCATTCGTATCGCCAACGAAGCCCGCCGTGGc           <  7:235218/144‑1 (MQ=255)                                                                                                                                        cTTCACCGAAGGCGCCCAGGTCAAGAACGGCGACCTTCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTgg            <  7:317541/147‑1 (MQ=255)                                                                                                                                           cACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGa         <  1:64447/147‑1 (MQ=255)                                                                                                                                           cACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAgc       >  2:223289/1‑149 (MQ=255)                                                                                                                                            aCCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAgc       <  1:69581/148‑1 (MQ=255)                                                                                                                                              cGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGcc    >  7:227404/1‑149 (MQ=255)                                                                                                                                                aaGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTg  >  5:335482/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                               GCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCG‑GCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTG  >  NC_002947/3877011‑3877301

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