Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,621,030	(A)10→9	coding (40/417 nt)	PP_5654 →	phosphotriesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,621,021	0	A	.	90.7%	277.8 / 22.7	75	coding (31/417 nt)	PP_5654	phosphotriesterase
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base . (33/35);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAAG  >  NC_002947/4620875‑4621164                                                                                                                                                     |                                                                                                                                                 aaCTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCaa                                                                                                                                                 >  7:103076/1‑148 (MQ=255)    tctctcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGCAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GACCaaaaa                                                                                                                                              >  3:176740/1‑148 (MQ=255)       ctcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCaaaaaaaa                                                                                                                                           >  3:254597/1‑148 (MQ=255)        tcGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAATATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCaaaaaaaa                                                                                                                                           <  1:225760/147‑1 (MQ=255)          gACCTTCCAAGAATCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCaaaaaaaaa                                                                                                                                          <  5:245937/146‑1 (MQ=255)          gACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCaaaaaaaaa                                                                                                                                          <  5:255107/146‑1 (MQ=255)              ttCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGAc                                                                                                                                  >  7:161204/1‑149 (MQ=255)                 cAAGAAGCTAAGTTAGAGCGAATTAATAAAAAAAAGCAAACAATAAAGCCCAATCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGTAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑c                                                                                                                              <  8:298247/149‑1 (MQ=255)                  aagaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAACAAAGACTa                                                                                                                                >  1:150116/1‑148 (MQ=255)                  aagaagCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑Ca                                                                                                                             >  6:134556/1‑149 (MQ=255)                       gCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAA‑AAAGACTA‑CATTaaa                                                                                                                        >  3:286870/1‑149 (MQ=255)                       gCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTaaa                                                                                                                        <  7:22108/149‑1 (MQ=255)                         gAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAg                                                                                                                       >  2:249677/1‑148 (MQ=255)                          aaGTTAGAGAGAATTATTTAAAATAAGCAAACAATAAAGACCAAGCAAAACTCACACAAACTCCCGTAAAACTCATAAGTGGAACATGCAAATGAACTTAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGaa                                                                                                                     <  1:289747/149‑1 (MQ=255)                             ttAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACa                                                                                                                   >  7:305676/1‑148 (MQ=255)                               agagCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCACGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATa                                                                                                                 <  1:86779/148‑1 (MQ=255)                                  gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATaaaa                                                                                                              <  1:249677/148‑1 (MQ=255)                                  gCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAg                                                                                                             >  4:281887/1‑149 (MQ=255)                                   cGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACACACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAACCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGc                                                                                                            <  6:103981/149‑1 (MQ=255)                                   cGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGc                                                                                                            <  7:34481/149‑1 (MQ=255)                                      attattTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAgg                                                                                                         <  3:281887/149‑1 (MQ=255)                                       ttattTAAAATAAGCAAACAATAAATCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGga                                                                                                        <  4:254597/149‑1 (MQ=255)                                         attAAAAATAACCAAAAAAAAAACCCCACGCAAAACACACACAAACTCCCGTAAAACTCAGAAGTGTAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGgaga                                                                                                      <  8:303806/149‑1 (MQ=255)                                          tttAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGgaga                                                                                                      >  3:220536/1‑148 (MQ=255)                                             aaaaTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATAtt                                                                                                  >  8:12956/1‑149 (MQ=255)                                                aTAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGAAGGAGATATTTTc                                                                                               >  4:111435/1‑149 (MQ=255)                                                     cAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGaa                                                                                                                     <  8:328610/122‑1 (MQ=255)                                                     cAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGaa                                                                                                                     >  7:328610/1‑122 (MQ=255)                                                      aaaCAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTc                                                                                         <  8:305676/149‑1 (MQ=255)                                                          aaTAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCAcc                                                                                      >  1:307796/1‑148 (MQ=255)                                                           aTAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCa                                                                                    <  8:103076/149‑1 (MQ=255)                                                             aaaGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTGCGAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑Ca                                                                                                                             <  1:35413/107‑1 (MQ=255)                                                             aaaGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAA‑AAAGACTA‑Ca                                                                                                                             >  2:35413/1‑106 (MQ=255)                                                              aaGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAgg                                                                                  <  2:10698/148‑1 (MQ=255)                                                                gCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTATGGTTACCCCAgg                                                                                  >  6:230596/1‑146 (MQ=255)                                                                 cccAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGa                                                                              >  7:29142/1‑149 (MQ=255)                                                                 cccAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGa                                                                              >  6:160231/1‑149 (MQ=255)                                                                 cccAAGCAAAACTCACACAAACTCCCGTAAAAATCAGAAGTGGAACATGCAAATGAAAAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGAATA‑CATTAAAGAACATAAAAGCAGGAGAAAATTTCTGCTTAACCCAGGAAGa                                                                              >  6:370061/1‑149 (MQ=255)                                                                   cAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACg                                                                            >  1:165681/1‑149 (MQ=255)                                                                      gCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGaa                                                                         <  7:57451/149‑1 (MQ=255)                                                                        aaaaCTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAAc                                                                        <  2:150116/149‑1 (MQ=255)                                                                        aaaaCTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACa                                                                       <  7:12956/149‑1 (MQ=255)                                                                          aaCTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATTCAAATGAACAAACCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAAGAACGAGTA‑CATTAAAGAACATATTAGCAGGTGATTTTTTCTGCTTCACCCAGGAAGACGGAACata                                                                     <  8:230725/149‑1 (MQ=255)                                                                                 acaAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTAtttt                                                              >  8:265192/1‑149 (MQ=255)                                                                                   aaaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                            >  7:26280/1‑149 (MQ=255)                                                                                   aaaCTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTgg                                                            <  2:165681/149‑1 (MQ=255)                                                                                         ccGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCaaa                                                       >  1:32919/1‑148 (MQ=255)                                                                                              aaaCTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCATGAAGACGGAACATATTATTTTGGGAAAataata                                                  >  5:232999/1‑148 (MQ=255)                                                                                               aaCTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGAAAAGCATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAAAATATTATTTTGGGCAAataata                                                  >  1:119861/1‑148 (MQ=255)                                                                                                 cTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCa                                               >  1:23685/1‑148 (MQ=255)                                                                                                  tCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATGCGTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCa                                               >  1:274986/1‑148 (MQ=255)                                                                                                      aaGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAaat                                          >  6:387111/1‑148 (MQ=255)                                                                                                      aaGTGGAACATGAAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAaata                                         <  3:111435/149‑1 (MQ=255)                                                                                                       aGTGGAACATTCAAATTACCAAAGCTTAGATATG‑GAGATG‑TTCC‑AAAAAAATACAAAGACTA‑CATTAATGAACATAATAGGAGGAGATATTTTCTGCTTCACCGAGTAATACTGTTCTTATTATTTTGGGCAAATAATATCaaaaaaaa                                        <  4:176740/149‑1 (MQ=255)                                                                                                        gTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGGTAtttt                                                                                                >  2:153462/1‑92 (MQ=255)                                                                                                        gTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGGTAtttt                                                                                                <  1:153462/92‑1 (MQ=255)                                                                                                        gTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAaataa                                        >  3:354980/1‑148 (MQ=255)                                                                                                         tGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATGTAAGAACATAAAAGCCGGAGATATTTTATGCATAACCCAGGAAGAAGGGACATTTTATTTTTGGCAAATAATATAAAAaattat                                       >  1:309564/1‑148 (MQ=255)                                                                                                           gAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCAc                                    <  2:32919/149‑1 (MQ=255)                                                                                                            aaCATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCTAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCAc                                    <  3:13068/148‑1 (MQ=255)                                                                                                                tGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGcc                               <  2:274986/149‑1 (MQ=255)                                                                                                                  cAAATGAAAAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAACAACATAAAAGCAGGAGATATTTTCTGCTTCACGCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCa                              <  5:370061/148‑1 (MQ=255)                                                                                                                   aaaTGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                            >  8:15450/1‑149 (MQ=255)                                                                                                                   aaaTGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                            >  5:35600/1‑149 (MQ=255)                                                                                                                    aaTGAACAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATg                            >  8:268578/1‑148 (MQ=255)                                                                                                                        aaCAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCt                       >  1:365143/1‑149 (MQ=255)                                                                                                                          cAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAATCATGAGATATTTTCTTATTCACACAAGAAGACGGAACATAATATTTTTGGCAAATAATAAGAAAATTAATAACAGGCCATGAAgat                       >  4:236866/1‑146 (MQ=255)                                                                                                                          cAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAataata                                                  >  3:293751/1‑120 (MQ=255)                                                                                                                          cAAAGCTAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAataata                                                  <  4:293751/120‑1 (MQ=255)                                                                                                                              gatAACATCTG‑GAGATT‑TTAAAAAAAAAAAAAAAATACTA‑CATTAAACAACATAAAAGAATGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAat                  <  3:236866/147‑1 (MQ=255)                                                                                                                                tAACATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATAtt               <  4:220536/149‑1 (MQ=255)                                                                                                                                 aaCATCTG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGAAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATAtt               >  5:81855/1‑148 (MQ=255)                                                                                                                                   cATCTG‑GAGATT‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAATAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAg            <  5:230596/149‑1 (MQ=255)                                                                                                                                       tG‑GAGATG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATg        >  3:64709/1‑149 (MQ=255)                                                                                                                                          gagaTG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGaa          <  4:286870/145‑1 (MQ=255)                                                                                                                                          gagaTG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAAt         <  2:309564/146‑1 (MQ=255)                                                                                                                                            gaTG‑GT‑CAAAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTa     >  5:176365/1‑148 (MQ=255)                                                                                                                                             aTG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa   <  4:354980/149‑1 (MQ=255)                                                                                                                                             aTG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa   <  2:365143/149‑1 (MQ=255)                                                                                                                                              tG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaag  <  5:134556/149‑1 (MQ=255)                                                                                                                                              tG‑GTCC‑AAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTaaa   <  7:15450/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AACTCTCTCGACCTTCCAAGAAGCGAAGTTAGAGCGAATTATTTAAAATAAGCAAACAATAAAGCCCAAGCAAAACTCACACAAACTCCCGTAAAACTCAGAAGTGGAACATGCAAATGAACAAAGCTAACATCTG‑GAGATG‑GTCCAAAAAAAAAACAAAGACTA‑CATTAAAGAACATAAAAGCAGGAGATATTTTCTGCTTCACCCAGGAAGACGGAACATATTATTTTGGGCAAATAATATCAAAAATAATCACTGGCCATGTAGCTGAAATATTCAGAATGGTAAAG  >  NC_002947/4620875‑4621164

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