Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,909,949	C→G	8.0%	intergenic (+342/+137)	PP_4318 → / ← PP_4319	transposase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,909,949	0	C	G	8.0%	174.1 / 3.4	61	intergenic (+342/+137)	PP_4318/PP_4319	transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (3/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACTGTGC  >  NC_002947/4909802‑4910094                                                                                                                                                       |                                                                                                                                                  tGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATc                                                                                                                                                   <  2:22562/148‑1 (MQ=32)     aaaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCAT‑CGCTAAGGTAATGGTCCCGGGA‑‑GTAAGCGCTGCCGGCCCCAACGAAGGGGCCTCCGCGCATGACGCGCATCCaggg                                                                                                                                              >  6:11841/1‑146 (MQ=11)       aCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                                            <  3:380532/149‑1 (MQ=32)            aCAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                       <  2:335579/149‑2 (MQ=17)                aaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTcc                                                                                                                                   <  4:418954/149‑1 (MQ=32)                   gTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTCAGCTACCCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTACCGCGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGATTCAAGGGAATCCATGCCCGTCACCTCCCt                                                                                                                                 >  5:442823/1‑148 (MQ=11)                       ccAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGc                                                                                                                             >  4:165491/1‑148 (MQ=34)                         aaaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAg                                                                                                                          <  3:165491/149‑1 (MQ=34)                          aaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCCTCACCACCCTTCGCCAGc                                                                                                                         >  1:44613/1‑149 (MQ=255)                                    gCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCg                                                                                                               <  6:261106/149‑1 (MQ=37)                                      tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGAAGGTACGCCCTGCAGGCCCCATCCATGGGGCCTCCGCGCTTGACGGGCAGCCATGCCCGGCCCCTCCCTTCGCAGGCCCTGCGGTCg                                                                                                               >  1:348233/1‑149 (MQ=11)                                      tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATACATGCCCGTCACCTCCCATCGCAAGCCCTGCGCACGGc                                                                                                             >  6:252341/1‑149 (MQ=255)                                        cTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcc                                                                                                            <  2:260911/148‑1 (MQ=35)                                         tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCCTGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCtt                                                                                                          >  5:182949/1‑149 (MQ=25)                                          atcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTACCGGGA‑‑GTACGCGCTGCAGGCCCCATCCCTGGGGCCTCAGCGCTTGACGGGCATGCATGGAGGTCTCCTCCCTTCTCAAGCCCTGCGCTCGGCCTTc                                                                                                         >  7:401820/1‑149 (MQ=11)                                          atcatcCGGCCCCTGCGCTACACTCCGGGGTTCCCTACGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCtt                                                                                                          >  1:72949/1‑149 (MQ=255)                                           tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCtt                                                                                                          <  6:90701/147‑1 (MQ=35)                                           tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTAAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCt                                                                                                        <  1:21228/149‑1 (MQ=37)                                                cGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAg                                                                                                    <  4:335501/148‑1 (MQ=37)                                                   ccccTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGc                                                                                                <  8:332381/149‑1 (MQ=37)                                                    cccTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGc                                                                                                >  5:202901/1‑148 (MQ=37)                                                     ccTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGATCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCAAGGGGCCTAATCGCTTGACGGGCATACAATCCCGTCACCTCCCATCGCAAGCCCAGCGCTCGGCCTTCTGAAGTCCCa                                                                                               >  6:239014/1‑148 (MQ=11)                                                      cTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAAt                                                                                             >  2:98902/1‑149 (MQ=37)                                                          gcTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                                                         <  3:172613/149‑1 (MQ=37)                                                              cGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTc                                                                                      <  4:150116/148‑1 (MQ=37)                                                               gCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                                                    >  7:407787/1‑149 (MQ=37)                                                                cTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGcc                                                                                   <  5:252341/149‑1 (MQ=37)                                                                     gggTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAAc                                                                              >  8:135349/1‑149 (MQ=37)                                                                     gggTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAAc                                                                              <  2:101455/149‑1 (MQ=37)                                                                        ttCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAAt                                                                                             >  8:160745/1‑131 (MQ=37)                                                                                gCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCgcgcgc                                                                    <  3:87702/149‑1 (MQ=40)                                                                                 cTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCt                                                                   >  6:71867/1‑149 (MQ=255)                                                                                    cGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCGATGGGGCCTCAGCGCTTGACGGGCGTGCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTa                                                                 >  7:133658/1‑148 (MQ=255)                                                                                            cTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCCCCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGAGTACAAGCaa                                                          >  5:234410/1‑147 (MQ=255)                                                                                            cTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCaaaa                                                        <  5:125997/149‑1 (MQ=40)                                                                                               ccGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaag                                                       >  8:99393/1‑147 (MQ=255)                                                                                               ccGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaag                                                       <  5:71867/147‑1 (MQ=255)                                                                                                 gggA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcaga                                                   <  5:239014/149‑1 (MQ=255)                                                                                                 gggA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcaga                                                   <  6:264915/149‑1 (MQ=255)                                                                                                   gA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaag                                                 <  7:52268/149‑1 (MQ=255)                                                                                                          cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcaga                                             <  3:367657/148‑1 (MQ=255)                                                                                                           gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                                           <  2:72949/149‑1 (MQ=255)                                                                                                           gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                                           >  5:383638/1‑149 (MQ=255)                                                                                                           gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                                           <  6:383638/149‑1 (MQ=255)                                                                                                             gcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                                           >  6:213114/1‑147 (MQ=255)                                                                                                             cgctgaGGCCCCATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTACTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                                           >  5:219281/6‑147 (MQ=255)                                                                                                              cTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcag                                        <  8:407787/149‑1 (MQ=255)                                                                                                                         atccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcagaa                             <  6:182949/149‑1 (MQ=255)                                                                                                                          tccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcagaag                            <  6:442823/149‑1 (MQ=255)                                                                                                                                     tCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGACGCAGAAGCGGCAGCAGCAGAGAcgg                   >  4:347948/1‑145 (MQ=255)                                                                                                                                           cTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAgtg           >  8:109612/1‑149 (MQ=255)                                                                                                                                            gtGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAACTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAACCAGAAGCAGAAGCAGAAGCAGCAGAAGCACAAACAGAAACAgtc           >  8:907/2‑147 (MQ=18)                                                                                                                                              gACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAgtgt          <  2:348233/147‑1 (MQ=255)                                                                                                                                                 gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtg     >  1:57005/1‑149 (MQ=255)                                                                                                                                                 gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtg     >  1:243282/1‑149 (MQ=255)                                                                                                                                                 gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtg     >  4:125140/1‑149 (MQ=255)                                                                                                                                                 gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtg     >  6:47123/1‑149 (MQ=255)                                                                                                                                                  ggCATGGATGCGCGTCACCTCCCTTCGCAAGCCGTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtgt    <  7:907/149‑1 (MQ=255)                                                                                                                                                  ggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtgt    <  2:243282/149‑1 (MQ=255)                                                                                                                                                  ggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtgt    <  7:135349/149‑1 (MQ=255)                                                                                                                                                  ggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACtgt    <  6:202901/149‑1 (MQ=255)                                                                                                                                                    cATCCATGCGCGTCACCTCCCTTCGCAAGCCCTGCGCTAGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACTGTGc  <  6:219281/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  TGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCT‑CGCTACGGTCTTGCTCCCGGGA‑‑GTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAGATACAGAAACAGTGTTACTGTGC  >  NC_002947/4909802‑4910094

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑