Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,632,097	A→G	8.3%	E284G (GAG→GGG)	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,632,097	0	A	G	8.3%	161.0 / 4.6	60	E284G (GAG→GGG)	lepA	elongation factor 4
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

ACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCG  >  NC_002947/1631950‑1632239                                                                                                                                                      |                                                                                                                                                    aCTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCAAACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGAATCTAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAg                                                                                                                                                    <  7:269770/149‑1 (MQ=255)     ggTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTc                                                                                                                                                 <  7:71701/148‑1 (MQ=255)       tCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAg                                                                                                                                              <  8:183066/149‑1 (MQ=255)        cGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAgg                                                                                                                                             >  3:307296/1‑149 (MQ=255)         gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCTCAGGGGTCGAg                                                                                                                                              <  8:84667/147‑1 (MQ=255)         gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCTCAGGGGTCGAg                                                                                                                                              >  7:84667/1‑147 (MQ=255)             gCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCTCAGGGGTCGAggacagg                                                                                                                                        <  6:116582/149‑6 (MQ=255)              cGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGc                                                                                                                                       <  1:145080/149‑1 (MQ=255)               gTCGGTGTCTTCACCCCGAAGCACACCCAGACCTGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGCCCCCTGAGGTCGAGgtgt                                                                                                                                          >  1:225304/1‑146 (MQ=255)               gTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGgtgt                                                                                                                                          <  2:225304/145‑1 (MQ=255)                   gtgtCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAACATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGc                                                                                                                                  <  7:193299/149‑1 (MQ=255)                        ttCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTc                                                                                                                               >  2:269066/1‑147 (MQ=255)                           aCCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                     >  3:170357/1‑138 (MQ=255)                           aCCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                     <  4:170357/138‑1 (MQ=255)                           aCCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaagaa                                                                                                                          <  4:343787/149‑1 (MQ=255)                             cccGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                            >  2:136671/1‑145 (MQ=255)                             cccGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                        <  1:269066/149‑1 (MQ=255)                             cccGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGGGACACCCTGACA‑CTGTCCACGCCCCGGGAGGTCGAGGTGTTGCCAGGCTTCAagatga                                                                                                                        >  4:33889/1‑149 (MQ=255)                              ccGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGa                                                                                                                                                     >  5:270053/1‑119 (MQ=255)                              ccGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGa                                                                                                                                                     <  6:270053/119‑1 (MQ=255)                               cGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAtt                                                                                                                      <  5:241602/149‑1 (MQ=255)                                         ccAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGc                                                                                                                                       >  7:86495/1‑122 (MQ=255)                                         ccAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGc                                                                                                                                       <  8:86495/122‑1 (MQ=255)                                                 gCTGATCTGATAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGTACATCCACTGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAATATTCAGCCGCTTGTCTACGc                                                                                                     <  2:31259/149‑1 (MQ=255)                                                     aTCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCg                                                                                                   <  4:307296/147‑1 (MQ=255)                                                        tGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTg                                                                                              >  1:36042/1‑149 (MQ=255)                                                         gAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                             >  2:175881/1‑149 (MQ=255)                                                          aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCt                                                                                               <  8:40979/146‑1 (MQ=255)                                                          aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                             <  5:349070/148‑1 (MQ=255)                                                          aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTGCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGtt                                                                                            >  3:154474/1‑149 (MQ=255)                                                            aGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTcc                                                                                          >  6:124676/1‑149 (MQ=255)                                                            aGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTcc                                                                                          <  1:175881/149‑1 (MQ=255)                                                              ccGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑g                                                                                       <  2:351958/149‑1 (MQ=255)                                                                ggCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑Gt                                                                                   <  2:36042/149‑1 (MQ=255)                                                                ggCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCACTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑Gt                                                                                   >  5:155797/1‑149 (MQ=255)                                                                   gAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAg                                                                                >  2:155832/1‑149 (MQ=255)                                                                     aGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACTGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTc                                                                                                          <  7:110636/124‑1 (MQ=255)                                                                     aGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTc                                                                                                          >  8:110636/1‑124 (MQ=255)                                                                      gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTc                                                                             >  8:319568/1‑149 (MQ=255)                                                                         ggCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCcga                                                                          <  5:124676/149‑1 (MQ=255)                                                                                  atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGTAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCg                                                                <  1:136671/148‑1 (MQ=255)                                                                                     gCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAgg                                                             <  1:233514/148‑1 (MQ=255)                                                                                      ccAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCGCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCcgac                                                                         >  1:79469/1‑138 (MQ=255)                                                                                      ccAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCcgac                                                                         <  2:79469/137‑1 (MQ=255)                                                                                                      ccACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTg                                            >  7:31572/1‑148 (MQ=255)                                                                                                               cccGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                  >  7:71661/1‑149 (MQ=255)                                                                                                                ccGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCTGAAGCTGAc                                 >  1:334600/1‑149 (MQ=255)                                                                                                                ccGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAc                                 >  1:125152/1‑149 (MQ=255)                                                                                                                   gTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCt                              <  8:71661/149‑1 (MQ=255)                                                                                                                          acacCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCTGAAGCTGACCCTGAACGa                        <  2:334600/148‑1 (MQ=255)                                                                                                                              cctgacc‑ctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcct                   >  7:330613/1‑149 (MQ=255)                                                                                                                               ctgacc‑ctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                  >  6:74377/1‑149 (MQ=255)                                                                                                                                    c‑ctgTCCTCGACCCGTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGaa                           <  3:351492/135‑1 (MQ=255)                                                                                                                                    c‑ctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGaa                           >  4:351492/1‑135 (MQ=255)                                                                                                                                    c‑ctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTg              >  5:212430/1‑148 (MQ=255)                                                                                                                                      ctgTCCTCGACCCCTGAGGTCGAGGGGTTGCCAGGCTGCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa            <  8:189937/149‑1 (MQ=255)                                                                                                                                       tgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAg           >  6:124577/1‑149 (MQ=255)                                                                                                                                         accTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTa         >  2:350274/2‑149 (MQ=255)                                                                                                                                           cTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCAGAAGGAAGCTCCGACAGAGCTACTAGGATTTACTCGACGATCAGCAGAAGATGAACCAGAACAACTCCTCgaag              >  1:111610/1‑144 (MQ=255)                                                                                                                                             cGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAccc                               <  8:286529/123‑1 (MQ=255)                                                                                                                                               aCCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGcc   <  1:155832/149‑1 (MQ=255)                                                                                                                                                ccccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCg  <  2:125152/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACC‑GCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACC‑CTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCC‑G‑‑GTCAGCTCCGAC‑GA‑CTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCG  >  NC_002947/1631950‑1632239

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