Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,959,995	A→T	8.7%	noncoding (74/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,959,995	0	A	T	8.7%	175.2 / 2.1	63	noncoding (74/103 nt)	PP_mr63	rnk_pseudo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (2/4);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTCGT  >  NC_002947/5959854‑5960143                                                                                                                                              |                                                                                                                                                     tctcGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGTATAGGCAGTTTTCTATGTGCAAATAAAGC‑AAACCCCGCCCTTCGGCggg                                                                                                                                              >  6:92443/1‑149 (MQ=255)   tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCgggg                                                                                                                                             >  8:32070/1‑149 (MQ=255)   tcgtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCgggg                                                                                                                                             >  8:262388/1‑149 (MQ=255)     gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGTGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCTGTCAGTCTCGTTTTAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCgggg                                                                                                                                             <  3:215717/147‑1 (MQ=255)     gtcgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGtt                                                                                                                                           <  3:327300/149‑1 (MQ=255)        gAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGCTCCAGTCAGACTCCTTTGAATAGGCAGTTTTCAATGTGCAAATAAAGC‑AAACCCCGCCCTTGGGCGGGGttttgc                                                                                                                                       <  5:92443/149‑4 (MQ=255)        gAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGtttt                                                                                                                                         <  3:57167/148‑2 (MQ=255)           gCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCGGGGttttgctt                                                                                                                                     <  1:187253/149‑6 (MQ=255)             tGTCGATCAGACGCTCCAGACGATGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAgg                                                                                                                                   >  5:23182/1‑149 (MQ=255)               tCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCt                                                                                                                                 <  7:262388/149‑1 (MQ=255)                 gATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCttt                                                                                                                               >  8:40905/1‑149 (MQ=255)                  aTCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                              <  2:246379/149‑1 (MQ=255)                     aGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCg                                                                                                                           >  3:157589/1‑149 (MQ=255)                      gACGCTCCAGACGGTGTAGGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTg                                                                                                                                         <  2:129770/134‑1 (MQ=255)                           tCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCAcc                                                                                                                      >  6:320956/1‑148 (MQ=255)                             cAGACGCTGTACGTCCAATCGGGTGAGGATGAGTGAAGGCTTTGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCATTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGAc                                                                                                                   <  4:80668/149‑1 (MQ=255)                              aGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACg                                                                                                                  >  1:191418/1‑149 (MQ=255)                                      tACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAAccc                                                                                                          >  2:330720/1‑149 (MQ=255)                                        cGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                        <  5:40616/149‑1 (MQ=255)                                        cGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                        <  5:320956/149‑1 (MQ=255)                                         gTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                        >  3:95156/1‑148 (MQ=255)                                         gTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGa                                                                                                        <  4:95156/148‑1 (MQ=255)                                            cAATCGGGTGAGGATGAGGGAAGGCTTTGTGCTCATGATCCAGTCTGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGGGTCACCGACGCAGATCCCGtcagt                                                                                                    <  8:108860/149‑3 (MQ=255)                                            cAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACt                                                                                                    <  7:40905/149‑1 (MQ=255)                                             aaTCGGGTGAGGATGAGGGACGGCTGGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTa                                                                                                   <  8:132472/149‑1 (MQ=255)                                              aTCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTAc                                                                                                  <  6:131335/149‑1 (MQ=255)                                               tCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTAc                                                                                                  >  1:78184/1‑148 (MQ=255)                                                        gATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCggg                                                                                        >  7:286819/1‑149 (MQ=255)                                                        gATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCggg                                                                                        >  4:289720/1‑149 (MQ=255)                                                             gggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCggg                                                                                        >  6:212653/1‑144 (MQ=255)                                                             gggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCggg                                                                                        <  5:212653/144‑1 (MQ=255)                                                             gggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCaaa                                                                                    >  8:22155/1‑148 (MQ=255)                                                              ggAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATa                                                                                  >  1:20812/1‑149 (MQ=255)                                                                 cgccgtGGTCCTCCTGCCCCAGCCACCCTCCTTTGCACAGCCCGTTTTCAATGCCCAAATCACCCAAACCCCCGCCCACGGGCGGGGTTTGAGAAGGCGTTGGCGCCCCCGACGCAGAACCCGACACTACCACAGCCGCGCAACTACTc                                                                               <  4:337665/144‑1 (MQ=255)                                                                 aGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTc                                                                               <  2:212349/149‑1 (MQ=255)                                                                   gCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCaa                                                                             <  8:223089/149‑1 (MQ=255)                                                                   gCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCaa                                                                             <  2:20812/149‑1 (MQ=255)                                                                    cTTGGTGCTCATGATCCTGTCAGACTCCTTTGAACTGGCAGTTTTCATTGTGCAAATAAAGCATAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCaa                                                                             <  6:23182/148‑1 (MQ=255)                                                                    cTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCaa                                                                             <  7:22155/148‑1 (MQ=255)                                                                    cTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAg                                                                            <  2:191418/149‑1 (MQ=255)                                                                        gTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGcccc                                                                        >  3:177661/1‑149 (MQ=255)                                                                                 aTCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCa                                                                <  2:78184/148‑1 (MQ=255)                                                                                 aTCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAg                                                               <  1:330720/149‑1 (MQ=255)                                                                                 aTCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAg                                                               >  8:253840/1‑149 (MQ=255)                                                                                  tCCAGTCAGACTCCTTTGAACAGGCAGTTTTCGATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                       <  2:36246/124‑1 (MQ=255)                                                                                  tCCAGTCAGACTCCTTTGAACAGGCAGTTTTCGATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                       >  1:36246/1‑124 (MQ=255)                                                                                     aGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAACTACTCAAGCCCCTGCGCTCAGtgtg                                                           >  4:68740/1‑149 (MQ=255)                                                                                      gTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGc                                                          <  4:157589/149‑1 (MQ=255)                                                                                          gACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCt                                                        >  6:73778/1‑147 (MQ=255)                                                                                            cTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAgt                                                              >  7:258008/1‑139 (MQ=255)                                                                                            cTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAgt                                                              <  8:258008/139‑1 (MQ=255)                                                                                                    aaCAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTcgc                                            >  4:78474/1‑149 (MQ=255)                                                                                                       aGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGCCACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCgcgc                                          >  7:137076/1‑148 (MQ=255)                                                                                                        ggCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCaa                                        <  7:32070/149‑1 (MQ=255)                                                                                                           aGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGAc                                     <  6:276096/149‑1 (MQ=255)                                                                                                                       gCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcgcg                         >  3:103990/1‑149 (MQ=255)                                                                                                                          aaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                        >  4:49702/1‑147 (MQ=255)                                                                                                                          aaTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCa                      <  7:253840/149‑1 (MQ=255)                                                                                                                              aaGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGAt                   <  3:49702/148‑1 (MQ=255)                                                                                                                               tctcaaaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCAGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt                 >  8:325812/5‑149 (MQ=255)                                                                                                                               aGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt                 <  3:289720/149‑1 (MQ=255)                                                                                                                                           cccAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGc                                                          <  3:126375/96‑1 (MQ=255)                                                                                                                                           cccAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGc                                                          >  4:126375/1‑96 (MQ=255)                                                                                                                                              aaGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTCGt  >  2:177994/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     TCTCGTCGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGTTCGT  >  NC_002947/5959854‑5960143

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