Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,458,892:1	+A	8.7%	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,892	1	.	A	8.7%	258.2 / 18.0	69	intergenic (+17/‑45)	leuS/PP_4795	leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (25/38);  new base A (3/3);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

CCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAA  >  NC_002947/5458745‑5459039                                                                                                                                                      |                                                                                                                                                    ccGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑a                                                                                                                                                    >  4:315138/1‑149 (MQ=255)  cGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑At                                                                                                                                                   <  4:3843/149‑1 (MQ=255)   gCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATc                                                                                                                                                  <  2:26470/149‑1 (MQ=255)    cTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATc                                                                                                                                                  >  3:134010/1‑148 (MQ=255)       ccAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGa                                                                                                                                                <  1:208395/147‑1 (MQ=255)        cAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTg                                                                                                                                             >  4:15258/1‑149 (MQ=255)          gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCg                                                                                                                                           <  5:277209/149‑1 (MQ=255)          gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCg                                                                                                                                           >  6:46532/1‑149 (MQ=255)             gcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGt                                                                                                                                        <  4:134010/149‑1 (MQ=255)             gcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGt                                                                                                                                        >  3:206039/1‑149 (MQ=255)               gaagaaGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTgg                                                                                                                                      <  1:349006/149‑1 (MQ=255)               gaagaaCTCGAAGCCGCTCCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATCGCCTGACCAACCGCAACGCCATCCTGGTGCCGGGCAACCTCCTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCCATTGCGCGtcg                                                                                                                                      >  4:158309/1‑147 (MQ=255)                 agaagTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTTACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTGGGc                                                                                                                                    >  6:263061/1‑149 (MQ=255)                  gaagTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCg                                                                                                                                   >  4:350318/1‑149 (MQ=255)                    agTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTCATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCgg                                                                                                                                  <  1:356231/148‑1 (MQ=255)                        gAAGCCGCTGCCCGCAGCAACGAGAATGGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGaa                                                                                                                                >  1:209837/1‑147 (MQ=255)                        gAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTTACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTGGGCgttttc                                                                                                                              <  5:263061/148‑6 (MQ=255)                              gCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCATGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCa                                                                                                                        >  7:149562/1‑148 (MQ=255)                                  cccGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCaaaa                                                                                                                     <  8:201958/147‑1 (MQ=255)                                   ccGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAgg                                                                                                                   >  6:277904/1‑148 (MQ=255)                                    cGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                 <  4:206039/149‑1 (MQ=255)                                    cGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                 <  5:330403/149‑1 (MQ=255)                                    cGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                 <  3:83994/149‑1 (MQ=255)                                    cGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                 <  7:95533/149‑1 (MQ=255)                                         caACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaac                                                                                                            <  2:103766/149‑1 (MQ=255)                                          aaCGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaac                                                                                                            >  7:259321/1‑148 (MQ=255)                                           aCGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGCCATAGTGGTGCCGGGCAAGCTGGTCAACATCGTCGCCAACTGATTAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaacaa                                                                                                          <  3:158309/149‑1 (MQ=255)                                                    gTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATc                                                                                                   >  3:315794/1‑148 (MQ=255)                                                       cTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaa                                                                                                 >  5:183111/1‑147 (MQ=255)                                                       cTGCGCTTCATCGATGGCATGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                                <  8:149562/148‑1 (MQ=255)                                                                aTCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTg                                                                                       <  5:277904/148‑1 (MQ=255)                                                                 tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                       <  5:198399/131‑1 (MQ=255)                                                                 tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                       >  6:198399/1‑131 (MQ=255)                                                                 tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                     >  6:148067/1‑149 (MQ=255)                                                                 tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                     <  4:315794/149‑1 (MQ=255)                                                                 tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCt                                                                                     <  3:303035/149‑1 (MQ=255)                                                                   gATTGCATGACCATACGAAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCATCTGATGTAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                    <  6:183111/148‑1 (MQ=255)                                                                   gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                    <  7:290170/148‑1 (MQ=255)                                                                   gATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                    <  8:34780/148‑1 (MQ=255)                                                                       gCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGGTCAAACGCAATTTGCTGGTaa                                                                                <  2:209837/148‑1 (MQ=255)                                                                             ccATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGcc                                                                          <  3:315138/148‑1 (MQ=255)                                                                               aTCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaac                                                                                                            >  4:63104/1‑112 (MQ=255)                                                                               aTCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGcaac                                                                                                            <  3:63104/112‑1 (MQ=255)                                                                               aTCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                       <  5:254394/149‑1 (MQ=255)                                                                                tCCGCAAGGACATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAAAGCAAGTTGCTGGTAATCGCCCTGGc                                                                      >  2:224823/1‑149 (MQ=255)                                                                                   gCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGt                                                                   >  5:34454/1‑149 (MQ=255)                                                                                   gCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGt                                                                   >  6:69760/1‑149 (MQ=255)                                                                                     aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCa                                                                 <  3:350318/149‑1 (MQ=255)                                                                                      aGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCAt                                                                <  8:259321/149‑1 (MQ=255)                                                                                       ggTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATg                                                               <  5:148067/149‑1 (MQ=255)                                                                                          cATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCt                                                             >  2:354836/1‑148 (MQ=255)                                                                                           aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTc                                                            <  5:69760/148‑1 (MQ=255)                                                                                            tCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAg                                                          >  8:239415/1‑149 (MQ=255)                                                                                              gTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                         <  5:46532/148‑1 (MQ=255)                                                                                                     cGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCgg                                                 <  6:34454/149‑1 (MQ=255)                                                                                                          aaGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTcc                                            <  4:177906/149‑1 (MQ=255)                                                                                                                 tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGTGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcg                                     >  6:359991/1‑149 (MQ=255)                                                                                                                     cATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGc                                 <  3:15258/149‑1 (MQ=255)                                                                                                                                    atgaAAACGCCCACCCGC‑ATCGATTGTGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAg                  <  5:359991/149‑1 (MQ=255)                                                                                                                                       gttccGCCCACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTcc                                            >  2:52572/5‑120 (MQ=38)                                                                                                                                       gttccGCCCACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTcc                                            <  1:52572/116‑1 (MQ=38)                                                                                                                                         aaCGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAg             >  8:206974/1‑149 (MQ=255)                                                                                                                                         aaCGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAg             <  1:224823/149‑1 (MQ=255)                                                                                                                                            gCCCACCCGCAATCGA‑TGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGt          >  8:180136/1‑149 (MQ=255)                                                                                                                                              ccACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGa                   >  8:189756/1‑138 (MQ=255)                                                                                                                                              ccACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGa                   <  7:189756/138‑1 (MQ=255)                                                                                                                                               cACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa       <  1:354836/149‑1 (MQ=255)                                                                                                                                                   cGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGa    >  2:289397/1‑148 (MQ=255)                                                                                                                                                    gC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGaaa  <  7:239415/149‑1 (MQ=255)                                                                                                                                                    gC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGaaa  >  3:308433/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAA  >  NC_002947/5458745‑5459039

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