Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,391,521	G→C	6.8%	W244S (TGG→TCG)	PP_2994 →	FMN‑binding oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,391,521	0	G	C	6.8%	182.8 / 3.0	59	W244S (TGG→TCG)	PP_2994	FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAC  >  NC_002947/3391376‑3391659                                                                                                                                                     |                                                                                                                                           gAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCCGAGGCTTAACAACCCAAGCCGCGAAAAAGGTCGTTGAGTTGCt                                                                                                                                         >  4:241341/1‑149 (MQ=255)  aaaaCCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCt                                                                                                                                         >  3:5142/1‑148 (MQ=255)     aCCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAAc                                                                                                                                     <  7:373466/149‑1 (MQ=255)      ccGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAAc                                                                                                                                     >  6:254067/1‑148 (MQ=255)        gggCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGAc                                                                                                                                  <  4:294914/149‑1 (MQ=255)         ggCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGAcc                                                                                                                                 <  2:279763/149‑1 (MQ=255)            gcgcCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTgg                                                                                                                              >  4:32105/1‑149 (MQ=255)            gcgcCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTgg                                                                                                                              <  4:5142/149‑1 (MQ=255)                      cTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAt                                                                                                                    <  2:16728/149‑1 (MQ=255)                          aaaTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAAc                                                                                                                                     <  5:356943/128‑1 (MQ=255)                          aaaTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAAc                                                                                                                                     >  6:356943/1‑128 (MQ=255)                            aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                              <  2:30839/149‑1 (MQ=255)                            aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                              <  3:262475/149‑1 (MQ=255)                            aTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGGCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGa                                                                                                              <  2:308482/149‑1 (MQ=255)                             tCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAg                                                                                                             >  1:333381/1‑149 (MQ=255)                              cGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAg                                                                                                             <  3:32105/148‑1 (MQ=255)                               gTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGtt                                                                                                           <  3:37415/149‑1 (MQ=255)                               gTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGt                                                                                                            >  4:37415/1‑148 (MQ=255)                                 cAAGGCGGTCCGGGCCGTGGTTTCTGAGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCGGGGCGTa                                                                                                                         >  3:74208/1‑133 (MQ=255)                                    ggCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTgg                                                                                                                              <  5:190274/125‑1 (MQ=255)                                    ggCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTgg                                                                                                                              >  6:190274/1‑125 (MQ=255)                                      gggTCCGGGGCGTGGTTTCTGTGGTTTTTGCAGTGGCGGTCAACCTCACTACTTCCGCCGATTTCCCGCGCGGTGGGTTCACCTCCGAAGACGCGAAAAAGGGCGGTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCg                                                                                                       <  1:359310/148‑1 (MQ=255)                                       ggTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAg                                                                                                   >  4:205959/1‑149 (MQ=255)                                         tCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGc                                                                                                 <  7:289986/149‑1 (MQ=255)                                         tCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGc                                                                                                 <  3:217638/149‑1 (MQ=255)                                         tCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGAc                                                                                                                                  >  2:245635/1‑116 (MQ=255)                                         tCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGCCACTCAACGAc                                                                                                                                  <  1:245635/116‑1 (MQ=255)                                              gCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGc                                                                                              <  5:211466/147‑1 (MQ=255)                                              gCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGc                                                                                              >  6:211466/1‑147 (MQ=255)                                              gCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTa                                                                                            >  8:179156/1‑149 (MQ=255)                                                  tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                        <  7:24923/149‑1 (MQ=255)                                                  tGGGTTCTGTGTATTTTGCAGTGGCGTTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCAACGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                        <  8:272698/149‑1 (MQ=255)                                                      ttCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTcc                                                                                    >  3:170193/1‑149 (MQ=255)                                                               ttttGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGc                                                                           <  5:254067/149‑1 (MQ=255)                                                                   gCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAgg                                                                        <  8:381273/148‑1 (MQ=255)                                                                       tGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTg                                                                                                                               >  6:310954/1‑89 (MQ=255)                                                                       tGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTg                                                                                                                               <  5:310954/89‑1 (MQ=255)                                                                          cGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGAcgcg                                                                <  4:170193/149‑1 (MQ=255)                                                                            gTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGAcgcgcg                                                              >  4:343458/1‑149 (MQ=255)                                                                             tCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGAcgcgcgc                                                             >  5:162680/1‑149 (MQ=255)                                                                                 actca‑‑ac‑tcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATg                                                         <  7:179156/149‑1 (MQ=255)                                                                                  ctca‑‑ac‑tcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATgg                                                        >  8:141844/1‑149 (MQ=255)                                                                                    ca‑‑ac‑tcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTc                                                      >  5:71121/1‑149 (MQ=255)                                                                                    ca‑‑ac‑tcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTc                                                      >  7:258069/1‑149 (MQ=255)                                                                                     a‑‑ac‑tcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTc                                                      <  8:258069/148‑1 (MQ=255)                                                                                         c‑tcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCa                                                   >  5:190163/1‑149 (MQ=255)                                                                                                cgATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTgg                                             <  3:343458/149‑1 (MQ=255)                                                                                                          cgcgGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCGGGGCGTAGACATGGTGTAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGaa                                     <  3:241341/147‑1 (MQ=255)                                                                                                                        caccGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGt                     >  2:424215/1‑149 (MQ=255)                                                                                                                             aaGACGCGAAAAAGGTCGTTGAGTGGCTCAATGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGcc                >  5:319855/1‑149 (MQ=255)                                                                                                                                aCGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGt             >  1:23674/1‑149 (MQ=255)                                                                                                                                  gcgAATAAGGACGTTCAGTGGGTGAAGCACCTGGGCGTAGACATGGTGCAGATATGGGGACGCAGGTACGAAGTCCCCGGCATGCACGGTCACGCGCGCGATGGTCGGTCCGAGGCTCGCGAAGCTTAATTCGTCGCCAACGCCCGTGa           >  1:214761/1‑149 (MQ=255)                                                                                                                                  gcgAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTc                                                      <  7:332310/106‑1 (MQ=255)                                                                                                                                  gcgAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTc                                                      >  8:332310/1‑106 (MQ=255)                                                                                                                                  gcgAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGa           >  5:171189/1‑149 (MQ=255)                                                                                                                                    gAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtc         >  7:114098/1‑149 (MQ=255)                                                                                                                                    gAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtc         >  6:280754/1‑149 (MQ=255)                                                                                                                                           gTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAc  <  6:171189/149‑1 (MQ=255)                                                                                                                                           gTCGTTGAGCCACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAg                      >  2:374114/1‑129 (MQ=255)                                                                                                                                                     |                                                                                                                                           GAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCA‑‑AC‑TCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAGGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAC  >  NC_002947/3391376‑3391659

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