Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,378,046	T→C	6.1%	F32S (TTC→TCC)	PP_2978 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,378,046	0	T	C	6.1%	187.8 / 2.4	65	F32S (TTC→TCC)	PP_2978	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (32/33)
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

CAGGGTATTTGGCTGGATTGAGGAGGGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAG  >  NC_002947/3377900‑3378191                                                                                                                                                      |                                                                                                                                                    caGGGTATTTGGCTGGATTGAGGAGGGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTc                                                                                                                                                    >  1:222276/1‑148 (MQ=255)    ggTATTTGGCTGGATTGAGGAGGGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAg                                                                                                                                                >  8:82108/1‑149 (MQ=255)        tttGGCTGGATTGAGGAGGGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGa                                                                                                                                            <  7:82108/149‑1 (MQ=255)                gtttGAGGAGGGGGGGGTGCTTCGTGAATAAGTCTGCGGGAAATTCTAGGTGGATACCGGCGCAGTGACAGCGTCTATGGTTGCGAAAGATTTGGTGGGTTGCGTAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTCATGTGaaa                                                                                                                                     <  4:133030/146‑1 (MQ=255)                  ttGAGGAGGGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGa                                                                                                                                   >  2:146796/1‑148 (MQ=255)                     aggaggGGCTGCTGCTTACTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGt                                                                                                                               >  8:2041/1‑149 (MQ=255)                      ggaggGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGt                                                                                                                               >  5:209369/1‑148 (MQ=255)                       gaggGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGa                                                                                                                             <  3:222683/149‑1 (MQ=255)                        aggGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGa                                                                                                                             <  1:146796/148‑1 (MQ=255)                                ctgctTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGc                                                                                                                      <  5:75195/147‑1 (MQ=255)                                    ttCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGagca                                                                                                                >  5:393632/1‑149 (MQ=255)                                    ttCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGagca                                                                                                                >  6:408497/1‑149 (MQ=255)                                     tCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGagca                                                                                                                <  6:393632/148‑1 (MQ=255)                                        tGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGaa                                                                                                                                           >  1:121226/1‑118 (MQ=255)                                        tGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGaa                                                                                                                                           <  2:121226/118‑1 (MQ=255)                                         gAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTCCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGagcagc                                                                                                              >  3:179244/1‑146 (MQ=255)                                         gAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTCCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGcagccag                                                                                                           >  3:185464/1‑149 (MQ=255)                                          aaaaaGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGaa                                                                                                                                           >  4:245352/1‑116 (MQ=255)                                          aaaaaGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGaa                                                                                                                                           <  3:245352/116‑1 (MQ=255)                                                   cGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTc                                                                                                 >  2:333152/1‑149 (MQ=255)                                                    ggTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATGAAGTTGCGGTAGTGCTAAGGTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCtt                                                                                                   >  1:337868/1‑149 (MQ=255)                                                    ggTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTc                                                                                                 >  2:234332/1‑148 (MQ=255)                                                          aaTCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑Cg                                                                                          >  1:72970/1‑147 (MQ=255)                                                           aTCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑Cg                                                                                          >  4:108434/1‑146 (MQ=255)                                                           aTCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGaa                                                                                        <  3:108434/148‑1 (MQ=255)                                                                gCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGt                                                                                                                               >  5:209674/1‑106 (MQ=255)                                                                gCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGt                                                                                                                               <  6:209674/106‑1 (MQ=255)                                                                 cTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGa                                                                                 <  5:313163/149‑1 (MQ=255)                                                                    gATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGc                                                                               <  6:209369/148‑1 (MQ=255)                                                                                cTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGt                                                                  <  5:408497/149‑1 (MQ=255)                                                                                 tGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGc                                                                               <  8:206071/135‑1 (MQ=255)                                                                                 tGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGc                                                                               >  7:206071/1‑135 (MQ=255)                                                                                 tGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTc                                                                 >  7:37600/1‑149 (MQ=255)                                                                                        gTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCgtggtg                                                          <  2:72970/149‑1 (MQ=255)                                                                                                  gCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAg                                                                                                                                                <  7:7675/55‑1 (MQ=255)                                                                                                  gCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAg                                                                                                                                                >  8:7675/1‑55 (MQ=255)                                                                                                   ccAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCCGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCt                                               >  1:393546/1‑149 (MQ=255)                                                                                                    cAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGc                                                                               >  4:364543/1‑116 (MQ=255)                                                                                                    cAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGc                                                                               <  3:364543/116‑1 (MQ=255)                                                                                                    cAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCt                                               >  5:33825/1‑148 (MQ=255)                                                                                                    cAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCt                                               >  3:392347/1‑148 (MQ=255)                                                                                                    cAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCAGGCGGCGCt                                               >  3:195568/1‑148 (MQ=255)                                                                                                     aaaGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCt                                               >  4:36313/1‑147 (MQ=255)                                                                                                     aaaGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCt                                               >  4:102045/1‑147 (MQ=255)                                                                                                      aaGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGt                                                                        <  3:148986/121‑1 (MQ=255)                                                                                                      aaGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGt                                                                        >  4:148986/1‑121 (MQ=255)                                                                                                         aTTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCAc                                         <  6:1947/149‑1 (MQ=255)                                                                                                          tttGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAATTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCAc                                         <  2:222276/148‑1 (MQ=255)                                                                                                            tGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGt                                                                        <  6:103435/115‑1 (MQ=255)                                                                                                            tGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGt                                                                        >  5:103435/1‑115 (MQ=255)                                                                                                                  gTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTCCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                <  4:179244/149‑1 (MQ=255)                                                                                                                  gTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTCCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                <  4:185464/149‑1 (MQ=255)                                                                                                                         aaGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                         >  5:65631/1‑149 (MQ=255)                                                                                                                         aaGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                         <  3:36313/149‑1 (MQ=255)                                                                                                                         aaGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                         <  6:65631/149‑1 (MQ=255)                                                                                                                         aaGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCAGTGGTAGtga                         >  4:156234/1‑149 (MQ=255)                                                                                                                           ggAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtgatg                       <  2:337868/149‑1 (MQ=255)                                                                                                                                   g‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTTAAAGACAGTGACACTGGAAGAGCAGCCAGCATCCCTTTCTAGCTCGCCTGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCAGCGCTGAGCACCGCCCCGTGGTAGTGATGAAg                    <  1:333152/148‑1 (MQ=255)                                                                                                                                     cGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgt                <  4:392347/149‑1 (MQ=255)                                                                                                                                     cGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgt                <  2:355241/149‑1 (MQ=255)                                                                                                                                        gaAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCa             <  3:102045/149‑1 (MQ=255)                                                                                                                                          aaTGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCCGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATTTGCAtt           <  4:195568/149‑1 (MQ=255)                                                                                                                                              ggTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCa       <  1:234332/149‑1 (MQ=255)                                                                                                                                                tGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAAc     <  7:2041/149‑1 (MQ=255)                                                                                                                                                   gTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAg  >  5:368074/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CAGGGTATTTGGCTGGATTGAGGAGGGGCTGCTGCTTCCTGAAAAAGTCTCGGTGAAAATCAAGCTGGATACCGGCGCACTGACATCGTCTATGGATGCCAAAGATTTGGAGCGTTTCGAAAGGAATG‑‑G‑CGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAG‑‑CGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAG  >  NC_002947/3377900‑3378191

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