Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,173,010	T→A	8.1%	Q165L (CAA→CTA) ‡	PP_1925 ←	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,173,010	0	T	A	8.1%	171.1 / 2.4	60	Q165L (CAA→CTA) ‡	PP_1925	monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base A (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAGC  >  NC_002947/2172862‑2173157                                                                                                                                                     |                                                                                                                                                    ccGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCtt                                                                                                                                                     >  7:142873/1‑149 (MQ=255)       gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAATGACCAGCACGCGCTTTCCTTGGaaa                                                                                                                                                <  2:59621/148‑1 (MQ=255)       gAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGaaa                                                                                                                                                <  1:273541/148‑1 (MQ=255)        aaCAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGc                                                                                                                                              <  1:406365/149‑1 (MQ=255)        aaCAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGAGGTTACCCACGTTGCGTCGGCGACCTCTGACAGCTCAGCATAGATTTGCGCCCCCGAGTAGCAACCACCAACCACCAGCACCCGCTATCCTTGGAAAGc                                                                                                                                              >  7:167614/1‑149 (MQ=255)          cAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgc                                                                                                                                            <  2:300278/149‑3 (MQ=255)              acacGGCCATCAACGGCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGGTGCCACCACCAACGACCAGCACGCGCTTTCCAAGGAAAGcgcgtgc                                                                                                                                        <  1:327679/149‑7 (MQ=255)                  ggCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGaa                                                                                                                                                 <  2:385829/136‑1 (MQ=255)                  ggCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGaa                                                                                                                                                 >  1:385829/1‑136 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATg                                                                                                                            <  8:365868/149‑1 (MQ=255)                          aCGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATg                                                                                                                            >  8:127049/1‑149 (MQ=255)                            gTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATggg                                                                                                                          <  3:164972/149‑1 (MQ=255)                                tcaGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCg                                                                                                                        >  2:317607/1‑147 (MQ=255)                                 caGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAg                                                                                                                     >  4:137436/1‑149 (MQ=255)                                  atgCAGAAACGCCGGTTCCTCCGTTGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCTCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAg                                                                                                                     <  4:63766/146‑1 (MQ=255)                                  aGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGc                                                                                                                         >  6:248489/1‑144 (MQ=255)                                  aGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGc                                                                                                                         <  5:248489/144‑1 (MQ=255)                                  aGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGt                                                                                                                    <  7:127049/149‑1 (MQ=255)                                      aGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTgg                                                                                                                                                   <  6:165061/114‑1 (MQ=255)                                      aGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTgg                                                                                                                                                   >  5:165061/1‑114 (MQ=255)                                         aaCGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATctg                                                                                                             <  6:51752/149‑1 (MQ=255)                                                gTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTccc                                                                                                      <  7:166830/149‑1 (MQ=255)                                                 ttCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCt                                                                                                     <  8:142873/149‑1 (MQ=255)                                                       cGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGaa                                                                                                 <  5:337214/147‑1 (MQ=255)                                                         tGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGc                                                                                             <  7:308119/149‑1 (MQ=255)                                                           gTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGCCCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGctc                                                                                           >  2:297235/1‑149 (MQ=255)                                                           gTTACCCACGTTGCGTCGGCGACCTAGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGACTCGGGTGAGTGATAATGGTCGGAGTGAATCTGCTGTCCCTGGAACAGctc                                                                                           >  2:297239/1‑149 (MQ=255)                                                               cccACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATc                                                                                       >  3:415922/1‑149 (MQ=255)                                                                ccACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATcc                                                                                      >  7:10334/1‑149 (MQ=255)                                                                    gTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCg                                                                                   >  4:420680/1‑148 (MQ=255)                                                                     ttGCGTCGGCGACCTCGTACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGaa                                                                                 <  8:164379/149‑1 (MQ=255)                                                                          tCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGt                                                                                                                    <  8:259713/109‑1 (MQ=255)                                                                          tCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGt                                                                                                                    >  7:259713/1‑109 (MQ=255)                                                                             gtgACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCg                                                                          >  3:127280/3‑148 (MQ=255)                                                                             gCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCgg                                                                         <  3:137436/149‑1 (MQ=255)                                                                                     ggACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGTTTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAg                                                                  <  6:270129/148‑1 (MQ=255)                                                                                           cTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCAGGAACATAGGGGTTGc                                                           >  4:247303/1‑149 (MQ=255)                                                                                            tCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCt                                                          >  1:422537/1‑149 (MQ=255)                                                                                                 aTAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATCGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATACGCCGGAACATAGGGGTTGCTCCAgg                                                     >  2:393418/1‑149 (MQ=255)                                                                                                      tttGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTcccc                                                >  8:346266/1‑149 (MQ=255)                                                                                                        tGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCg                                               <  1:317607/148‑1 (MQ=255)                                                                                                          cgcCCCCGAGTTTCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGtt                                             <  3:209339/148‑1 (MQ=255)                                                                                                             ccccGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTgcgc                                         >  6:372892/1‑149 (MQ=255)                                                                                                             ccccGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTgcgc                                         >  1:101195/1‑149 (MQ=255)                                                                                                             ccacGAGTTGCCACCACTAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGTGGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTgctc                                         <  8:167614/146‑3 (MQ=255)                                                                                                                       ccaccaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGCACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATgacgac                               >  4:244714/1‑149 (MQ=255)                                                                                                                       caaccaccaACAACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGGACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATgacga                                >  7:151823/3‑148 (MQ=255)                                                                                                                        caccaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGa                                                                                  <  1:162091/97‑1 (MQ=255)                                                                                                                        caccaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGa                                                                                  >  2:162091/1‑97 (MQ=255)                                                                                                                        caccaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGaca                              <  6:349524/149‑1 (MQ=255)                                                                                                                        caccaccaACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGaca                              >  5:349524/1‑149 (MQ=255)                                                                                                                                aCGACCAGCACTCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCACTGGAACAGCTATATCCCCGAATAGGACGGAACATAGCGGTTGCTCCAGGTACACGTTGCGCTGATGACGACACGAGCCt                       >  8:289874/1‑148 (MQ=255)                                                                                                                                     cAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGt                 >  8:23751/1‑149 (MQ=255)                                                                                                                                         cgaggcTTTCCAAGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAAAAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTcgcg             >  7:36195/5‑149 (MQ=255)                                                                                                                                         aCGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTcgcg             >  1:376981/1‑149 (MQ=255)                                                                                                                                          taggcTTTCCAATGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                                <  2:170610/94‑1 (MQ=255)                                                                                                                                          gaggcTTTCCAAGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAggg                                                                >  1:170610/4‑97 (MQ=255)                                                                                                                                          ctcgagttCCTTGGAAAGCCTCGGGTGATTTATAATGTGCGGATTGAATCTGCTGTCCCTGGGACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            <  8:151823/143‑1 (MQ=255)                                                                                                                                          cgcgCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCgg            >  8:164862/1‑149 (MQ=255)                                                                                                                                               tttCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCtt       >  3:196593/1‑149 (MQ=255)                                                                                                                                                   cTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAg   >  1:234489/1‑149 (MQ=255)                                                                                                                                                    ttGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAGc  >  5:222432/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCGTTCGAACAGCACACGGCCATCAACGTCATCAGGCAGAAACGCCGGTTCCTCCGTGGTTACCCACGTTGCGTCGGCGACCTCGGACAGCTCCGCATAGATTTGCGCCCCCGAGTTGCCACCACCAACGACCAGCACGCGCTTTCCTTGGAAAGCCTCGGGTGATTGATAATGGGCGGAGTGAATCTGCTGTCCCTGGAACAGCTCTATCCCCGAATAGGCCGGAACATAGGGGTTGCTCCAGGTCCCCGTTGCGCTGATGACGACACGAGCCTCCCAGTCGCGGTCCTTCGAGC  >  NC_002947/2172862‑2173157

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