Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,562,071	C→T	6.1%	S154N (AGC→AAC)	PP_3145 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,562,071	0	C	T	6.1%	195.2 / 2.4	65	S154N (AGC→AAC)	PP_3145	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base T (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTG  >  NC_002947/3561937‑3562217                                                                                                                                        |                                                                                                                                                   tCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCgcgg                                                                                                                                      >  7:5055/1‑149 (MQ=255)   acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                           <  6:2842/142‑4 (MQ=255)   acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                           >  5:2842/1‑139 (MQ=255)          ccGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGcc                                                                                                                             >  3:93573/1‑149 (MQ=255)             cgATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCg                                                                                                                          >  2:159385/1‑149 (MQ=255)                tAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGtt                                                                                                                        >  7:131758/1‑148 (MQ=255)                  aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                     <  6:162985/149‑1 (MQ=255)                  aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                     <  7:20793/149‑1 (MQ=255)                  aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGa                                                                                                                      <  8:132144/148‑1 (MQ=255)                   gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                                <  8:40139/137‑1 (MQ=255)                   gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                                >  7:40139/1‑137 (MQ=255)                       gATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                                <  1:159385‑M1/149‑1 (MQ=255)                       gATAGGGATCACCTCGTAAGGCCGGAGGCCCTCGTTACTGGGGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                                <  7:143815‑M1/149‑1 (MQ=255)                            ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                          >  7:170960/1‑114 (MQ=255)                            ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                          <  8:170960/118‑5 (MQ=255)                             gATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttat                                                                                                          >  5:1028‑M1/1‑145 (MQ=255)                              aTCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGAt                                                                                                         <  4:93573/149‑1 (MQ=255)                                 aCCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑CCCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTc                                                                                                      >  2:318708/1‑149 (MQ=255)                                  ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                            <  4:271667/142‑1 (MQ=255)                                  ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                            >  3:271667/1‑142 (MQ=255)                                   cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                                    <  7:11254/149‑1 (MQ=255)                                   cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                                    >  1:216897/1‑149 (MQ=255)                                       tAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCGACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTg                                                                                                   >  1:241985/1‑147 (MQ=255)                                        aaGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                               >  6:285762/1‑149 (MQ=255)                                         aGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCGCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                               <  5:285762/148‑1 (MQ=255)                                           gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                     <  7:224920/124‑1 (MQ=255)                                           gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                     >  8:224920/1‑124 (MQ=255)                                                  gCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCtgctgc                                                                                      >  7:28615/1‑148 (MQ=255)                                                     cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                               >  4:252782‑M1/1‑120 (MQ=255)                                                     cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                               <  3:252782‑M1/120‑1 (MQ=255)                                                        gTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATa                                                                                <  8:5055/148‑1 (MQ=255)                                                           aCTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAAccg                                                                            <  5:247132‑M2/149‑4 (MQ=255)                                                             tggtggACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGAGAAccgcg                                                                          <  8:80284‑M2/149‑6 (MQ=255)                                                               gtggACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGt                                                                        <  3:287060/149‑1 (MQ=255)                                                                tggACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGtt                                                                       <  2:241985/149‑1 (MQ=255)                                                                  gACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTgg                                                                     <  3:153293/149‑1 (MQ=255)                                                                   acacTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcacccgctacttca                                                                    <  6:1028‑M1/149‑43 (MQ=255)                                                                       tGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                               >  1:92087‑M1/1‑103 (MQ=255)                                                                       tGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                               <  2:92087‑M1/134‑32 (MQ=255)                                                                         ccGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCgg                                                               <  1:309204/148‑1 (MQ=255)                                                                            tGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                           >  2:228375/1‑149 (MQ=255)                                                                                 gCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTa                                                      <  8:28615/149‑1 (MQ=255)                                                                                  ccGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                                        >  7:43944/1‑146 (MQ=255)                                                                                  ccGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                                        >  7:183866/1‑146 (MQ=255)                                                                                    gggCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGaa                                                   >  6:253517/1‑149 (MQ=255)                                                                                     ggCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAAc                                                  <  2:216897/149‑1 (MQ=255)                                                                                        gc‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                           >  8:101645/1‑137 (MQ=255)                                                                                        gc‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                           <  7:101645/137‑1 (MQ=255)                                                                                                        ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGaa                                 >  5:6853/1‑148 (MQ=255)                                                                                                        ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAac                                <  8:131758/149‑1 (MQ=255)                                                                                                          gggCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAacac                              <  8:183866/149‑1 (MQ=255)                                                                                                              cAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTg                          <  1:228375/149‑1 (MQ=255)                                                                                                                   cccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTa                     >  8:193783/1‑149 (MQ=255)                                                                                                                    ccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTAc                    <  6:247132/149‑1 (MQ=255)                                                                                                                          ccAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTg              <  8:43944/149‑1 (MQ=255)                                                                                                                             aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata           <  1:318708/149‑1 (MQ=255)                                                                                                                             aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata           <  6:6853/149‑1 (MQ=255)                                                                                                                                 cATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa       >  7:196228/1‑149 (MQ=255)                                                                                                                                    cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTTATATTGa       >  2:217183/1‑146 (MQ=255)                                                                                                                                    cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa       <  1:212564/146‑1 (MQ=255)                                                                                                                                    cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa       >  2:212564/1‑146 (MQ=255)                                                                                                                                    cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa       <  1:217183/146‑1 (MQ=255)                                                                                                                                      cGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAacac                              >  7:95514/1‑121 (MQ=255)                                                                                                                                      cGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAacac                              <  8:95514/121‑1 (MQ=255)                                                                                                                                      cGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTg  <  8:196228/149‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                   TCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGC‑ACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTG  >  NC_002947/3561937‑3562217

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