Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	557,065	A→T	7.9%	P38P (CCA→CCT)	rpsN →	30S ribosomal protein S14

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	557,065	0	A	T	7.9%	162.8 / 2.1	59	P38P (CCA→CCT)	rpsN	30S ribosomal protein S14
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base T (2/3);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTCAAG  >  NC_002947/556920‑557210                                                                                                                                                   |                                                                                                                                                      gTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCagaa                                                                                                                                                    <  7:51698/149‑1 (MQ=255)   tCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCagaa                                                                                                                                                    >  1:107894/1‑147 (MQ=255)   tCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCagaa                                                                                                                                                    >  8:175107/1‑147 (MQ=255)        ttCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGtt                                                                                                                                             >  6:103854/1‑149 (MQ=255)            gCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGc                                                                                                                                         >  1:317059/1‑149 (MQ=255)                    ttGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTg                                                                                                                                 <  7:69855/149‑1 (MQ=255)                       gAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCg                                                                                                                              <  5:279851/149‑1 (MQ=255)                          tAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGc‑tc                                                                                                                          >  1:148401/1‑149 (MQ=255)                              aaaaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCa                                                                                                                      <  5:103854/149‑1 (MQ=255)                               aaaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCa                                                                                                                      <  2:148401/148‑1 (MQ=255)                                    gCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAa‑‑‑g                                                                                                               >  7:243785/1‑147 (MQ=255)                                      cAATAAGAGCATGAAAAACCGCTAGCTGAAGCGTCAGCTCACCACCAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAa‑‑‑gcag                                                                                                            <  1:110523/149‑1 (MQ=255)                                       aagaagAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGcc                                                                                                          <  2:50430/149‑1 (MQ=255)                                         gaagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCAc                                                                                                        >  2:263779/1‑149 (MQ=255)                                         gaagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCAc                                                                                                        <  2:317059/149‑1 (MQ=255)                                                  gAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGc                                                                                                                       >  1:56323/1‑128 (MQ=255)                                                  gAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGc                                                                                                                       <  2:56323/128‑1 (MQ=255)                                                  gAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGc                                                                                                 >  4:124932/1‑147 (MQ=255)                                                  gAAAAACCGCGAGCTGAAGCGACAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCAGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGACCGTTTCGCTGCCGTTGTCGC‑TCTTCAGAA‑‑‑GCAGCCACGCGACGcaa                                                                                               >  7:73509/1‑147 (MQ=255)                                                   aaaaaCCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTa                                                                                               <  2:107894/148‑1 (MQ=255)                                                    aaaaCCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAgc                                                                                             <  2:138026/149‑1 (MQ=255)                                                     aaaCCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAgcg                                                                                            >  1:210143/1‑149 (MQ=255)                                                         cgtgAGCTGACGCTTCAGCTCG‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTc                                                                                        >  1:282142/4‑149 (MQ=255)                                                                 gAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTgc                                                                                >  1:195314/1‑149 (MQ=255)                                                                                     gCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCTGAAGAGCGTTTCGCTGCCGTTGTCGc‑tct                                                                                                                         >  8:268212/1‑92 (MQ=255)                                                                                     gCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCTGAAGAGCGTTTCGCTGCCGTTGTCGc‑tct                                                                                                                         <  7:268212/92‑1 (MQ=255)                                                                                     gCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCTGAAGAGCGTTTCGCTGCCGTTGTCGc‑tct                                                                                                                         <  7:263543/92‑1 (MQ=255)                                                                                     gCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCTGAAGAGCGTTTCGCTGCCGTTGTCGc‑tat                                                                                                                         <  7:268215/92‑3 (MQ=255)                                                                                     gCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGc‑tct                                                                                                                         >  8:263543/1‑92 (MQ=255)                                                                                     gCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCAGCGGCAACCTGAACGCCTCGCCTGAAGAGCGTTTCGCTGCCGTTTTCGc‑tct                                                                                                                         >  8:268215/1‑92 (MQ=255)                                                                                      cAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGccc                                                            >  5:255989/1‑149 (MQ=255)                                                                                       aaaGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCg                                                           >  7:235868/1‑149 (MQ=255)                                                                                       aaaGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCg                                                           <  2:282142/149‑1 (MQ=255)                                                                                                 aagaagCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCa                                                 >  5:52051/1‑149 (MQ=255)                                                                                                  agaagCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCAc                                                >  7:150648/1‑149 (MQ=255)                                                                                                  agaagCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTACCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCAc                                                <  6:51657/149‑1 (MQ=255)                                                                                                      gCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGgtg                                            <  7:175107/149‑1 (MQ=255)                                                                                                       cGTGCGGAGCGGAAAGCGACCATCGTCACCCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCATAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGgtgt                                           <  8:243785/149‑1 (MQ=255)                                                                                                        gTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAACCCTCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGgt                                             <  1:263779/149‑1 (MQ=255)                                                                                                          gctgagctgaAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGtata                                        <  6:52051/149‑1 (MQ=255)                                                                                                                 tgaAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAg                                 >  5:327561/1‑149 (MQ=255)                                                                                                                  gaAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCATAAGt                                >  5:77268/1‑149 (MQ=255)                                                                                                                    aaGCGACCATCCTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGtt                               <  8:226090/148‑1 (MQ=255)                                                                                                                           cATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTgg                       >  1:236978/1‑149 (MQ=255)                                                                                                                              cGTCAACCTGAACGCCTCTCCAGTAGAGAGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGTCCGGTCGTCCTCACGGTGTTTACCGTAAGTTCGGCCTGGGc                     <  8:73509/148‑1 (MQ=255)                                                                                                                              cGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGcc                    >  6:37404/1‑149 (MQ=255)                                                                                                                                tCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCg                   >  5:297718/1‑148 (MQ=255)                                                                                                                                 cAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTa                 <  8:150648/149‑1 (MQ=255)                                                                                                                                   aCCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTaa                >  8:108471/1‑148 (MQ=255)                                                                                                                                         aCGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACAtgctg         <  5:37404/149‑1 (MQ=255)                                                                                                                                         aCGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACAtgct          <  8:235868/148‑1 (MQ=255)                                                                                                                                            ccTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGt      >  2:128021/1‑149 (MQ=255)                                                                                                                                             ctctCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGgt                                             <  2:186925/109‑1 (MQ=255)                                                                                                                                             ctctCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGgt                                             >  1:186925/1‑109 (MQ=255)                                                                                                                                             ctctCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGt      >  4:224754/1‑148 (MQ=255)                                                                                                                                             ctctCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGt      >  3:12829/1‑148 (MQ=255)                                                                                                                                               ctcCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTc     >  4:114255/1‑147 (MQ=255)                                                                                                                                               ctcCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTc     >  8:195916/1‑147 (MQ=255)                                                                                                                                                tcCAGAAGAGCGTTTCGCTGCCGTTGTCGCATCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTc     >  1:28983/1‑147 (MQ=255)                                                                                                                                                tcCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTCAAg  >  8:199500/1‑149 (MQ=255)                                                                                                                                                  cAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTCaa   <  3:124932/146‑1 (MQ=255)                                                                                                                                                  cAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCATAAGTTCGGCCTGGGCTGTAACATGCTGCGTCAAg  <  6:77268/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCA‑CGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGC‑TCTGCAGAA‑‑‑GCAGCCACGCGACGCTAGCGCCTCGCGCCTGCGCAACCGTTGCCGCCTGACCGGTCGTCCTCACGGTGTATACCGTAAGTTCGGCCTGGGCCGTAACATGCTGCGTCAAG  >  NC_002947/556920‑557210

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