Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,039,653	A→C	7.5%	I284S (ATC→AGC) ‡	PP_4446 ←	group II intron‑encoding maturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,039,653	0	A	C	7.5%	191.3 / 4.5	64	I284S (ATC→AGC) ‡	PP_4446	group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base C (3/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAACGTC  >  NC_002947/5039505‑5039799                                                                                                                                                     |                                                                                                                                                    aTTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGa                                                                                                                                                     <  5:260232/149‑1 (MQ=255) aTTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGa                                                                                                                                                     <  3:174611/149‑1 (MQ=255)  ttACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGAt                                                                                                                                                    <  4:172054/149‑1 (MQ=255)          gCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCAt                                                                                                                                            <  2:137114/149‑1 (MQ=255)           cTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCag                                                                                                                                            >  3:279058/1‑147 (MQ=255)               gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACGACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGt                                                                                                                                       >  8:43760/1‑149 (MQ=255)               gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGt                                                                                                                                                >  3:138300/1‑140 (MQ=255)               gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGt                                                                                                                                                <  4:138300/140‑1 (MQ=255)                 cgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCg                                                                                                                                     <  7:128308/149‑1 (MQ=255)                       aGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGttt                                                                                                                               <  7:310016/149‑1 (MQ=255)                              gccgccATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTc                                                                                                                         <  3:97502/148‑1 (MQ=255)                                    aTAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGttt                                                                                                                               <  8:33936/136‑1 (MQ=255)                                    aTAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGttt                                                                                                                               >  7:33936/1‑136 (MQ=255)                                      aaaCGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCa                                                                                                                             <  6:246403/136‑1 (MQ=255)                                      aaaCGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACAGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCa                                                                                                                             >  5:246403/1‑136 (MQ=255)                                       aaCGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCAt                                                                                                                <  1:104436/148‑1 (MQ=255)                                         cGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGgt                                                                                                                                   >  3:98537/1‑127 (MQ=255)                                         cGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGgt                                                                                                                                   <  4:98537/127‑1 (MQ=255)                                          ggCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATg                                                                                                               >  1:109697/1‑146 (MQ=255)                                          ggCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATg                                                                                                               <  2:109697/146‑1 (MQ=255)                                          ggCAGACCGCCGAACGAAAACCACACAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACGAAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGc                                                                                                             >  3:276085/1‑148 (MQ=255)                                            cAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTTGTTCCGCTCATCCCGCAGGGCcag                                                                                                          >  6:251147/1‑149 (MQ=255)                                             aGACCGCCGAACGAAAACCACGCAGACGTGTCAGGTTTCCCGTCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCATGAACCCATCTTGGTTCaag                                                                                                                                           <  8:153394/115‑3 (MQ=255)                                             aGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCcagc                                                                                                         <  4:328511/149‑1 (MQ=255)                                             aGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCATGAACCCATCTGGGTTCaag                                                                                                                                           <  8:153392/115‑3 (MQ=255)                                             aGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCATGAACCCATCTGGGTTCaag                                                                                                                                           >  7:153392/1‑113 (MQ=255)                                             aGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTATCCCGGCACCGCGAGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATAGCCCATGAACCCATCTCGGTTCaac                                                                                                                                           >  7:153394/1‑113 (MQ=255)                                                  gacgaacgaaAACCACACAGACGTGTCAGGTTTCCAGCCACCTCGTGGAAATTGAATTACAAGCTCACACCCCTTGTGATCCATTGCCATTCAACCCAGTAGGGTTCATGACGTCGGTGTTTCAGCTAATCCCGCATGCACAGCAGCAt                                                                                                    >  4:307882/3‑149 (MQ=255)                                                   ccgAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATc                                                                                                   <  2:304890/149‑1 (MQ=255)                                                     gaacgaaAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATcg                                                                                                  >  8:110175/1‑148 (MQ=255)                                                             accaccCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTc                                                                                         <  3:227354/149‑1 (MQ=255)                                                               caccCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑Gc                                                                                      >  3:124047/1‑149 (MQ=255)                                                                   cAGACGTCTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑Gcttct                                                                                  <  6:120233/149‑1 (MQ=255)                                                                    aGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑Gcttctt                                                                                 <  4:124047/149‑1 (MQ=255)                                                                                tttCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCg                                                                     >  2:203294/1‑149 (MQ=255)                                                                                          ccGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCTGCTTCTTGACCGTCAGTCGCAGTATTTg                                                            <  1:203294/149‑1 (MQ=255)                                                                                           cgcgTGGTAATTGAAGTAACCGCTGACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGa                                                           >  2:338743/1‑148 (MQ=255)                                                                                             cgTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAAc                                                        <  7:123976/149‑1 (MQ=255)                                                                                              gTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAAcc                                                       >  4:93133/1‑149 (MQ=255)                                                                                                 gTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCAc                                                    >  1:38161/1‑149 (MQ=255)                                                                                                 gTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTAGCAGTATTTGAAACCCAc                                                    >  8:38215/1‑149 (MQ=255)                                                                                                  tAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCAcc                                                   <  7:38215/149‑1 (MQ=255)                                                                                                    aTTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCg                                                  <  5:251147/148‑1 (MQ=255)                                                                                                        aaGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCa                                                                   <  6:183186/127‑1 (MQ=255)                                                                                                        aaGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCa                                                                   >  5:183186/1‑127 (MQ=255)                                                                                                          gTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGt                                           >  3:301161/1‑149 (MQ=255)                                                                                                           tAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGtt                                          <  4:301161/149‑1 (MQ=255)                                                                                                             aCCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTgg                                        <  2:38161/149‑1 (MQ=255)                                                                                                              ccGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTgg                                        >  7:329442/1‑148 (MQ=255)                                                                                                               cGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGta                                      <  8:329442/149‑1 (MQ=255)                                                                                                                       accCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTacaaca                              <  7:110175/149‑1 (MQ=255)                                                                                                                           tggtgAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGc                          <  6:90167/147‑1 (MQ=255)                                                                                                                                  ccATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGc                   >  3:160943/1‑149 (MQ=255)                                                                                                                                  ccATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACAGTCAGTCGCAGTATTTGAAACCCACGGCTTCTGTTTGTACTACAACAGTGCGTGAAGc                   >  4:9699/1‑149 (MQ=255)                                                                                                                                    aTTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGcc                  >  6:34412/1‑148 (MQ=255)                                                                                                                                        cccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAgg              >  2:16320/1‑148 (MQ=255)                                                                                                                                         ccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGa                                                                               >  1:78153/1‑82 (MQ=255)                                                                                                                                         ccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGa                                                                               <  2:78153/82‑1 (MQ=255)                                                                                                                                           atGAACCCATCTGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAACCCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAg           >  3:254342/2‑148 (MQ=255)                                                                                                                                             gAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGt          <  2:135469/147‑1 (MQ=255)                                                                                                                                             gAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTc         <  7:43760/148‑1 (MQ=255)                                                                                                                                             gAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTc         <  1:338743/148‑1 (MQ=255)                                                                                                                                              aaCCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGGTGGGTCTACGACGGTGCGGGAAGCCCAGGAAGTCa        >  8:280815/1‑148 (MQ=255)                                                                                                                                                cccAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGc                          <  5:99842/128‑1 (MQ=255)                                                                                                                                                cccAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGc                          >  6:99842/1‑128 (MQ=255)                                                                                                                                                   aGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAACGtc  >  5:147117/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTC‑GCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAACGTC  >  NC_002947/5039505‑5039799

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