Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	262,752	T→A	6.4%	L90Q (CTG→CAG)	PP_0211 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	262,752	0	T	A	6.4%	189.0 / 3.2	62	L90Q (CTG→CAG)	PP_0211	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base A (2/2);  total (28/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 = 8.97e-01

CCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGC  >  NC_002947/262606‑262894                                                                                                                                                     |                                                                                                                                                ccGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTg                                                                                                                                                <  4:100944/148‑1 (MQ=255)    tGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGt                                                                                                                                              >  3:144245/1‑147 (MQ=255)       ggcgCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCa                                                                                                                                        <  4:43737/149‑1 (MQ=255)        gcgcgcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCaa                                                                                                                                       <  4:144245/149‑1 (MQ=255)        gcgcgcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCaa                                                                                                                                       <  6:98559/149‑1 (MQ=255)         cgcgcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGAGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAg                                                                                                                                      <  5:68473/149‑1 (MQ=255)             ccTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTg                                                                                                                                    >  5:208281/1‑147 (MQ=255)              cTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAg                                                                                                                                 <  6:159453/149‑1 (MQ=255)              cTGGCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGC‑T‑GACCAGGTCACGCAAGtaagc                                                                                                                                 <  2:366855/149‑5 (MQ=255)                gTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCt                                                                                                                               >  6:130514/1‑149 (MQ=255)                   cGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                            <  6:208281/149‑1 (MQ=255)                    ggAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATc                                                                                                                           >  5:369400/1‑149 (MQ=255)                     gAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                            >  2:280971/1‑147 (MQ=255)                     gAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                            <  1:280971/147‑1 (MQ=255)                      aGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATccc                                                                                                                         <  7:243829/149‑1 (MQ=255)                                cTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGcc                                                                                                               <  5:130514/149‑1 (MQ=255)                                 tGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCg                                                                                                              >  6:291549/1‑149 (MQ=255)                                   gTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAAGTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCgg                                                                                                             <  3:206983/148‑1 (MQ=255)                                   gTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCaa                                                                                                                                       <  6:119217/122‑1 (MQ=255)                                   gTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCgg                                                                                                             >  5:131201/1‑148 (MQ=255)                                      cgcgACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAg                                                                                                          <  2:72903/148‑1 (MQ=255)                                      cgcgACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAg                                                                                                          <  6:10576/148‑1 (MQ=255)                                          acgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGTTCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGt                                                                                                                                     <  8:154158/117‑1 (MQ=255)                                          acgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGt                                                                                                                                     >  7:154158/1‑117 (MQ=255)                                                tgCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGaa                                                                                                <  2:176355/148‑1 (MQ=255)                                                 gcgcGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGTATGGGCATGAGCGCGGGATCTACCATTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGATGACCCCTTCAGTTGCCGGAAGGGGTGCAg                                                                                                  >  1:282801/1‑146 (MQ=255)                                                   gcgTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAg                                                                                            <  4:316748/149‑1 (MQ=255)                                                   gcgTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAg                                                                                            >  3:303355/1‑149 (MQ=255)                                                       ggtgCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGgac                                                                                         >  5:47507/1‑148 (MQ=255)                                                            aaaTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGc                                                                                   >  6:317961/1‑149 (MQ=255)                                                                       cAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACg                                                                         >  1:148356/1‑148 (MQ=255)                                                                       cAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACg                                                                         >  4:141982/1‑148 (MQ=255)                                                                         aGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAAc                                                                      <  6:47507/149‑1 (MQ=255)                                                                             tATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGc                                                                  >  2:44128/1‑149 (MQ=255)                                                                              aTGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGcc                                                                 >  5:207927/1‑149 (MQ=255)                                                                               tGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGcc                                                                 <  6:131201/148‑1 (MQ=255)                                                                                 gCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGa                                                               <  1:192760/148‑1 (MQ=255)                                                                                 gCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGa                                                               <  5:291549/148‑1 (MQ=255)                                                                                 gCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAg                                                              >  2:211170/1‑149 (MQ=255)                                                                                  cGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGACCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAgg                                                             <  2:282801/149‑1 (MQ=255)                                                                                   gTACAACTGGTGTGCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGCCCGCCTAGCCAACGAACTGGCCGAggg                                                            >  4:107636/1‑149 (MQ=255)                                                                                    tACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGt                                                           <  4:303355/149‑1 (MQ=255)                                                                                      cAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTaa                                                         >  4:162045/1‑149 (MQ=255)                                                                                                 cAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCttt                                               >  1:82292/1‑149 (MQ=255)                                                                                                   gCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGTTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGc                                            <  2:82292/149‑1 (MQ=255)                                                                                                   gCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGTACTGGCCGAGGGTAACCTGCTCTTTGGc                                            <  1:211170/149‑1 (MQ=255)                                                                                                    cG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCg                                           <  6:207927/149‑1 (MQ=255)                                                                                                            cATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACAcccc                                    >  6:304312/1‑149 (MQ=255)                                                                                                                 gcgcGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAAc                                                        >  7:219110/1‑124 (MQ=255)                                                                                                                 gcgcGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAAc                                                        <  8:219110/124‑1 (MQ=255)                                                                                                                    cgGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                            <  1:44128/149‑1 (MQ=255)                                                                                                                    cgGGAGCTACCCTTGCGCTTACTTGCGT‑GACCTGGGCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                            <  6:369400/149‑1 (MQ=255)                                                                                                                            aCCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAgg                                                             <  8:137857/108‑1 (MQ=255)                                                                                                                            aCCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAgg                                                             >  7:137857/1‑108 (MQ=255)                                                                                                                             cccTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGgc                   >  5:272181/1‑149 (MQ=255)                                                                                                                                   gcttcaCTTGC‑T‑GACCAGGTCACGCAAGT‑AAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg            >  8:22392/6‑149 (MQ=255)                                                                                                                                    gCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg            >  8:188966/1‑149 (MQ=255)                                                                                                                                    cttcaCTTGC‑T‑GACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg            >  5:322592/5‑149 (MQ=255)                                                                                                                                         cTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCaa       >  8:39509/1‑149 (MQ=255)                                                                                                                                             cGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGt    <  7:188966/148‑1 (MQ=255)                                                                                                                                              gT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGc  <  5:304312/149‑1 (MQ=255)                                                                                                                                              ct‑GACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGt    <  6:322592/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGC  >  NC_002947/262606‑262894

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