Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,624,544	T→C	6.9%	intergenic (+3137/+1427)	PP_5654 → / ← PP_4091	phosphotriesterase/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,624,544	0	T	C	6.9%	200.6 / 5.1	73	intergenic (+3137/+1427)	PP_5654/PP_4091	phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base C (3/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

CAAATTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTC  >  NC_002947/4624397‑4624685                                                                                                                                                     |                                                                                                                                              cAAATTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGAAATGACAAACTAGCTAGCAGGGAAATTCt                                                                                                                                               >  6:62307/1‑148 (MQ=255)  aaaTTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTaa                                                                                                                                             >  2:202838/1‑149 (MQ=255)  aaaTTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTaa                                                                                                                                             >  2:241009/1‑149 (MQ=255)           gTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACA‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTg                                                                                                                                    >  7:43015/1‑149 (MQ=255)               tAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGttt                                                                                                                                 <  4:296255/148‑1 (MQ=255)                aGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTc                                                                                                                                >  5:137101/1‑148 (MQ=255)                 gATTATGTCCCCTTAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCaa                                                                                                                              <  5:206964/149‑1 (MQ=255)                  ataaagTCCCCAGAGAAAAAAATAATCGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAg                                                                                                                             <  6:306129/144‑1 (MQ=255)                  aTTATGTCCCCTGAGAAAAAAATAATTGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCaa                                                                                                                              <  3:212643/148‑1 (MQ=255)                   ttATGTCCCCTGAGAAAAAAATAATTGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGc                                                                                                                            <  4:303222/149‑1 (MQ=255)                         ccccTGAGAAAAAGATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaa                                                                                                                       >  8:228545/1‑148 (MQ=255)                          cccTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGGTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaac                                                                                                                      >  3:56296/1‑148 (MQ=255)                          cccTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaac                                                                                                                      >  6:51547/1‑148 (MQ=255)                          cccTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCAAAAACC‑TCATGACATTGAAATGCTGAATAGAGAAATTTACAAATGGCTTTGAAATAAAAAAAAAGAAAGCAGGAAAATAAAAATGAAATTTTTTAAAAGcttaaaa                                                                                                                      >  8:239692/1‑142 (MQ=255)                           ccTGAGAAAAAAATAATGGCTACCATTAGAAAACGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTaacaa                                                                                                                    <  4:86043/149‑1 (MQ=255)                            cTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACCAGCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGTAATTTGTTTCAAGCGTTaaca                                                                                                                     >  1:281157/1‑148 (MQ=255)                               agaAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACt                                                                                                                  <  2:186635/147‑1 (MQ=255)                                gaAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCa                                                                                                                <  5:62307/148‑1 (MQ=255)                                   aaaatTAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACAATGAAATGCTGCCTAGAGAAAGTTTAACTTGGCTTTCCAAAGACAATCTAGCTAGCAGGGAATTTGTAATGGATTTTGTTTCAAGCGTTAACAACTCAAAAt                                                                                                            >  2:89545/1‑149 (MQ=255)                                    aaaaTAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAAtt                                                                                                           >  8:245711/1‑149 (MQ=255)                                        taatGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTATCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCt                                                                                                       <  8:258443/149‑1 (MQ=255)                                              cTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTTGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGt                                                                                                  >  3:358513/1‑148 (MQ=255)                                                 ccATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCAt                                                                                               <  7:245711/148‑1 (MQ=255)                                                      aGAAACCGATACCGAAGACC‑GCATGACATTGAAATGCTGGCTCGCGAACTTTCCACATGGCTTTGCAATGACAAACTAGCTCGCAGGTAAATTCTAATGGAATTTGTTTCAAGCGTTAACTACTCAAAATTCCCTGACGTCTTTCTAc                                                                                          <  4:269732/148‑1 (MQ=255)                                                        aaaaGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTaa                                                                                       >  7:285541/1‑149 (MQ=255)                                                        aaaaGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTaa                                                                                       >  4:240956/1‑149 (MQ=255)                                                        aaaaGGATACCGAAGACC‑TAATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCACTGACGTCATTCAACTaa                                                                                       >  8:49998/1‑149 (MQ=255)                                                         aaaGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAAc                                                                                      <  4:282192/149‑1 (MQ=255)                                                             gATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGATATTTACTCATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCttt                                                                                  <  4:56296/149‑1 (MQ=255)                                                              aTACCGAAGACC‑TAATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCttt                                                                                  <  7:239692/148‑1 (MQ=255)                                                               tACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTg                                                                                 >  4:346372/1‑148 (MQ=255)                                                               tACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                <  5:51547/149‑1 (MQ=255)                                                               tACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                <  3:346372/149‑1 (MQ=255)                                                               tACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGa                                                                                <  8:285541/149‑1 (MQ=255)                                                                  cGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGtt                                                                                                                         >  4:304843/1‑105 (MQ=255)                                                                  cGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGtt                                                                                                                         <  3:304843/105‑1 (MQ=255)                                                                    aaGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACt                                                                           <  2:80875/149‑1 (MQ=255)                                                                      gACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTg                                                                         <  1:189454/149‑1 (MQ=255)                                                                         c‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAg                                                                      >  8:185119/1‑149 (MQ=255)                                                                                  aTTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATc                                                              >  4:9106/1‑149 (MQ=255)                                                                                       aaTGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTc                                                         <  7:49998/149‑1 (MQ=255)                                                                                        aTGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCa                                                        <  4:279410/149‑1 (MQ=255)                                                                                            tGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGa                                                     >  5:34020/1‑148 (MQ=255)                                                                                            tGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGa                                                     >  5:28435/1‑148 (MQ=255)                                                                                               cTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGTGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGc                                                 <  2:281157/149‑1 (MQ=255)                                                                                                    gAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGt                                             <  1:241009/148‑1 (MQ=255)                                                                                                    gAAATTTACACATGGATTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCGCTGAAGTCATTCAACTAACCTTTGAATACTATAAGCTCCTATCCAGTCATGAAACCTAGt                                             >  4:313360/1‑148 (MQ=255)                                                                                                     aaaTTTACACATGGCGTTGGAAAGACAAACTAGCTAGCAGGGAAATTGTTATGGAATTTGTTTGAAGCGGTTACAACTCAAAATTACCGGACGTCATTCAACTTACCTTTGAATACTTGAAGCGCCTTTCCACTACAGCAAGCGAGttt                                           >  4:368466/1‑148 (MQ=255)                                                                                                       aTTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACt                                         <  7:185119/149‑1 (MQ=255)                                                                                                         tgagggaTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACtat                                       <  1:79198/143‑1 (MQ=255)                                                                                                               aTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACTTCATTCAACTAAACTTTGAATACTTTAAGCGACTATCCACTCATGAAAGAGAGTTACTATACgat                                  >  5:389461/1‑147 (MQ=255)                                                                                                               aTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGAGCCTATCCACTCATGAAAGCGAGTTACTATACgag                                  >  4:368529/1‑148 (MQ=255)                                                                                                               aTGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTAATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTAATATACgag                                  >  5:195568/1‑148 (MQ=255)                                                                                                               aTGGC‑TTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgag                                  >  6:69262/1‑147 (MQ=255)                                                                                                                tGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgag                                  >  6:370172/1‑147 (MQ=255)                                                                                                                  gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagga                                >  6:290681/1‑147 (MQ=255)                                                                                                                  gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagga                                <  5:290681/147‑1 (MQ=255)                                                                                                                  gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGgaga                              <  1:202838/149‑1 (MQ=255)                                                                                                                  gCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGgaga                              <  5:370172/149‑1 (MQ=255)                                                                                                                   gcttGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGgaga                              <  5:69262/146‑1 (MQ=255)                                                                                                                    tttGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACtt                                                                          >  5:329486/1‑103 (MQ=255)                                                                                                                    tttGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACtt                                                                          <  6:329486/103‑1 (MQ=255)                                                                                                                               aaaCTAGCTACCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAgg                 >  4:80045/1‑149 (MQ=255)                                                                                                                                     gatagcAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCTCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCTaaaa                       >  6:286685/3‑137 (MQ=255)                                                                                                                                        ttgaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCatttatt        >  5:58609/8‑149 (MQ=255)                                                                                                                                         gcAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCatttatt        <  7:228545/148‑1 (MQ=255)                                                                                                                                          gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                 <  4:241516/117‑1 (MQ=255)                                                                                                                                          gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                 <  4:236662/117‑1 (MQ=255)                                                                                                                                          gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                 >  3:241516/6‑122 (MQ=255)                                                                                                                                          gaaacaaaTTCCATTAGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACgagg                                 >  3:236662/6‑122 (MQ=255)                                                                                                                                            gggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGaa                         >  7:171642/1‑128 (MQ=255)                                                                                                                                            gggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGaa                         <  8:171642/128‑1 (MQ=255)                                                                                                                                             ggAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGAtt   <  4:358513/149‑1 (MQ=255)                                                                                                                                              gAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTc  <  6:137101/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              CAAATTAGAGGTTTTAGATTATGTCCCCTGAGAAAAAAATAATGGCTACCATTAGAAAAGGATACCGAAGACC‑TCATGACATTGAAATGCTGGCTAGCGAAATTTACACATGGCTTTGCAATGACAAACTAGCTAGCAGGGAAATTCTAATGGAATTTGTTTCAAGCGTTAACAACTCAAAATTCCCTGACGTCATTCAACTAACCTTTGAATACTTGAAGCGCCTATCCACTCATGAAAGCGAGTTACTATACGAGGAGAGCGAAAAAATAGGGCATTTATTTGATTC  >  NC_002947/4624397‑4624685

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