Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,219,674	C→T	6.5%	D56N (GAC→AAC)	PP_2821 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,219,674	0	C	T	6.5%	179.8 / 2.0	61	D56N (GAC→AAC)	PP_2821	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (2/2);  total (32/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 = 1.65e-01

ATCACCGCGGCGTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTGGCCAGGGTTGCGCCGAACATTGC  >  NC_002947/3219526‑3219821                                                                                                                                                      |                                                                                                                                                    aTCACCGCGGCGTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTc                                                                                                                                                     <  3:246472/149‑1 (MQ=255) aTCACCCCGCCCTCAGCGTCCACGCTCTGCCGCCGAGGCCCCTTGCCCGTGACCCCCACGCCTGACCAGTCCG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTc                                                                                                                                                     <  4:281390/149‑1 (MQ=255)    aCCGCGGCTTCATCGGCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCATCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAg                                                                                                                                                  <  8:171501/149‑1 (MQ=255)    aCCGCGGCGTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAg                                                                                                                                                  >  3:318341/1‑149 (MQ=255)     ccGCGGCGTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAgg                                                                                                                                                 <  8:153865/149‑1 (MQ=255)          gcgTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTACCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCtt                                                                                                                                            <  4:123885/149‑1 (MQ=255)            gTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGa                                                                                                                                          <  8:203037/149‑1 (MQ=255)                 cGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTg                                                                                                                                     >  3:313863/1‑149 (MQ=255)                  gTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTgg                                                                                                                                    >  1:173653/1‑149 (MQ=255)                   tCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGgt                                                                                                                                   >  5:185831/1‑149 (MQ=255)                     cACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGgtgt                                                                                                                                 >  2:203966/1‑149 (MQ=255)                        gCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACGTACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTcc                                                                                                                               >  1:92318/1‑149 (MQ=255)                          tctGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTg                                                                                                                            >  4:203789/1‑149 (MQ=255)                            tGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTTGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGa                                                                                                                                          <  5:160377/133‑1 (MQ=255)                            tGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGa                                                                                                                                          >  6:160377/1‑133 (MQ=255)                              cAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGcc                                                                                                                          <  2:92318/147‑1 (MQ=255)                                gCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAgg                                                                                                                      <  6:2923/149‑1 (MQ=255)                                         ccTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTg                                                                                                             <  4:318341/149‑1 (MQ=255)                                               ccGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATAGTTACTTGGTGGGGTTCAAGACTTTCTGTTCGAGGCCAGTGAGGATGGTGTCCTTGACAAGGAATATTTTGGACAcc                                                                                                       >  8:73726/1‑149 (MQ=255)                                                 gTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTc                                                                                                     >  1:35111/1‑149 (MQ=255)                                                  tGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCg                                                                                                    >  2:302041/1‑149 (MQ=255)                                                      cacCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCt                                                                                                <  6:199998/149‑1 (MQ=255)                                                       acCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTg                                                                                               <  8:250258/149‑1 (MQ=255)                                                         cAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTggg                                                                                             >  8:196979/1‑149 (MQ=255)                                                         cAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTgg                                                                                              >  2:113034/1‑148 (MQ=255)                                                             ccTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAg                                                                                                                                                  >  8:238708/1‑92 (MQ=255)                                                             ccTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAg                                                                                                                                                  <  7:238708/92‑1 (MQ=255)                                                                   cAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTgcggc                                                                                   >  8:6917/1‑149 (MQ=255)                                                                     gTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGgcgc                                                                                 >  1:294970/1‑149 (MQ=255)                                                                        cG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGgcgcca                                                                               <  1:203966/148‑3 (MQ=255)                                                                            catGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGgcgc                                                                           >  2:123891/3‑149 (MQ=255)                                                                               gAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCg                                                                         <  3:275801/148‑1 (MQ=255)                                                                                       aTTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCAt                                                                 >  7:66907/1‑148 (MQ=255)                                                                                       aTTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCAt                                                                 >  5:24820/1‑148 (MQ=255)                                                                                           gTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTTGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTgg                                                                                              >  5:58542/1‑115 (MQ=255)                                                                                           gTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTTGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTgg                                                                                              >  5:52750/1‑115 (MQ=255)                                                                                           gTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTTGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTgg                                                                                              <  6:52750/115‑1 (MQ=255)                                                                                           gTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTTGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTgg                                                                                              <  6:58542/115‑1 (MQ=255)                                                                                           gTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCa                                                            >  7:72972/1‑149 (MQ=255)                                                                                           gTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCa                                                            >  4:239239/1‑149 (MQ=255)                                                                                            ttCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCAc                                                           >  1:148886/1‑149 (MQ=255)                                                                                            ttCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCAc                                                           >  7:199045/1‑149 (MQ=255)                                                                                                tgttgGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGc                                                        <  8:199045/148‑1 (MQ=255)                                                                                                 gttgGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGCTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGc                                                        <  1:302041/147‑1 (MQ=255)                                                                                                     gTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTc                                                                   >  7:22710/1‑132 (MQ=255)                                                                                                     gTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTc                                                                   <  8:22710/132‑1 (MQ=255)                                                                                                     gTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCt                                                   <  3:203789/148‑1 (MQ=255)                                                                                                              cTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGc                                         <  1:203241/149‑1 (MQ=255)                                                                                                              cTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGc                                         <  7:6917/149‑1 (MQ=255)                                                                                                                gTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCa                                       <  6:185831/149‑1 (MQ=255)                                                                                                                 tAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAg                                      >  2:178101/1‑149 (MQ=255)                                                                                                                   gAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAggg                                    <  8:72972/149‑1 (MQ=255)                                                                                                                   gAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAgg                                     <  2:173653/148‑1 (MQ=255)                                                                                                                         ttCCTTGGTGGGGGTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTgctgg                              <  1:113034/149‑1 (MQ=255)                                                                                                                            cTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTgctggctg                           >  8:177580/1‑149 (MQ=255)                                                                                                                               ggtggGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCtgt                          <  6:56427/147‑1 (MQ=255)                                                                                                                                gtggGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTgg                        <  3:239239/148‑1 (MQ=255)                                                                                                                                   gggTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTc                                                                   >  6:103374/1‑102 (MQ=255)                                                                                                                                   gggTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTc                                                                   <  5:103374/102‑1 (MQ=255)                                                                                                                                           gCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTGGCCAg                    >  8:213951/1‑141 (MQ=255)                                                                                                                                           gCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTGGCCAg                    <  7:213951/141‑1 (MQ=255)                                                                                                                                                    gTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTGGCCAGGGTTGCGCCGAACATTg   >  1:79453/1‑149 (MQ=255)                                                                                                                                                     tCGAGGCCCTTGAGGCTGGTGTCCATGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTGGCCAGGGTTGCGCCGAACAtgcg  >  7:230703/1‑146 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ATCACCGCGGCGTCATCGTCCACGCTCTGCAGCGGTTGCCCCTTGCCCGTGAACACCAGGCCTGACCAGTACG‑ACTTGAGCTGGGATTCGTTCTTGTTGGTCAGCGCGCTGTAGAATCGTTCCTTGGTGGGGTTCCAGCCTTTCTGGTCGAGGCCCTTGAGGCTGGTGTCCTTGCCCAGGAATATGTTGGCCACCTCGGCCTGGGACGGTGCGGCGCTGGCGCCGGCATTCATGATCACGGCCACCTCGGCACTGGCCAGGGTGCTGGCTGTGGCCAGGGTTGCGCCGAACATTGC  >  NC_002947/3219526‑3219821

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