Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	464,323	G→C	11.9%	H259Q (CAC→CAG) ‡	davB ←	lysine 2‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	464,323	0	G	C	11.9%	190.1 / 15.8	67	H259Q (CAC→CAG) ‡	davB	lysine 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base C (3/5);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CGGCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGC  >  NC_002947/464175‑464467                                                                                                                                                     |                                                                                                                                                 cggcggCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGt                                                                                                                                                 <  1:668762/150‑1 (MQ=255)   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                                                                                                                                     <  1:1426952/112‑2 (MQ=255)                                                   aGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGCACCTGGTGgtcgg                                                                                                                                     >  2:1426952/1‑111 (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                                                                                    >  2:1708390/1‑150 (MQ=255)                                                              gCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGaa                                                                                    <  1:616545/150‑1 (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:1450821/1‑150 (MQ=255)                                                                                    aCGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTg                                                              >  2:1470557/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                                                 <  1:982891/115‑1 (MQ=255)                                                                                                                                    aCGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCt                                                 >  2:982891/1‑115 (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)                                                                                                                                                     |                                                                                                                                                 CGGCGGCGCGGGCAATGCGCTTGACCCCGGTACGCGGTGCGCCGCCATGCAGCGTGCTCAGGCTGGTGCCCTCTGGCCAATGCACGCAGCGTTCCGGTACGTCGCGCCAGATGCCTTGTGGCACTTGTTCCACGCCCCCGACCACCAGGTGCTGGTGGTCGTCGCAGTTGGTCATCACCACGCGGAAGATTTCCAGCATCGAGTTGGGGAAGTCCGAGTCCCAACCGCCGGTGCCGAAACCGACCTGGCCGAACACTTCGCGGTGCTGGAAGCTCAGCTTGGCAAAAGAGCGC  >  NC_002947/464175‑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: ‑