Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,367,609	A→T	6.1%	H44L (CAT→CTT)	PP_4720 →	RNA‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,367,609	0	A	T	6.1%	199.4 / 3.0	65	H44L (CAT→CTT)	PP_4720	RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base T (2/2);  total (36/29)
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.06e-01

GACGGAAAGGGTCAGATTATGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAA  >  NC_002947/5367461‑5367756                                                                                                                                                     |                                                                                                                                                     gACGGAAAGGGTCAGATTATGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCa                                                                                                                                                      <  4:330924/149‑1 (MQ=255)    ggAAAGGGTCAGATTATGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGa                                                                                                                                                   >  5:44394/1‑149 (MQ=255)            tCAGATTATGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATc                                                                                                                                            >  5:368842/1‑148 (MQ=255)              aGATTATGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAg                                                                                                                                         <  3:290356/149‑1 (MQ=255)                    tGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTCAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCTTGATCAGCTCATGGTCgac                                                                                                                                   <  5:368852/149‑2 (MQ=255)                     gCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTc                                                                                                                                      >  6:25559/1‑145 (MQ=255)                       cGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTc                                                                                                                                      >  1:188988/1‑143 (MQ=255)                       cGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTc                                                                                                                                      <  2:188988/143‑1 (MQ=255)                       cGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGAtt                                                                                                                                >  4:175875/1‑149 (MQ=255)                        gCTCAATAACGAGCAGAAGAAGCAATTCATGTCCTTTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGGTGATCAAGGTCGAGATTc                                                                                                                               <  6:327726/149‑1 (MQ=255)                        gCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGGCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTc                                                                                                                               <  4:384723/149‑1 (MQ=255)                                         agaagCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCg‑aa                                                                                                             <  2:320378/149‑1 (MQ=255)                                         agaagCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCg‑aa                                                                                                             <  4:336501/149‑1 (MQ=255)                                             gCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACcgc                                                                                                                 <  5:306507/142‑1 (MQ=255)                                             gCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACcgc                                                                                                                 >  6:306507/1‑142 (MQ=255)                                               aaTACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGGAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCGCTGGACGTCCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGCCCGCG‑AAGAAcgc                                                                                                       >  2:292758/1‑149 (MQ=255)                                                 tACAAGTCAATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACgcgc                                                                                                     <  6:44394/149‑1 (MQ=255)                                                      gTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCgaga                                                                                                <  6:368842/149‑1 (MQ=255)                                                        ccATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGAc                                                                                               >  7:133738/1‑148 (MQ=255)                                                         cATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCa                                                                                             >  1:241547/1‑149 (MQ=255)                                                          aTTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAAcgc                                                                                                       >  1:306629/1‑138 (MQ=255)                                                          aTTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAAcgc                                                                                                       <  2:306629/138‑1 (MQ=255)                                                              gTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGc                                                                                         <  3:175875/148‑1 (MQ=255)                                                                  tGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCGCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGc                                                                                    <  4:354873/149‑1 (MQ=255)                                                                     ccTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAgg                                                                                                                                        >  1:249917/1‑95 (MQ=255)                                                                     ccTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAgg                                                                                                                                        <  2:249917/95‑1 (MQ=255)                                                                      cTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCt                                                                                                                            <  3:207893/106‑1 (MQ=255)                                                                      cTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCt                                                                                                                            >  4:207893/1‑106 (MQ=255)                                                                      cTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACgcgc                                                                                                     >  7:251134/1‑128 (MQ=255)                                                                      cTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACgcgc                                                                                                     <  8:251134/128‑1 (MQ=255)                                                                      cTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCACTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCtgtg                                                                                >  4:126907/1‑149 (MQ=255)                                                                        gAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGc                                                                               >  1:171684/1‑148 (MQ=255)                                                                           gCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAgg                                                                           >  7:308123/1‑149 (MQ=255)                                                                                    gATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTg                                                                  >  2:106967/1‑149 (MQ=255)                                                                                      tCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGcc                                                                <  5:25559/149‑1 (MQ=255)                                                                                        gTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGcc                                                                >  6:387421/1‑147 (MQ=255)                                                                                         ttGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAg                                                             <  8:308123/149‑1 (MQ=255)                                                                                             tGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTgg                                                         >  3:8879/1‑149 (MQ=255)                                                                                             tGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTgg                                                         >  1:86367/1‑149 (MQ=255)                                                                                                 aaCGGGTTGAATGAATGCGTGATCGCCGTACTGGAGCGCCCTCTGGTCGCCCATGAGCTGATCAACGTCGAGATTCGCTCGGAAGCCCGCGAAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGc                                                      <  1:106967/149‑1 (MQ=255)                                                                                                          aatgaaGGCGTGATCGCCGAACTGGAGCGCGCTCTGGGCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCgg                                            <  8:275024/149‑1 (MQ=255)                                                                                                           atgaaGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCgg                                            >  5:365966/1‑148 (MQ=255)                                                                                                                 gCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATc                                              >  5:341235/1‑140 (MQ=255)                                                                                                                 gCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATc                                              <  6:341235/140‑1 (MQ=255)                                                                                                                   gTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGc                                   >  6:110883/1‑149 (MQ=255)                                                                                                                   gTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACAGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGc                                   >  7:63359/1‑149 (MQ=255)                                                                                                                     gATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGCAAGccc                                 >  8:295259/1‑149 (MQ=255)                                                                                                                        cGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGa                              >  1:263654/1‑149 (MQ=255)                                                                                                                            gAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTa                          >  3:294722/1‑149 (MQ=255)                                                                                                                                  gAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCaa                    <  2:86367/149‑1 (MQ=255)                                                                                                                                    gcgcgcTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGa                  <  4:294722/149‑1 (MQ=255)                                                                                                                                      gcgcTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAAc                <  2:263654/149‑1 (MQ=255)                                                                                                                                        gcgaatctCGACCTTGATCAGCTCATGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACTGCAAGAAccc              >  8:344461/8‑149 (MQ=255)                                                                                                                                          tctGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCcac            >  4:32274/1‑149 (MQ=255)                                                                                                                                            tGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAg          >  7:266979/1‑149 (MQ=255)                                                                                                                                             ggTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAAc                <  7:311830/142‑1 (MQ=255)                                                                                                                                             ggTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAAc                >  8:311830/1‑142 (MQ=255)                                                                                                                                                cGACCTTGATCAGCTCATGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCa                                                     <  8:90330/102‑1 (MQ=255)                                                                                                                                                cGACCTTGATCAGCTCATGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCa                                                     >  7:90330/1‑102 (MQ=255)                                                                                                                                                 gACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCt                           >  2:126070/1‑127 (MQ=255)                                                                                                                                                 gACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACCCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCt                           <  1:126070/127‑1 (MQ=255)                                                                                                                                                  aCCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGaa     >  3:159837/1‑148 (MQ=255)                                                                                                                                                   ccATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAAc    >  7:213869/1‑148 (MQ=255)                                                                                                                                                   ccATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGAGCTACCGCAAGCACCCACCGCCGAAc    >  8:129384/1‑148 (MQ=255)                                                                                                                                                    cATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACaa  <  7:38073/149‑1 (MQ=255)                                                                                                                                                    cATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACaa  <  3:32274/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GACGGAAAGGGTCAGATTATGCCGCTCAATAACGAGCAGAAGAAGCAATACAAGTCCATTGGTCATGACCTGAAGCCGGTCCTGATCGTTGCTGGCAACGGTTTGAATGAAGGCGTGATCGCCGAACTGGAGCGCGCTCTGGTCGACCATGAGCTGATCAAGGTCGAGATTCGCTCGGAAGACCGCG‑AAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAA  >  NC_002947/5367461‑5367756

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