Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,960,052	A→T	9.1%	noncoding (73/103 nt)	PP_mr63 ←	rnk_pseudo ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,960,052	0	A	T	9.1%	222.9 / 10.9	77	noncoding (73/103 nt)	PP_mr63	rnk_pseudo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base T (5/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

GTACTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTC  >  NC_002947/5959906‑5960197                                                                                                                                                   |                                                                                                                                                  gTACTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAggg                                                                                                                                                >  1:428887/1‑150 (MQ=255)   aCTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCg                                                                                                                                              >  1:800115/1‑150 (MQ=255)     tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                            >  1:451632/1‑150 (MQ=255)     tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                            <  1:1638445/150‑1 (MQ=255)     tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                            <  2:1102912/150‑1 (MQ=255)        cgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGtt                                                                                                                                         <  1:188537/150‑1 (MQ=255)         gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGCCGGAAGAAAGCAAAACCCCGCCCAAGGGCGGGGttt                                                                                                                                        <  1:1893845/150‑1 (MQ=255)          tcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTg                                                                                                                                       <  2:451632/150‑1 (MQ=255)          tcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTg                                                                                                                                       >  1:689333/1‑150 (MQ=255)           cgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGttttg                                                                                                                                      <  1:269265/150‑3 (MQ=255)             aGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTgaga                                                                                                                                    >  2:1319413/1‑150 (MQ=255)              ggCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGaa                                                                                                                                   <  1:1524966/150‑1 (MQ=255)              ggCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGaa                                                                                                                                   <  1:154239/150‑1 (MQ=255)                 tGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGtttt                                                                                                                                       <  1:1519027/143‑2 (MQ=255)                 tGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGtttt                                                                                                                                       >  2:1519027/1‑142 (MQ=255)                  gTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCt                                                                                                                               >  1:360006/1‑150 (MQ=255)                    cGATCAGACGCTCCAGACGCTTTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCGTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCttt                                                                                                                             <  2:1873172/150‑1 (MQ=255)                    cGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCttt                                                                                                                             <  1:818326/150‑1 (MQ=255)                    cGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCttt                                                                                                                             <  1:573396/150‑1 (MQ=255)                     gATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                            <  2:800115/150‑1 (MQ=255)                       tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                          >  1:761818/1‑150 (MQ=255)                       tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                          <  2:360006/150‑1 (MQ=255)                            cGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCAc                                                                                                                     <  2:283017/150‑1 (MQ=255)                                 cAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACg                                                                                                                >  1:1799326/1‑150 (MQ=255)                                  aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGc                                                                                                               <  2:1138496/150‑1 (MQ=255)                                          tACGTCCAATCGGGTGAGGATGAGGGAAGGCTT‑GTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                      >  1:556234/1‑150 (MQ=255)                                              tCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGacac                                                                                                   >  2:1234520/1‑150 (MQ=255)                                                 aaTCGGGTGATGATGAGGGCAGGCTTGGTGCTGATGCTCCAGTCAGACTCCTTTGAACACGCCGTTTTCACTGTGCAAATAAAGTAAAACCCCGCCCAAGGGCGCGGTTTGAGAAGGCTTTGGCTTCACCGCCGCAGAACCCGACACTAc                                                                                                <  2:480308/150‑1 (MQ=255)                                                 aaTCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGCCACTAc                                                                                                >  1:411439/1‑150 (MQ=255)                                                       gTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGc                                                                                          >  1:1497201/1‑150 (MQ=255)                                                          aGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCgg                                                                                       <  1:1319413/150‑1 (MQ=255)                                                              tGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCaa                                                                                   >  1:1527264/1‑150 (MQ=255)                                                               gAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCaaa                                                                                  >  1:415397/1‑150 (MQ=255)                                                                  ggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATAc                                                                               >  2:1868845/1‑150 (MQ=255)                                                                   gAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                      <  1:1502304/126‑1 (MQ=255)                                                                   gAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                      >  2:1502304/1‑126 (MQ=255)                                                                    aaGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTc                                                                             >  2:230104/1‑150 (MQ=255)                                                                    aaGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTc                                                                             >  1:1678271/1‑150 (MQ=255)                                                                      ggCTTGGTGCTCATGACCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAACCAAAACCCCGCCCAAGGGCGGGGTTTGACAAGGCTTTGGCGTCACCGACGCAGAACCCGCCACTACCACAGCCGGGCACACACTCaa                                                                           <  2:862242/150‑1 (MQ=255)                                                                         ttGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGcc                                                                        >  1:621368/1‑150 (MQ=255)                                                                         ttGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGcc                                                                        >  2:1664084/1‑150 (MQ=255)                                                                            gTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCt                                                                     <  1:1378496/150‑1 (MQ=255)                                                                                   tGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCa                                                              <  2:556234/150‑1 (MQ=255)                                                                                    gATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAg                                                             >  2:1488460/1‑150 (MQ=255)                                                                                              gACTCCTTTGAACAGGCAGTTTTAAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAg                                                   >  1:1825507/1‑150 (MQ=255)                                                                                               aCTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGc                                                  <  2:331472/150‑1 (MQ=255)                                                                                                cTCCTTTGATCAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGaca                                                                                                    >  2:395031/1‑99 (MQ=255)                                                                                                cTCCTTTGATCAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGaca                                                                                                    <  1:395031/99‑1 (MQ=255)                                                                                                cTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGcc                                                 >  1:577734/1‑150 (MQ=255)                                                                                                   cTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGc                                              >  1:404177/1‑150 (MQ=255)                                                                                                         aCAGGCAGTTTTCAATGTGCAAATAAAGCAAGAAAACGCATAAGAGAGTAGTGTGAGAAAGCTTTTGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCACATACTCAAGCACCTGCGCTCAGTGTGCCTTAGCCTGCTCgcgc                                        <  2:415397/150‑1 (MQ=255)                                                                                                           aGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCaa                                      <  2:777570/150‑1 (MQ=255)                                                                                                            ggCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAg                                     <  2:1825507/150‑1 (MQ=255)                                                                                                             gCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCg                                                                                        <  1:55625/98‑1 (MQ=255)                                                                                                             gCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCg                                                                                        >  2:55625/1‑98 (MQ=255)                                                                                                             gCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGa                                    >  2:1889330/1‑150 (MQ=255)                                                                                                             gCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACAATACCACAGCCGAGCAAGGACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGa                                    >  1:1870493/1‑150 (MQ=255)                                                                                                               aGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGAcc                                  >  1:1660539/1‑150 (MQ=255)                                                                                                                 ttttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTg                                <  2:761818/150‑1 (MQ=255)                                                                                                                 ttttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGAAACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTg                                >  1:1466783/1‑150 (MQ=255)                                                                                                                        tgtgCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGTGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCAGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcg                         <  2:1176155/150‑1 (MQ=255)                                                                                                                         gtgCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcgc                        >  1:1627037/1‑150 (MQ=255)                                                                                                                           gCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                      <  2:123429/150‑1 (MQ=255)                                                                                                                             aaaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                    >  1:916720/1‑150 (MQ=255)                                                                                                                             aaaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                    <  2:1627037/150‑1 (MQ=255)                                                                                                                                ccttctcaaaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGAt                 >  1:215392/8‑150 (MQ=255)                                                                                                                                 cttctcaaaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATc                >  1:335928/7‑150 (MQ=255)                                                                                                                                  aaGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt               <  2:709868/150‑1 (MQ=255)                                                                                                                                   tctcaaaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTg              >  1:1496608/5‑150 (MQ=255)                                                                                                                                     tcaaaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGcc            >  1:1219629/3‑150 (MQ=255)                                                                                                                                     cAAAACCCCGCCCAAGGGCGGGGTTTGAGCAGGATTTCAAACAAGTGATGTAAAACTCTACACTACCACAGCAGGACAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGcc            <  1:1868845/150‑1 (MQ=255)                                                                                                                                     cAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGcc            >  1:659752/1‑150 (MQ=255)                                                                                                                                          ccccGCCCCAGGGCCGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGCACCCGACACTGCCACAGCCGGTCAAATCCTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCgcgag       <  1:1664084/150‑5 (MQ=255)                                                                                                                                          ccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcgc       >  2:1407654/1‑150 (MQ=255)                                                                                                                                           cccGCCCAAGGGCGGGGAATGACAAGGCATTCGCGTCACCGACGCTGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGcc      >  2:20342/1‑150 (MQ=255)                                                                                                                                             cGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGCTACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGt    >  1:1079346/1‑150 (MQ=255)                                                                                                                                               cccAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTc  <  2:1527264/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTACTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTC  >  NC_002947/5959906‑5960197

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑