Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,960,052	A→T	7.7%	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	7.7%	195.9 / 6.0	65	noncoding (73/103 nt)	PP_mr63	rnk_pseudo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base T (2/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTCG  >  NC_002947/5959909‑5960198                                                                                                                                                 |                                                                                                                                                   ctctCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCgg                                                                                                                                              <  1:782935/150‑1 (MQ=255)  tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                             >  2:1043298/1‑150 (MQ=255)  tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCggg                                                                                                                                             <  2:417182/150‑1 (MQ=255)     cgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGtt                                                                                                                                          <  2:957487/150‑1 (MQ=255)      gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGTTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGttt                                                                                                                                         <  1:566520/150‑1 (MQ=255)       tcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGtttt                                                                                                                                        <  1:255118/150‑2 (MQ=255)        cgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                       <  2:400019/150‑1 (MQ=255)          aGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTgaga                                                                                                                                     <  1:1070994/150‑1 (MQ=255)                  gATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                             <  1:686088/150‑1 (MQ=255)                    tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGCTGAGGGACGGCTTGGTGCTCA‑TGATCCAGTAAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCCAGTGCGCGGTTTGAGAAGGCTTTGGc                                                                                                                           >  2:362503/1‑150 (MQ=255)                    tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                           >  1:955595/1‑150 (MQ=255)                    tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                           <  2:955595/150‑1 (MQ=255)                       gACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTc                                                                                                                        >  2:863715/1‑150 (MQ=255)                           cTCCATGCGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCg                                                                                                                    <  1:362503/150‑1 (MQ=255)                             ccAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGAc                                                                                                                  <  1:71524/150‑1 (MQ=255)                              cAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACg                                                                                                                 >  2:53698/1‑150 (MQ=255)                                        aCGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                       <  2:316303/150‑1 (MQ=255)                                           tCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGCACCCGacac                                                                                                    >  1:126395/1‑150 (MQ=255)                                                   ggTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAg                                                                                            >  2:1053854/1‑150 (MQ=255)                                                   ggTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAg                                                                                            >  1:81709/1‑150 (MQ=255)                                                         gATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                      >  1:54843/1‑150 (MQ=255)                                                         gATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                      >  2:746040/1‑150 (MQ=255)                                                             aGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAAt                                                                                  >  1:594715/1‑150 (MQ=255)                                                              gggAAGGCTTGGTGCTCA‑TGAGCCAGTCAGACTCCGTTGAACAGGCAGTTTTCAATGTGCAAATCCTCATCTACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATa                                                                                 <  1:69297/150‑1 (MQ=255)                                                                gAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCgg                                                                                        >  1:1025361/1‑141 (MQ=255)                                                                gAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCgg                                                                                        <  2:1025361/141‑1 (MQ=255)                                                                 aaGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGAc                                                                                                                  >  2:1000981/1‑114 (MQ=255)                                                                 aaGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGAc                                                                                                                  <  1:1000981/114‑1 (MQ=255)                                                                  aGGCTTGGTGATCACTGCT‑CAGGCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCa                                                                             <  2:758255/150‑1 (MQ=255)                                                                    gCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAg                                                                           >  1:89328/1‑150 (MQ=255)                                                                            cTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTgcg                                                                   >  2:120086/1‑150 (MQ=255)                                                                             tCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                      <  2:1016497/130‑1 (MQ=255)                                                                             tCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                      >  1:1016497/1‑130 (MQ=255)                                                                                    tCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACCGCCGGGCAAATACTCAAGCCCCTGCGCTCAgtg                                                            <  1:413221/150‑1 (MQ=255)                                                                                    tCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCACTGCGCTCAgtg                                                            <  1:746040/150‑1 (MQ=255)                                                                                     ccAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAgtgt                                                           <  2:1047745/150‑1 (MQ=255)                                                                                      cAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGCACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGtgtg                                                          >  2:579002/1‑150 (MQ=255)                                                                                       aGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGCACCCCCCACTACCACAGCCGTTCAAATACTCAAGCCCCTGCGCTCACTGTGc                                                         >  2:851527/1‑150 (MQ=255)                                                                                           aGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATACTGGAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCACTGTGCCTTa                                                     <  1:851527/150‑1 (MQ=255)                                                                                              cTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGcc                                                  <  1:611195/150‑1 (MQ=255)                                                                                                 cTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGc                                               >  1:176507/1‑150 (MQ=255)                                                                                                        cAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCa                                        >  2:846697/1‑150 (MQ=255)                                                                                                          ggCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACAACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAg                                      >  2:1113396/1‑150 (MQ=255)                                                                                                           gCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGc                                                   <  2:835152/136‑1 (MQ=255)                                                                                                           gCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGc                                                   >  1:835152/1‑136 (MQ=255)                                                                                                             aGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTGTGAGAAGGCTTTTGCGTCACCGACGCAGAACCCGACACTACCAGAGCCGGGCAAATGCGCTAGCCCCTGCGCTCAGTGGGCCTTAGGCGGCTGGCGCAAGccc                                   >  2:505006/1‑150 (MQ=255)                                                                                                                tttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTgg                                <  1:811771/150‑1 (MQ=255)                                                                                                                tttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTgg                                <  1:552017/150‑1 (MQ=255)                                                                                                                tttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTgg                                >  1:823253/1‑150 (MQ=255)                                                                                                                   cAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGcc                             <  2:54843/150‑1 (MQ=255)                                                                                                                         gCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                       >  1:319355/1‑150 (MQ=255)                                                                                                                           aaaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                     >  2:432985/1‑150 (MQ=255)                                                                                                                           aaaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                     >  1:648242/1‑150 (MQ=255)                                                                                                                              ccttctcaaaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGAt                  >  2:233043/8‑150 (MQ=255)                                                                                                                                    aaaaCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCg            <  1:120086/150‑1 (MQ=255)                                                                                                                                       aCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcg         >  1:470784/1‑150 (MQ=255)                                                                                                                                       aCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                       >  2:688727/1‑56 (MQ=255)                                                                                                                                       aCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                       <  1:688727/56‑1 (MQ=255)                                                                                                                                       aCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcg         >  1:214050/1‑150 (MQ=255)                                                                                                                                        ccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcgc        <  1:1043298/150‑1 (MQ=255)                                                                                                                                         cccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGTAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGcc       >  1:113799/1‑150 (MQ=255)                                                                                                                                         cccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCACGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGcc       >  1:1082253/1‑150 (MQ=255)                                                                                                                                          ccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCg      <  2:594715/150‑1 (MQ=255)                                                                                                                                          ccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCg      >  1:807040/1‑150 (MQ=255)                                                                                                                                              ccAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTCg  >  1:575431/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CTCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCA‑TGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTCG  >  NC_002947/5959909‑5960198

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