Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	101,802	A→C	6.3%	K208Q (AAG→CAG) ‡	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	101,802	0	A	C	6.3%	196.3 / 4.0	63	K208Q (AAG→CAG) ‡	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base C (2/2);  total (30/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 = 9.48e-01

AGAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGA  >  NC_002947/101656‑101950                                                                                                                                                   |                                                                                                                                                     agagCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAg                                                                                                                                                    <  6:194422/149‑1 (MQ=255)  gagCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                   <  3:54393/149‑1 (MQ=255)      aaGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCt                                                                                                                                               <  4:311787/149‑1 (MQ=255)      aaGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCt                                                                                                                                               <  8:61505/149‑1 (MQ=255)       aGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTg                                                                                                                                              >  5:373413/1‑149 (MQ=255)         cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                            <  1:170249/149‑2 (MQ=255)            aGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGACTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACgg                                                                                                                                         >  1:306307/1‑149 (MQ=255)              cGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCt                                                                                                                                       <  8:114956/149‑1 (MQ=255)                aGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTgg                                                                                                                                     <  6:373413/149‑1 (MQ=255)                  cTGGGCAGCGGTTTCTCCAATCAGCTGCTGGACGCCACCCTGGCTTGGACCAAGCATGTTACCTACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCATGGCGCATTTGGCCGCTGCAGCCCAGGCCATGGGCCTGGACGGCTGGt                                                                                                                                    <  8:202577/148‑1 (MQ=255)                  cTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGtt                                                                                                                                   <  6:334394/149‑1 (MQ=255)                     ggCAGCCGTTTCTCAAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACTAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGa                                                                                                                                 <  3:214138/148‑1 (MQ=255)                          ccGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGTCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAccc                                                                                                                           >  5:84026/1‑149 (MQ=255)                            gTTTCTCCAATCCGGTGCGGGCCGCCCCCCCGCCGTGCACCAAGCTTGTTACCGCCCACGCGGCGCTGGCCGGCCTTACCGATTCGCCCAAGGCGCAGATGGCCGCTGCCCCCCCGGCCAAGGGCCTGGACGGCTGGTTGATCACCCtt                                                                                                                         <  2:168536/149‑1 (MQ=255)                                  tataTCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACCAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGGCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAAtt                                                                                                                    <  3:56079/145‑1 (MQ=255)                                   cAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTcc                                                                                                                  >  4:106443/1‑149 (MQ=255)                                    aaTCAGCTGCTGCACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCGGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTccc                                                                                                                 <  7:301615/149‑1 (MQ=255)                                     aTCTGCTGCTGGACGCCACCCAGGCCTGGACCAAGAATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCATGGCGCAGATGGCCGCTGCAGCCCTGGCCATGGGCCTGGACGGCTGGTTGATCACCCTTGAATTcccc                                                                                                                <  8:309325/149‑1 (MQ=255)                                      tCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAAtt                                                                                                                    >  4:56079/1‑144 (MQ=255)                                         gctgctGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGc                                                                                                             >  4:321242/1‑148 (MQ=255)                                                gACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGc                                                                                                             <  6:87356/141‑1 (MQ=255)                                                gACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGc                                                                                                             >  5:87356/1‑141 (MQ=255)                                                  cGCCCCCCCGGACTGGACCATGAATGTTACCTACTAAGCGTGGCGGGCCGGCGTGACGGATTGGGCCATGGCGCAGATGGCCGCTGCAGCCCAGGCCATGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCgg                                                                                                   <  8:251146/149‑1 (MQ=255)                                                   gCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                  <  2:142091/149‑1 (MQ=255)                                                      aCCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTccc                                                                                                                 <  6:28816/131‑1 (MQ=255)                                                      aCCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTccc                                                                                                                 >  5:28816/1‑131 (MQ=255)                                                      aCCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgat                                                                                               <  4:83718/149‑1 (MQ=255)                                                       cccAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgatg                                                                                              <  1:349809/149‑1 (MQ=255)                                                          aGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtcgctggc                                                                                           >  7:139475/1‑142 (MQ=255)                                                          aGGCCTGGACAAAGCATGTTCCCGACGAAGCTGCGCTGGCCGGCCTGACCTATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGAcc                                                                                           <  5:31220/149‑1 (MQ=255)                                                               ggtaCCACGCATGTTACCGACGACGGTGCGCGGGCCGGCCTGACCGTTTCGGCCAAGGCGCTTTGGGCCGCTGCCCCCCAGGCCATGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGc                                                                                      <  2:374729/146‑1 (MQ=255)                                                                    cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCtactac                                                                                                       <  1:47175/127‑1 (MQ=255)                                                                    cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCtactac                                                                                                       >  2:47175/1‑127 (MQ=255)                                                                    cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                  >  7:344807/1‑132 (MQ=255)                                                                    cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                  <  8:344807/132‑1 (MQ=255)                                                                    cAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCg                                                                                 >  1:344871/1‑149 (MQ=255)                                                                                  acgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcgc                                                                   <  2:309005/149‑1 (MQ=255)                                                                                     aaGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGaa                                                                <  4:210735/149‑1 (MQ=255)                                                                                           cgcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCta                                                          >  8:31202/1‑149 (MQ=255)                                                                                                gccggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgc                                                     >  2:214956/1‑149 (MQ=255)                                                                                                       tGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTAc                                              >  5:57544/1‑149 (MQ=255)                                                                                                        gACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                             >  4:148013/1‑149 (MQ=255)                                                                                                        gACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                             >  1:239217/1‑149 (MQ=255)                                                                                                         aCCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCgg                                                                                                   <  5:25748/94‑1 (MQ=255)                                                                                                         aCCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCgg                                                                                                   >  6:25748/1‑94 (MQ=255)                                                                                                           cGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGACTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCa                                          >  8:9809/1‑149 (MQ=255)                                                                                                             aTTAGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGTCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGAGACCGCGCCCTGCGCGAAGAGCTATACGACTCCTACTGCAcc                                        >  7:146916/1‑149 (MQ=255)                                                                                                                         cgcAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGAc                            <  2:306307/149‑1 (MQ=255)                                                                                                                            aGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAg                         >  3:329155/1‑149 (MQ=255)                                                                                                                             gATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAgg                        <  3:321242/149‑1 (MQ=255)                                                                                                                              aTGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAggg                       >  8:100784/1‑149 (MQ=255)                                                                                                                              aTGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAggg                       >  6:95558/1‑149 (MQ=255)                                                                                                                               tGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGc                      >  4:359592/1‑149 (MQ=255)                                                                                                                                  ccGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCg                   <  5:181465/149‑1 (MQ=255)                                                                                                                                      agccgtccAGGCCCTTGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGcc                     <  5:186959/137‑1 (MQ=255)                                                                                                                                      agccgtccAGGCCCTTGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGcc                     >  6:186959/7‑143 (MQ=255)                                                                                                                                       gCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCACTGCGCGAAGAGCTATACGCCGGCTACAGCACCCGCGCCTCCGACCAGGGCCCGAATga              >  4:246963/1‑148 (MQ=255)                                                                                                                                        cAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGcc             >  6:336052/1‑149 (MQ=255)                                                                                                                                           tccAGGCCCTTGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGc          >  3:151783/2‑149 (MQ=255)                                                                                                                                             cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGCGCCCGAATGCCGGc          >  8:155493/1‑147 (MQ=255)                                                                                                                                                gCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGtt     <  3:148013/149‑1 (MQ=255)                                                                                                                                                  cAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc    >  6:257492/1‑148 (MQ=255)                                                                                                                                                  cAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCg   <  6:84026/149‑1 (MQ=255)                                                                                                                                                   aaGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGa  >  3:131186/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AGAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGA  >  NC_002947/101656‑101950

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑