Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,612,119	A→T	8.0%	I522N (ATT→AAT)	mcoA ←	Mn(II) copper oxidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,612,119	0	A	T	8.0%	180.1 / 4.6	62	I522N (ATT→AAT)	mcoA	Mn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (2/3);  total (31/31)
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.89e-01

GGACGGGCCTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAGGTTGAC  >  NC_002947/3611971‑3612267                                                                                                                                                     |                                                                                                                                                     ggACGGGCCTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGaa                                                                                                                                                      >  7:208690/1‑149 (MQ=255)   aCGGGCCTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATc                                                                                                                                                    >  3:31959/1‑149 (MQ=255)       gCCTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAAGTTTAAGTAGCGCGACACCGAGCCGTTGATGAGACGGAAACGGTAACTGCGCGCAGGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGt                                                                                                                                                <  1:266014/149‑1 (MQ=255)       gCCTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGt                                                                                                                                                <  2:327804/149‑1 (MQ=255)        ccTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGt                                                                                                                                                >  8:232446/1‑148 (MQ=255)          tttGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAg                                                                                                                                             >  4:167961/1‑149 (MQ=255)            tGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGATTCTGGTCaac                                                                                                                                            <  4:195364/148‑3 (MQ=255)            tGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGc                                                                                                                                            <  7:232446/148‑1 (MQ=255)            tGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGc                                                                                                                                            <  1:6968/148‑1 (MQ=255)             gAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGATTCTGGTCaact                                                                                                                                           <  6:240013/148‑4 (MQ=255)             gAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGc                                                                                                                                            >  5:100561/1‑147 (MQ=255)               aCTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAgg                                                                                                                                        <  2:157110/149‑1 (MQ=255)                 tCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAgg                                                                                                                                        >  1:366505/1‑147 (MQ=255)                   aCCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAg                                                                                                                                     <  4:31959/148‑1 (MQ=255)                   aCCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGc                                                                                                                                    >  4:46193/1‑149 (MQ=255)                             ccGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGtt                                                                                                                          <  8:7214/149‑1 (MQ=255)                              cGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGtt                                                                                                                          >  2:227440/1‑148 (MQ=255)                                      tcCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGtt                                                                                                                 >  1:231565/1‑149 (MQ=255)                                        ccGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGt                                                                                                                  <  2:175228/146‑1 (MQ=255)                                          gCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAAc                                                                                                             <  8:208690/149‑1 (MQ=255)                                              cacGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGTCCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACcaca                                                                                                         <  1:279010/149‑1 (MQ=255)                                              cacGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACcaca                                                                                                         >  7:30737/1‑149 (MQ=255)                                                    gAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGc                                                                                                   >  8:170644/1‑149 (MQ=255)                                                       cTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGTCGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAg                                                                                                                                         <  5:185365/108‑1 (MQ=255)                                                       cTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGGCGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAg                                                                                                                                         >  6:185365/1‑108 (MQ=255)                                                        ttGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCAtt                                                                                               >  8:96461/1‑149 (MQ=255)                                                          gAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTg                                                                                              <  1:226183/148‑1 (MQ=255)                                                          gAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTg                                                                                              <  5:258872/148‑1 (MQ=255)                                                          gAAGTAGCGCGACACCGAACCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGcc                                                                                                                                   >  1:324481/1‑111 (MQ=255)                                                          gAAGTAGCGCGACACCGAACCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGcc                                                                                                                                   <  2:324481/111‑1 (MQ=255)                                                                    gACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGc                                                                                   <  8:30737/149‑1 (MQ=255)                                                                             ccGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCgg                                                                           >  5:258093/1‑148 (MQ=255)                                                                             ccGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGc                                                                          <  3:139491/149‑1 (MQ=255)                                                                                ttGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAAt                                                                       <  6:258093/149‑1 (MQ=255)                                                                                   aggagACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCAc                                                                    <  1:227440/149‑1 (MQ=255)                                                                                     gagaCGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCa                                                                  <  7:91944/149‑1 (MQ=255)                                                                                         cGGAAACGGTAACTGCGCGCACGCACTTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGtt                                                              <  6:100561/149‑1 (MQ=255)                                                                                          ggAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGtt                                                              >  3:181752/1‑148 (MQ=255)                                                                                          ggAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGtt                                                              >  3:370089/1‑148 (MQ=255)                                                                                                          cgcacgcaCCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGt                                             <  7:170644/149‑1 (MQ=255)                                                                                                           gcacgcacGTTGAGCCTGGGCGTGAACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGATGAACGGGTTGCACCACAGCTGGCCATTGGCATCCGAGGCTTTGTCGGCAATCACCAGGTTGACGTCATACTCACGGt                                             >  5:66447/1‑148 (MQ=255)                                                                                                             acgcacCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCACTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTAc                                          >  1:319588/1‑149 (MQ=255)                                                                                                               gcacCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCa                                                       <  8:269985/134‑1 (MQ=255)                                                                                                               gcacCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCa                                                       >  7:269985/1‑134 (MQ=255)                                                                                                                  ccTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCaa                                     <  4:181752/149‑1 (MQ=255)                                                                                                                    ttGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCaa                                     >  3:261351/1‑147 (MQ=255)                                                                                                                    ttGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCaa                                     <  4:261351/147‑1 (MQ=255)                                                                                                                         ccTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATgg                              >  5:219366/1‑149 (MQ=255)                                                                                                                               cTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTAc                                          >  5:43491/1‑131 (MQ=255)                                                                                                                               cTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTAc                                          <  6:43491/131‑1 (MQ=255)                                                                                                                                  gTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTAc                                          <  6:302217/128‑1 (MQ=255)                                                                                                                                  gTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTAc                                          >  5:302217/1‑128 (MQ=255)                                                                                                                                     cTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTg                  <  2:366505/149‑1 (MQ=255)                                                                                                                                      tGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCGCAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTgg                 >  1:107964/1‑149 (MQ=255)                                                                                                                                          aGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAg             >  1:190854/1‑149 (MQ=255)                                                                                                                                           gctGACCAGATTCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAgc            >  3:366973/3‑149 (MQ=255)                                                                                                                                            ctGACCAGATTCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACgg                                              <  7:70041/113‑1 (MQ=255)                                                                                                                                            ctGACCAGATTCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACgg                                              >  8:70041/2‑114 (MQ=255)                                                                                                                                              gACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCTCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAgcg           >  2:151568/1‑147 (MQ=255)                                                                                                                                               aCCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAg        >  1:238169/1‑149 (MQ=255)                                                                                                                                               aCCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAg        <  3:167961/149‑1 (MQ=255)                                                                                                                                                ccAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAgg       >  5:26910/1‑149 (MQ=255)                                                                                                                                                 cAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAGGt      >  6:35348/1‑149 (MQ=255)                                                                                                                                                     aTCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCTGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAGGTTGAc  <  7:96461/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGACGGGCCTTTGAACTCACCGCTGTTACCGGCAATCTCCCGCACCACGGCGAACTTGAAGTAGCGCGACACCGAGCCGTTGAGGAGACGGAAACGGTAACTGCGCGCACGCACCTTGAGCCTGGGCTTGTACTGCCAGTTGACCAGAATCTGGTCAGCGAGGAAGCCGTCCGTGTTGAACGGGTTGAACCACAGCTGGCCATTGGCATCCCAGGCTTTGTCGGCAATCACCAGGTTGACGTCATAGTCACGGTTACCCCAAGGCATGGCTGAACCGCTTGGGAAGCGCAGGTTGAC  >  NC_002947/3611971‑3612267

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