Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,408,458	G→T	6.3%	H94N (CAC→AAC)	PP_4752 ←	Xaa‑Pro aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,408,458	0	G	T	6.3%	200.5 / 3.3	63	H94N (CAC→AAC)	PP_4752	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base T (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GGGGGCCGGCGAGGAAATACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCG  >  NC_002947/5408311‑5408596                                                                                                                                                      |                                                                                                                                           gggggCCGGCGAGGAAATACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGc                                                                                                                                           <  6:99007/149‑1 (MQ=255) gggggCCGGCGAGGAAATACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGc                                                                                                                                           <  5:241726/149‑1 (MQ=255)              gAAATACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAAc                                                                                                                              <  4:269179/149‑1 (MQ=255)              gAAATACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAAc                                                                                                                              >  8:261448/1‑149 (MQ=255)                aaTACAGCCAGGGAAACGCCGGGCCAAATTGATCGATGGTTTCAACGTGCACGCTGATGTGGTGGCAGCTGGTGAACTCGAACA‑GAACCACCGGGCCGTCG‑GTGGCGCCCCACACATAGCGGCTGGGCGAGTGCAACACCCACAAAcc                                                                                                                             >  8:122832/1‑148 (MQ=255)                 aTACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCaa                                                                                                                           <  3:299473/149‑1 (MQ=255)                       ccAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGTCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTAt                                                                                                                     <  7:261448/149‑1 (MQ=255)                                 gCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCGAGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCa                                                                                                                                       <  1:53684/122‑1 (MQ=255)                                 gCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCa                                                                                                                                       >  2:53684/1‑121 (MQ=255)                                 gCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGt                                                                                                           >  6:305306/1‑149 (MQ=255)                                   cGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGt                                                                                                           <  2:262018/147‑1 (MQ=255)                                    gggCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGc                                                                                                        <  7:122832/149‑1 (MQ=255)                                           tttCATCGATCGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGa                                                                                                 <  6:267038/149‑1 (MQ=255)                                             tCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGaa                                                                                                <  4:262097/148‑1 (MQ=255)                                              cATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAata                                                                                              >  2:232008/1‑149 (MQ=255)                                                   aTGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgat                                                                                         >  5:83222/1‑149 (MQ=255)                                                   aTGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgat                                                                                         <  1:232008/149‑1 (MQ=255)                                                    tgggtgcAACGTGCGCGCGGTTGTTCTGGCTGCTGGGGAACGCGCACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACCACCCACGGGTATTGGTGTCGGTGGCATAGCGAATATTgatc                                                                                        <  7:265875/143‑1 (MQ=255)                                                    tGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgatc                                                                                        <  1:142515/149‑1 (MQ=255)                                                    tGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgatc                                                                                        <  4:236929/149‑1 (MQ=255)                                                          cAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgatcggatc                                                                                   <  5:305306/148‑1 (MQ=255)                                                            aCGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGc                                                                                >  6:252600/1‑149 (MQ=255)                                                            aCGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGc                                                                                <  7:75434/149‑1 (MQ=255)                                                                   cGCTGTTGTGCTGGCTGCTGGTGAACTCGAACATGCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgat                                                                                         >  1:41113/1‑134 (MQ=255)                                                                   cGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAat                                                                                               >  6:291476/1‑127 (MQ=255)                                                                   cGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAat                                                                                               <  5:291476/127‑1 (MQ=255)                                                                   cGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTgat                                                                                         <  2:41113/133‑1 (MQ=255)                                                                               ggctgCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTCGTGTCGGTGGCATAGCGAAAATTGATCGGATCGGCCAACAGAAAGCCGGCGTa                                                              >  3:114080/1‑148 (MQ=255)                                                                                   gctgGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCt                                                         >  5:146298/1‑149 (MQ=255)                                                                                     tgGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTgg                                                       >  5:215646/1‑149 (MQ=255)                                                                                       gTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACTAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGgc                                                     <  5:52420/149‑1 (MQ=255)                                                                                       gTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAAc                                                                            >  3:200415/1‑126 (MQ=255)                                                                                       gTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAAc                                                                            <  4:200415/126‑1 (MQ=255)                                                                                        tGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGgcg                                                    >  1:258070/1‑149 (MQ=255)                                                                                         gAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGgcgc                                                   >  4:33774/1‑149 (MQ=255)                                                                                             tCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAg                                               >  5:243244/1‑149 (MQ=255)                                                                                             tCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAg                                               >  5:207555/1‑149 (MQ=255)                                                                                             tCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAg                                               >  7:25987/1‑149 (MQ=255)                                                                                             tCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAg                                               >  6:168557/1‑149 (MQ=255)                                                                                                 aCA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGc                                           >  7:8386/1‑149 (MQ=255)                                                                                                       accaccGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTaa                                                             <  7:49590/126‑1 (MQ=255)                                                                                                       accaccGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTaa                                                             >  8:49590/1‑126 (MQ=255)                                                                                                       accaccGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTcgc                                       <  6:83222/148‑1 (MQ=255)                                                                                                         caccGCGCCTGCG‑GGGGAGACGAACACATAGCTGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGaa                                    <  6:215646/149‑1 (MQ=255)                                                                                                          accGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAAt                                   >  7:285071/1‑149 (MQ=255)                                                                                                              ggCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATgcgc                               >  5:220763/1‑149 (MQ=255)                                                                                                                  gTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCa                           <  5:168557/149‑1 (MQ=255)                                                                                                                   tCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAg                          >  4:167156/1‑149 (MQ=255)                                                                                                                          gCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGt                     >  2:278427/1‑148 (MQ=255)                                                                                                                           cgacgaACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTaa                   >  6:260220/1‑149 (MQ=255)                                                                                                                             acgaACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATg                 <  8:96440/149‑1 (MQ=255)                                                                                                                                aaCACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGgcg              <  8:25987/149‑1 (MQ=255)                                                                                                                                     aTAGCGGCTGGGCGAGTTCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAgg         >  1:303105/1‑149 (MQ=255)                                                                                                                                     aTAGCGGCTGGGCGAGTTCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAgg         >  1:298339/1‑149 (MQ=255)                                                                                                                                     aTAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAgg         >  5:57890/1‑149 (MQ=255)                                                                                                                                        gCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGcgt      >  8:257154/1‑147 (MQ=255)                                                                                                                                         cGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCttt     <  5:252600/149‑1 (MQ=255)                                                                                                                                         cGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCttt     >  3:144677/1‑149 (MQ=255)                                                                                                                                          ggCTGGGCGAGTTCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTg    <  2:303105/149‑1 (MQ=255)                                                                                                                                          ggCTGGGCGAGTTCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTg    <  2:298339/149‑1 (MQ=255)                                                                                                                                          ggCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTg    <  6:220763/149‑1 (MQ=255)                                                                                                                                           gCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGc   >  8:226581/1‑149 (MQ=255)                                                                                                                                            cTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGcg  <  5:260220/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                           GGGGGCCGGCGAGGAAATACAGCCAGGGAATCGCCGGGCGAATTTCATCGATGGTTTCAACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACA‑GCACCACCGGGCCGTCG‑GTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCCAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCG  >  NC_002947/5408311‑5408596

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