Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	464,325	G→C	11.8%	H259D (CAC→GAC) ‡	davB ←	lysine 2‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	464,325	0	G	C	11.8%	191.1 / 13.8	66	H259D (CAC→GAC) ‡	davB	lysine 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (3/5);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGC  >  NC_002947/464177‑464467                                                                                                                                                     |                                                                                                                                               gcggcgCGGACAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGc                                                                                                                                               <  2:1283159/150‑1 (MQ=255)    gcgcGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGCACGTCGAGCCAGATGCCTTGTGGCACTTGTCCCACGCCCCCGACAACCAGGTGCggg                                                                                                                                            <  2:1817378/150‑1 (MQ=255)      gcgGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtggtg                                                                                                                                          <  1:457435/150‑1 (MQ=255)      gcgGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtggtg                                                                                                                                          <  2:1713627/150‑1 (MQ=255)      gcgGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCAGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCtggtg                                                                                                                                          <  2:1582308/150‑1 (MQ=255)         ggCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtc                                                                                                                                       <  2:1121298/150‑1 (MQ=255)          gCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcg                                                                                                                                      <  2:364463/150‑1 (MQ=255)          gCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcg                                                                                                                                      >  2:436339/1‑150 (MQ=255)           cAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcgg                                                                                                                                     <  2:481071/150‑2 (MQ=255)            aaTGCGCTTGACCCCGGTACGCGGTGAGCCTCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCGATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGgtcgtc                                                                                                                                    >  1:349064/1‑150 (MQ=255)                cgcTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAg                                                                                                                                <  2:128152/150‑1 (MQ=255)                 gcTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGt                                                                                                                               >  2:1233058/1‑150 (MQ=255)                 gcTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgggggcgt                                                                                                                               <  1:439033/150‑8 (MQ=255)                        cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCCACGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTCCGGACACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtcat                                                                                                                        >  2:730592/1‑150 (MQ=255)                        cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtcat                                                                                                                        >  1:752546/1‑150 (MQ=255)                        cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                     >  2:1762161/1‑136 (MQ=255)                        cccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                     <  1:1762161/137‑2 (MQ=255)                         ccGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGtcatc                                                                                                                       <  1:878684/150‑1 (MQ=255)                           ggTACGCGGTGCGCCGCCATGCAGCGAGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCATTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATcac                                                                                                                     >  2:962139/1‑150 (MQ=255)                                gcgGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcgc                                                                                                                >  2:616545/1‑150 (MQ=255)                                      cgccgcCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGa                                                                                                          <  1:1890097/150‑1 (MQ=255)                                        ccgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGAtt                                                                                                        >  1:1273251/1‑150 (MQ=255)                                        ccgccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTAGCAGTTGGTCATCACCACGCGGAAGAtt                                                                                                        <  1:956431/150‑1 (MQ=255)                                           ccATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTcc                                                                                                     >  2:1700812/1‑150 (MQ=255)                                            cATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCa                                                                                                    <  1:1700812/150‑1 (MQ=255)                                                 aGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                     >  2:1426952/1‑111 (MQ=255)                                                 aGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                     <  1:1426952/112‑2 (MQ=255)                                                    gTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAg                                                                                            >  2:1101805/1‑150 (MQ=255)                                                    gTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAg                                                                                            <  1:1733517/150‑1 (MQ=255)                                                        tCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTgg                                                                                        <  1:436339/150‑1 (MQ=255)                                                        tCAGGCTGGTGCCCTCTGGCCAATACACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTgg                                                                                        >  1:1424982/1‑150 (MQ=255)                                                            gCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGaa                                                                                    <  1:616545/150‑1 (MQ=255)                                                            gCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGaa                                                                                    >  2:1708390/1‑150 (MQ=255)                                                                gTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTcc                                                                                <  1:1708390/150‑1 (MQ=255)                                                                gTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCAACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTcc                                                                                <  2:1300757/150‑1 (MQ=255)                                                                     ctctGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCATGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGGTGGGGAAGTCCGAGTc                                                                           >  2:768313/1‑150 (MQ=255)                                                                        tGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCa                                                                        >  1:1120446/1‑150 (MQ=255)                                                                                  aCGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTg                                                              >  1:296201/1‑150 (MQ=255)                                                                                  aCGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTg                                                              >  2:1470557/1‑150 (MQ=255)                                                                                  aCGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTg                                                              >  2:1450821/1‑150 (MQ=255)                                                                                         gTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCATGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGCGAAGTCCGAGTCCCAACCGCCGGTGCCGAAAc                                                       <  2:1070676/150‑1 (MQ=255)                                                                                               gTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCt                                                 >  1:606635/1‑150 (MQ=255)                                                                                                      gcgcCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCCGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGaa                                          <  1:1562282/150‑1 (MQ=255)                                                                                                       cgcCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAac                                         <  1:342298/150‑1 (MQ=255)                                                                                                               gCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTcgcg                                 <  2:1120446/150‑1 (MQ=255)                                                                                                                 cTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGt                               <  1:962139/150‑1 (MQ=255)                                                                                                                  ttGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcgc                                                                                                                >  1:23367/1‑68 (MQ=255)                                                                                                                  ttGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCAcgc                                                                                                                <  2:23367/68‑1 (MQ=255)                                                                                                                  ttGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTg                              >  1:250611/1‑150 (MQ=255)                                                                                                                   tgtgTCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGc                             <  1:876063/150‑1 (MQ=255)                                                                                                                   tgtgGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGc                             <  1:588087/150‑1 (MQ=255)                                                                                                                          cTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCgc                                                                                             >  1:223984/1‑78 (MQ=255)                                                                                                                          cTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGa                                                                                             <  2:223984/79‑1 (MQ=255)                                                                                                                               tCCACGCCCCCGTCCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGTGGAAGTCGGAGTCCCAACCGCCGGTGCCGAAACCGCCCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGc                 <  2:1056788/150‑1 (MQ=255)                                                                                                                               tCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGc                 >  2:1670879/1‑150 (MQ=255)                                                                                                                                  aCGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCt                                                 >  2:982891/1‑115 (MQ=255)                                                                                                                                  aCGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCt                                                 <  1:982891/115‑1 (MQ=255)                                                                                                                                  aCGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTg              >  1:1246790/1‑150 (MQ=255)                                                                                                                                   cGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCGGTTACCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTgg             <  1:1482748/150‑1 (MQ=255)                                                                                                                                          gACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTg              <  2:563918/142‑1 (MQ=255)                                                                                                                                          gACCACCAGCACCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTg              >  1:563918/1‑142 (MQ=255)                                                                                                                                            ccaccaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgc    <  2:296201/150‑1 (MQ=255)                                                                                                                                              accaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgcgc  >  2:1443267/1‑150 (MQ=255)                                                                                                                                              accaGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAgcgc  >  1:151775/1‑150 (MQ=255)                                                                                                                                                 aGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCaa                                                                       >  2:545560/1‑78 (MQ=255)                                                                                                                                                 aGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCaa                                                                       <  1:545560/78‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGC  >  NC_002947/464177‑464467

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