Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,358,604	A→T	5.8%	L451Q (CTG→CAG)	PP_3834 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,358,604	0	A	T	5.8%	213.4 / 2.8	69	L451Q (CTG→CAG)	PP_3834	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base T (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

CAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGATTTC  >  NC_002947/4358463‑4358749                                                                                                                                              |                                                                                                                                                   cAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTg                                                                                                                                             >  6:370883/1‑149 (MQ=255)   ggTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGc                                                                                                                                            >  7:216130/1‑148 (MQ=255)    gTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGc                                                                                                                                            >  1:80779/1‑147 (MQ=255)      aTTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAgg                                                                                                                                        >  3:37756/1‑149 (MQ=255)         cGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCg                                                                                                                                     >  2:110388/1‑149 (MQ=255)          gTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCgg                                                                                                                                    <  3:181218/149‑1 (MQ=255)            cATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCgg                                                                                                                                    <  3:179103/147‑1 (MQ=255)                gACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCTGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCGTGCGCAAGGTAATGCCTAGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATg                                                                                                                              >  6:244991/1‑149 (MQ=255)                  cAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGt                                                                                                                            >  6:279289/1‑149 (MQ=255)                  cAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGt                                                                                                                            <  5:145128/149‑1 (MQ=255)                     tCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGAc                                                                                                                         >  2:191346/1‑149 (MQ=255)                     tCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCACAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGAc                                                                                                                         >  3:21757/1‑149 (MQ=255)                         gATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTc                                                                                                                     <  5:35406/149‑1 (MQ=255)                             gtgtTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACg                                                                                                                  >  2:248869/1‑148 (MQ=255)                              tgtTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCa                                                                                                                >  2:325990/1‑149 (MQ=255)                                ttCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCAcc                                                                                                              >  5:280956/1‑149 (MQ=255)                                ttCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCAc                                                                                                               <  5:370883/148‑1 (MQ=255)                                 tCAGCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGGCGGCGTTCCCTAAAGCACACGGCAGGAACCTTCTAAGGCGTTGCGATCGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCAcc                                                                                                              <  1:110388/148‑1 (MQ=255)                                   aTCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCtt                                                                                                            <  7:211693/148‑1 (MQ=255)                                        cTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGAc                                                                                                      <  6:139139/149‑1 (MQ=255)                                            tccgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGa                                                                                                   >  4:327924/1‑148 (MQ=255)                                             ccgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCaa                                                                                                                                          >  6:313525/1‑108 (MQ=255)                                             ccgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCaa                                                                                                                                          <  5:313525/108‑1 (MQ=255)                                              cgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCgg                                                                                                                                    <  2:169369/113‑1 (MQ=255)                                               gtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAg                                                                                                >  8:338976/1‑148 (MQ=255)                                                  aCCCGCTCTCCTGCGTAGTAGTCGTTCCCGAAAGCACACCGCAGGCACATTCGAACGCCTTGCGCAATGTAATGCCTTTAGTCACGATCAGCATCTCTTTAAAGGCAGGGTCATGGTTACTGTCACGCACCTTTATTACCGACATCATg                                                                                            <  5:76085/149‑1 (MQ=255)                                                    ccGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATcagca                                                                                                                                                    >  6:226065/1‑91 (MQ=255)                                                    ccGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATcagca                                                                                                                                                    <  5:226065/91‑1 (MQ=255)                                                     cGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATgcgc                                                                                         <  1:325990/149‑1 (MQ=255)                                                       ctctCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCg                                                                                       <  8:216130/149‑1 (MQ=255)                                                            cTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGtt                                                                                  <  1:248869/149‑1 (MQ=255)                                                                ggcgtcggcgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGa                                                                              <  6:280956/149‑1 (MQ=255)                                                                   gtcggcgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTAGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTc                                                                           <  5:244991/149‑1 (MQ=255)                                                                       gcgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGt                                                                       >  8:244870/1‑149 (MQ=255)                                                                        cgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGtt                                                                      <  5:279289/149‑1 (MQ=255)                                                                             ccGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTc                                                                           <  2:52381/139‑1 (MQ=255)                                                                             ccGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTc                                                                           >  1:52381/1‑139 (MQ=255)                                                                             ccGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGAc                                                                  >  6:43034/1‑148 (MQ=255)                                                                               gAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACa                                                                 >  3:28061/1‑147 (MQ=255)                                                                               gAAAGCACACCGCAGACACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACaaa                                                               <  1:11913/149‑1 (MQ=255)                                                                                     acacCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCAt                                                         >  3:204668/1‑149 (MQ=255)                                                                                           cAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCa                                                    >  5:311328/1‑148 (MQ=255)                                                                                            aGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATa                                                  <  3:327924/149‑1 (MQ=255)                                                                                              gCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCAt                                                         >  7:339565/1‑140 (MQ=255)                                                                                              gCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCAt                                                         <  8:339565/140‑1 (MQ=255)                                                                                                aCCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGt                                              >  6:193651/1‑149 (MQ=255)                                                                                                  cTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGc                                                                                <  6:355878/113‑1 (MQ=255)                                                                                                  cTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGc                                                                                >  5:355878/1‑113 (MQ=255)                                                                                                   ttCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAATGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGt                                           >  5:347624/1‑149 (MQ=255)                                                                                                         cGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cgacga                                    <  5:2663/149‑1 (MQ=255)                                                                                                          gCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cgacga                                    >  5:24764/1‑148 (MQ=255)                                                                                                            cTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATAc                                 >  3:68414/1‑149 (MQ=255)                                                                                                             ttGCGGAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTAACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCTCATCAGCATAATGTTGT‑cgacgcttct                                >  3:191391/1‑144 (MQ=255)                                                                                                               gcgcAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAACGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACTATACTc                               <  6:347624/148‑1 (MQ=255)                                                                                                                cgcAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGa                             <  4:68414/149‑1 (MQ=255)                                                                                                                          tGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGgc                   >  3:330120/1‑149 (MQ=255)                                                                                                                             cTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATTCGCCGTAGTTGCGAGTCAAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGagcg                >  2:174/1‑149 (MQ=255)                                                                                                                             cTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGcgcg                >  3:68311/1‑149 (MQ=255)                                                                                                                             cTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGcgcg                >  4:297340/1‑149 (MQ=255)                                                                                                                                gCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGTCCGACGATACTCGACAGGTCCGGCGCgg               <  1:191346/148‑1 (MQ=255)                                                                                                                                gCGGCACGATCACCAGCTCGTGCAAGGCCGGGTCATGGGGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCggg              <  4:37756/148‑1 (MQ=255)                                                                                                                                 ttgCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg             >  7:59157/3‑148 (MQ=255)                                                                                                                                 cGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg             <  6:193778/148‑1 (MQ=255)                                                                                                                                 cGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg             <  5:43034/148‑1 (MQ=255)                                                                                                                                  ggCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGgc           <  2:80779/149‑1 (MQ=255)                                                                                                                                   gCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCTGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGgcg          >  5:189434/1‑149 (MQ=255)                                                                                                                                   gCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cga                                       >  6:286927/1‑120 (MQ=255)                                                                                                                                   gCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cga                                       <  5:286927/120‑1 (MQ=255)                                                                                                                                   gCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg             <  8:59157/146‑1 (MQ=255)                                                                                                                                      cGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCTGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTg       <  6:189434/149‑1 (MQ=255)                                                                                                                                         tCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGAtt    >  4:351755/1‑149 (MQ=255)                                                                                                                                           agcagcTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCCACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGATTTc  <  7:338976/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   CAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGATTTC  >  NC_002947/4358463‑4358749

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