Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,458,898	Δ1 bp	9.2%	intergenic (+23/‑39)	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,898	0	T	.	9.2%	196.3 / 14.5	65	intergenic (+23/‑39)	leuS/PP_4795	leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (24/35);  new base . (3/3);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAA  >  NC_002947/5458752‑5459039                                                                                                                                                     |                                                                                                                                              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                                                                                                                 <  7:95533/149‑1 (MQ=255)                             cGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                 <  3:83994/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                                                                                                       >  6:198399/1‑131 (MQ=255)                                                          tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                       <  5:198399/131‑1 (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)                                                                                                                                                     |                                                                                                                                              CAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAAT‑GTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAA  >  NC_002947/5458752‑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: ‑