Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,475,315	C→T	6.4%	Y72Y (TAC→TAT)	PP_4812 →	3‑methyladenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,475,315	0	C	T	6.4%	182.6 / 2.2	62	Y72Y (TAC→TAT)	PP_4812	3‑methyladenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ATCCGTCACCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGATGCAA  >  NC_002947/5475172‑5475463                                                                                                                                                 |                                                                                                                                                     aTCCTTCACCGCCACGGCGAACTCTGGCTGGCCGCGCGGATCATCCAGACAGAGGCCTACTACCTATGCGACAAAGGTAGCCACGCTTCGCTCGGCAACACCGAAAAACTCAAGGCGCTA‑GTCGTGGATGGGGGGCACATCAACaagg                                                                                                                                                  >  8:83417/1‑145 (MQ=255) aTCAGTCACCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACGACCTTTCCGACAAACGTAGCCACGCTTCGCTCGGCTCCACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCGACAGGtt                                                                                                                                                 >  8:83421/1‑148 (MQ=255)  tCCGTCACCGCCACGGCCACCTGTGGCTGGCCGCGCGGATCATCGAGTCCGTGGCCTTCTACCTTTCCGACAAAGGTAGCCACGGTTCGCTCGGCTTCACCGAAAAACGCAAGGCGCTACTTCCTGGATGGGGGGCACATCTACATGta                                                                                                                                                 <  1:198445/149‑1 (MQ=255)    cGTCACCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGTCCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGtact                                                                                                                                               <  6:247088/148‑1 (MQ=255)     gTCACCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACCTCTGCATGtactac                                                                                                                                             >  3:176499/1‑149 (MQ=255)     gTCACCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGtactac                                                                                                                                             <  1:173060/149‑1 (MQ=255)        aCCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGta                                                                                                                                                 <  1:256584/142‑1 (MQ=255)        aCCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGta                                                                                                                                                 >  2:256584/1‑142 (MQ=255)                         tggctggcCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGc                                                                                                                          <  3:320859/148‑1 (MQ=255)                                 cgcgcgGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCa                                                                                                                 <  3:277116/149‑1 (MQ=255)                                     cgGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTATATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAgcgc                                                                                                             >  4:265766/1‑149 (MQ=255)                                     cgGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTATATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAgcgc                                                                                                             >  4:270204/1‑149 (MQ=255)                                      ggATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGcc                                                                                                            <  3:48090/149‑1 (MQ=255)                                            tcGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCAcggg                                                                                                      >  8:272677/1‑149 (MQ=255)                                            acGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCAcggg                                                                                                      >  6:55886/2‑149 (MQ=255)                                             cGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCAcgggc                                                                                                     >  5:311857/1‑149 (MQ=255)                                               agaCCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCAcgggccg                                                                                                   <  6:311857/149‑1 (MQ=255)                                                 aCCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCAcgggccggg                                                                                                 <  3:38085/149‑1 (MQ=255)                                                  ccGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCAcgggccggg                                                                                                 >  8:5759/1‑148 (MQ=255)                                                        ccTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAAcgcg                                                                                          <  7:42679/149‑1 (MQ=255)                                                          tactacCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGt                                                                                        >  2:306006/1‑149 (MQ=255)                                                          tactacCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGt                                                                                        >  3:155087/1‑149 (MQ=255)                                                            ctacCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTg                                                                                       <  7:272677/148‑1 (MQ=255)                                                            ctacCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGc                                                                                      <  1:245240/149‑1 (MQ=255)                                                                 tttCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGAt                                                                                  >  1:8730/1‑148 (MQ=255)                                                                    ccGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCaa                                                                               <  5:329678/148‑1 (MQ=255)                                                                    ccGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCaa                                                                               <  4:155087/148‑1 (MQ=255)                                                                       aCAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCt                                                                           >  1:268415/1‑149 (MQ=255)                                                                         aaaGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGc                                                                         >  7:183293/1‑149 (MQ=255)                                                                         aaaGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGc                                                                         <  8:183293/149‑1 (MQ=255)                                                                         aaaGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGc                                                                         >  2:230652/1‑149 (MQ=255)                                                                                gCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCCCGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCt                                                                  >  8:312520/1‑149 (MQ=255)                                                                                              ggCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTATATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACAcc                                                      <  3:270204/147‑1 (MQ=255)                                                                                              ggCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTATATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACAcc                                                      <  3:265766/147‑1 (MQ=255)                                                                                                  acacCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCg                                                >  6:254709/1‑149 (MQ=255)                                                                                                   cacCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCtctc                                                 <  6:21020/147‑1 (MQ=255)                                                                                                   cacCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCgg                                               >  7:213388/1‑149 (MQ=255)                                                                                                    acCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCggg                                              >  7:68609/1‑149 (MQ=255)                                                                                                     ccGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTgg                                                               >  4:151078/1‑131 (MQ=255)                                                                                                     ccGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTgg                                                               <  3:151078/131‑1 (MQ=255)                                                                                                     ccGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCgggg                                             <  1:294992/149‑1 (MQ=255)                                                                                                      cGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGc                                            >  3:170268/1‑149 (MQ=255)                                                                                                      cGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGc                                            >  8:194422/1‑149 (MQ=255)                                                                                                       gAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGcc                                           <  8:213388/149‑1 (MQ=255)                                                                                                            aCGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACa                                      <  5:55886/149‑1 (MQ=255)                                                                                                              gCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGc                                    <  7:5759/149‑1 (MQ=255)                                                                                                               cAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCctct                                                  >  5:261291/1‑134 (MQ=255)                                                                                                               cAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCctct                                                  <  6:261291/134‑1 (MQ=255)                                                                                                                    cgcTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGGCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGgcg                              <  4:170268/149‑1 (MQ=255)                                                                                                                      cTA‑TTCCTGGCTGGGGTGCACATCTACATGTACTACGCCCGCGGCGGCGCTTCGCGGAACTTCAGCGCCCACGGTCCTGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGTCGCa                            <  4:176499/149‑1 (MQ=255)                                                                                                                          ttCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGAt                         <  1:230652/149‑1 (MQ=255)                                                                                                                             cTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACAcc                                                      >  1:140346/1‑117 (MQ=255)                                                                                                                             cTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACAcc                                                      <  2:140346/117‑1 (MQ=255)                                                                                                                                      gggCACATCTACCTGTACTACCCCCGCGGCGGCGATGCGCTGAACTTCAGCTCCCACGGGCCGGGCAACGCGGGGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGaacac             <  7:83421/149‑2 (MQ=255)                                                                                                                                      gggCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGATCTTCAGCGGCCACGGGAGGGGCAACGCGGTGCTGATCAAGTCTGCCTATGCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCTATTGaacaa             <  7:83417/149‑1 (MQ=255)                                                                                                                                      gggCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGaacaa             >  8:23614/1‑149 (MQ=255)                                                                                                                                            aTCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGc                                                                         >  6:30943/1‑83 (MQ=255)                                                                                                                                            aTCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGc                                                                         <  5:30943/83‑1 (MQ=255)                                                                                                                                             tCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGAt      >  8:76889/1‑149 (MQ=255)                                                                                                                                               tACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGATg     >  1:206086/1‑148 (MQ=255)                                                                                                                                                aCATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTgg                                <  6:75628/120‑1 (MQ=255)                                                                                                                                                aCATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTgg                                >  5:75628/1‑120 (MQ=255)                                                                                                                                                 cATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGAt      <  7:194422/145‑1 (MQ=255)                                                                                                                                                 cATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGATGc    >  1:165045/1‑147 (MQ=255)                                                                                                                                                 cATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGATGCaa  <  2:268415/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ATCCGTCACCGCCACGGCGACCTGTGGCTGGCCGCGCGGATCATCGAGACCGAGGCCTACTACCTTTCCGACAAAGGTAGCCACGCTTCGCTCGGCTACACCGAAAAACGCAAGGCGCTA‑TTCCTGGATGGGGGGCACATCTACATGTACTACGCCCGCGGCGGCGATTCGCTGAACTTCAGCGCCCACGGGCCGGGCAACGCGGTGCTGATCAAGTCTGCCTATCCTTGGCAGGACACCCTCTCGGGGCCCGACAGCCTGGCGCAGATGCAATTGAACAACCCCGATGCAA  >  NC_002947/5475172‑5475463

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