Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,315,516	Δ1 bp	6.0%	coding (255/501 nt)	yciF ←	putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,315,516	0	A	.	6.0%	214.2 / 6.8	67	coding (255/501 nt)	yciF	putative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base . (2/2);  total (34/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.09e-01

CGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCA  >  NC_000913/1315372‑1315660                                                                                                                                                 |                                                                                                                                                 cGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGACATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGc                                                                                                                                                <  1:370884/147‑1 (MQ=255)         tACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAAt                                                                                                                                       >  1:321840/1‑148 (MQ=255)          aCGGTAACCTAATTGTTCAGCCAGCGTCGCCAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATa                                                                                                                                      <  8:122173/148‑1 (MQ=255)            ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGcttctt                                                                                                                                           >  6:264668/1‑141 (MQ=255)            ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGcttctt                                                                                                                                           <  5:264668/141‑1 (MQ=255)            ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGa                                                                                                                                   <  4:343516/149‑1 (MQ=255)                  cTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCCCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTc                                                                                                                              >  5:244842/1‑148 (MQ=255)                     aTTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCG‑TTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTg                                                                                                                         >  5:19819/1‑149 (MQ=255)                            agccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcaca                                                                                                                    >  5:70856/1‑148 (MQ=255)                            agccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACGTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcaca                                                                                                                    >  8:92274/1‑148 (MQ=255)                                   gTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCAt                                                                                                            <  7:114992/149‑1 (MQ=255)                                     cGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATgc                                                                                                          >  3:92937/1‑149 (MQ=255)                                              cccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATT‑GCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                                                                >  1:373858/1‑149 (MQ=255)                                              cccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATT‑GCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                                                                >  1:379333/1‑149 (MQ=255)                                                cGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTGCTTCAATAAGCCCTTCCATTGCCACACATTTCATGCGCTTAATTTTc                                                                                               <  3:74575/149‑1 (MQ=255)                                                 gTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                                                                >  5:360078/1‑147 (MQ=255)                                                        gCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                                                       >  1:191203/1‑149 (MQ=255)                                                          aaTCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGAtt                                                                                     <  6:70856/149‑1 (MQ=255)                                                            tctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                                                   >  7:261548/1‑149 (MQ=255)                                                            tctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCAGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                                                   >  1:208520/1‑149 (MQ=255)                                                              tcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCAATTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTGCTTCAATAAGACGTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGa                                                                                 >  4:18425/1‑149 (MQ=255)                                                                aTAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                                                   <  2:321840/145‑1 (MQ=255)                                                                aTAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACt                                                                               >  8:299542/1‑149 (MQ=255)                                                                 tAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCAGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTc                                                                              <  2:208520/149‑1 (MQ=255)                                                                       tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                                                       >  2:146930/1‑134 (MQ=255)                                                                       tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                                                       <  1:146930/134‑1 (MQ=255)                                                                       tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTAGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATATCATGCGCTTAATTTTCAGTTTCGATTCGGACTCCACAAc                                                                        >  8:126476/1‑149 (MQ=255)                                                                        cGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                                                        >  5:307136/1‑148 (MQ=255)                                                                        cGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                                                        <  6:307136/148‑1 (MQ=255)                                                                         gACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACtt                                                                      <  8:97461/149‑1 (MQ=255)                                                                              tCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGAt                                                                   <  7:299542/147‑1 (MQ=255)                                                                                tgtgCTGTGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCaa                                                                >  6:79555/1‑148 (MQ=255)                                                                                  tgctgcGGCAATCCGTCCGGCATCACGCCCTTCGTGTTGCGCGGTACTATCGAGGACCTCATACGCTTCTTCAATCAGACCTTCCCTTGCCACACATTTCATCCGCATAATATTCACGCTCGATTCGGACTCCACAACTTGATCAATAc                                                             <  6:301948/149‑1 (MQ=255)                                                                                    ctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                                                                  <  6:45614/142‑1 (MQ=255)                                                                                    ctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                                                                  >  5:45614/1‑142 (MQ=255)                                                                                        ggCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCATCCATTGCCACACATTTCATGCGATTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCaa                                                       <  7:92274/149‑1 (MQ=255)                                                                                         gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCTGACTCCACAACTTGATCAATACGTTCagt                                                      >  7:342632/1‑147 (MQ=255)                                                                                          cAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                                                     >  2:140566/1‑149 (MQ=255)                                                                                             tCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTTAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACTACTTGATCAATACGTTCAATCTgg                                                  <  1:110994/149‑2 (MQ=255)                                                                                               aGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTcc                                                <  6:244842/149‑1 (MQ=255)                                                                                                   cGGCATCACTCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATg                                             <  7:255030/148‑1 (MQ=255)                                                                                                   cGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGa                                            <  3:241516/149‑1 (MQ=255)                                                                                                   cGGCATCACGCACTTCG‑TTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGa                                            <  6:19819/148‑1 (MQ=255)                                                                                                    ggCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCAt                                              >  4:316342/1‑146 (MQ=255)                                                                                                    ggCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGa                                            <  3:316342/148‑1 (MQ=255)                                                                                                       aTCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATT‑GCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTc                                       <  2:379333/149‑1 (MQ=255)                                                                                                       aTCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATT‑GCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTc                                       <  2:373858/149‑1 (MQ=255)                                                                                                          ccGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAATACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGAtt                                                                                     <  2:200918/100‑1 (MQ=255)                                                                                                          aCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCt                                     <  1:204184/149‑1 (MQ=255)                                                                                                            gCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTc                                    >  5:279696/1‑148 (MQ=255)                                                                                                              aCTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAg                                 >  2:259351/1‑149 (MQ=255)                                                                                                                ttCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGt                               >  1:174709/1‑149 (MQ=255)                                                                                                                ttCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGt                               >  2:302922/1‑149 (MQ=255)                                                                                                                   gTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTGCGATTCGCACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTgcg                            >  8:320451/1‑149 (MQ=255)                                                                                                                      tttCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCa                          <  1:259351/148‑1 (MQ=255)                                                                                                                          tcGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATTCGCTTAATTTTCAGGTTCGATTCGGCCTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGaa                      >  5:167996/1‑148 (MQ=255)                                                                                                                          tcGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGaaa                     >  2:92275/1‑149 (MQ=255)                                                                                                                            ggTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAg                   >  1:253782/1‑149 (MQ=255)                                                                                                                             gTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGc                  <  1:302922/149‑1 (MQ=255)                                                                                                                                cTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCt                <  6:360078/148‑1 (MQ=255)                                                                                                                                cTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCt                <  8:261548/148‑1 (MQ=255)                                                                                                                                     gATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACtt           >  2:46173/1‑148 (MQ=255)                                                                                                                                       tGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAAt        <  2:174709/149‑1 (MQ=255)                                                                                                                                           cTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAAttttt    >  3:288411/1‑149 (MQ=255)                                                                                                                                           cTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAAttttt    <  4:92937/149‑1 (MQ=255)                                                                                                                                             cATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCa  >  2:324325/1‑149 (MQ=255)                                                                                                                                              aTTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTc   <  1:140566/147‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCA  >  NC_000913/1315372‑1315660

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