Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,458,940:1	+A	14.5%	intergenic (+17/‑45)	leuS → / → PP_4795	leucine‑‑tRNA ligase/lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,458,940	1	.	A	14.5%	219.3 / 37.2	69	intergenic (+17/‑45)	leuS/PP_4795	leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (29/30);  new base A (6/4);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAT  >  NC_002947/5458793‑5459088                                                                                                                                                     |                                                                                                                                                     ccGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGAAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGCGCAAGCTGGTCAACATCGACTCTAACTGATGAAAACGCCCACCCGC‑At                                                                                                                                                    <  2:412637/150‑1 (MQ=255)          gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCgg                                                                                                                                           <  2:396312/150‑1 (MQ=255)           ccGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACTCCCACCCGC‑ATCGATTGCggg                                                                                                                                          >  2:330710/1‑150 (MQ=255)             gcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTg                                                                                                                                        <  2:800509/150‑1 (MQ=255)              cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTgg                                                                                                                                       <  2:367032/150‑1 (MQ=255)               gaagaaGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTggg                                                                                                                                      >  2:607442/1‑150 (MQ=255)                 agaagTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCTGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCg                                                                                                                                    >  2:1408755/1‑150 (MQ=255)                    agTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGaa                                                                                                                                 <  1:764205/150‑1 (MQ=255)                    agTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGaa                                                                                                                                 <  1:291657/150‑1 (MQ=255)                      tCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACa                                                                                                                               <  1:773539/150‑1 (MQ=255)                       cGAAGCCGCGGCCCGCAGCAACGAGAATGTCATGCGCTTCATCGATGGCCTGACCATCCGCAAGGGCATCGTGGTAACGGGCAAGCTGGCCAACACCGTCGCCAACTGATGAAAACGCCCACCCCC‑ATCGATTGCTGGTGGGCGGaaacg                                                                                                                              <  1:1120617/150‑4 (MQ=255)                          aGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGcc                                                                                                                           >  2:1386096/1‑150 (MQ=255)                                          aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTg                                                                                                                                        <  2:549146/121‑1 (MQ=255)                                          aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTg                                                                                                                                        >  1:549146/1‑121 (MQ=255)                                          aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCaaaa                                                                                                                      <  2:634666/139‑1 (MQ=255)                                          aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCaaaa                                                                                                                      >  1:634666/1‑139 (MQ=255)                                          aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaacaa                                                                                                           <  1:607442/150‑1 (MQ=255)                                             gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGGGCAACAACAt                                                                                                        <  2:537027/150‑1 (MQ=255)                                             gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGGGCAACAACAt                                                                                                        >  1:537027/1‑150 (MQ=255)                                                aaTGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGCCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGAt                                                                                                     >  1:1355027/1‑150 (MQ=255)                                                  tGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                  <  2:141266/135‑1 (MQ=255)                                                  tGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                  >  1:141266/1‑135 (MQ=255)                                                   gTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaa                                                                                                  <  2:165160/150‑1 (MQ=255)                                                     ccTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCACGCTAGTGCACCTCGGCGGTGACTCCTGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAAc                                                                                                <  1:460871/150‑1 (MQ=255)                                                            ttCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGCGCAACAACATGATCAAACGCAAttt                                                                                         >  2:884358/1‑150 (MQ=255)                                                               aTCGATGGCCTGACCATCCGCAAGGTCATAGTGGTACCGGGCAAGCTGGTCAACATCGTCACCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTg                                                                                                                                        >  1:705678/1‑100 (MQ=255)                                                               aTCGATGGCCTGACCATCCGCAAGGTCATAGTGGTACCGGGCAAGCTGGTCAACATCGTCACCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTg                                                                                                                                        <  2:705678/100‑1 (MQ=255)                                                                  gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaa                                                                                                              >  2:1167075/1‑123 (MQ=255)                                                                  gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaa                                                                                                              <  1:1167075/123‑1 (MQ=255)                                                                        cTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCgg                                                                             >  2:682032/1‑150 (MQ=255)                                                                        cTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCg                                                                              >  1:425680/1‑149 (MQ=255)                                                                        cTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCg                                                                              <  2:425680/149‑1 (MQ=255)                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCggg                                                                                                                                          <  1:1382602/84‑1 (MQ=255)                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCggg                                                                                                                                          >  2:1382602/1‑84 (MQ=255)                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAgg                                                                                                                    <  1:1057163/106‑1 (MQ=255)                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAgg                                                                                                                    >  2:1057163/1‑106 (MQ=255)                                                                               tCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCg                                                                      >  1:76819/1‑150 (MQ=255)                                                                                  gCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTc                                                                   >  2:913651/1‑150 (MQ=255)                                                                                          aTCGTGGTGCCGTGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCaaaa                                                                                                                      <  1:602605/91‑1 (MQ=255)                                                                                          aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCaaaa                                                                                                                      >  2:602605/1‑91 (MQ=255)                                                                                          aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTg                                                                         <  2:8455/136‑1 (MQ=255)                                                                                          aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTg                                                                         >  1:8455/1‑136 (MQ=255)                                                                                          aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAg                                                           >  1:64557/1‑150 (MQ=255)                                                                                          aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAg                                                           >  1:172435/1‑150 (MQ=255)                                                                                                     gggCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGtt                                                <  2:76819/150‑1 (MQ=255)                                                                                                         aaGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCa                                            >  1:899941/1‑150 (MQ=255)                                                                                                          aGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAg                                           <  2:1355027/150‑1 (MQ=255)                                                                                                             tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                  <  2:943812/76‑1 (MQ=255)                                                                                                             tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                  >  1:943812/1‑76 (MQ=255)                                                                                                                tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcgc                                     <  2:172435/150‑1 (MQ=255)                                                                                                                  aaCATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTTGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCgg                                   >  1:784833/1‑150 (MQ=255)                                                                                                                              aaCTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGACCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCAc                       >  2:903444/1‑150 (MQ=255)                                                                                                                               aCTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGCCCGCTGCGGACGAATGGTGATCAACTACATGATCCAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCAcc                      <  1:1386096/150‑1 (MQ=255)                                                                                                                                 tgatgaAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCACGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGa                    <  1:330710/150‑1 (MQ=255)                                                                                                                                   atgaAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                     >  1:72663/1‑83 (MQ=255)                                                                                                                                   atgaAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                     <  2:72663/83‑1 (MQ=255)                                                                                                                                     gaAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTg                >  1:898223/1‑150 (MQ=255)                                                                                                                                       ttccGCCCACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCgg                                                                     >  2:876995/4‑95 (MQ=38)                                                                                                                                       ttccGCCCACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCgg                                                                     <  1:876995/92‑1 (MQ=38)                                                                                                                                        aaCGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGc             <  1:682032/150‑1 (MQ=255)                                                                                                                                        aaCGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGc             >  1:116005/1‑150 (MQ=255)                                                                                                                                          cGCCCACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGACGGTAATCCTAACCGCCAGAGTTTTGCTGCTGCCTAGCACTGACTCCACCGGGCTGAGCGt           >  2:514333/1‑150 (MQ=255)                                                                                                                                           gCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTg          >  1:1340781/1‑150 (MQ=255)                                                                                                                                             ccACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa        >  2:1216614/1‑150 (MQ=255)                                                                                                                                             ccACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa        >  1:315340/1‑150 (MQ=255)                                                                                                                                             ccACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa        <  2:898223/150‑1 (MQ=255)                                                                                                                                              cACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCAcc                      <  2:1083190/135‑1 (MQ=255)                                                                                                                                              cACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCAcc                      >  1:1083190/1‑135 (MQ=255)                                                                                                                                                   gC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAt  <  2:315340/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAT  >  NC_002947/5458793‑5459088

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