Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,456,854	Δ1 bp	6.8%	intergenic (‑42/+161)	sixA ← / ← fadJ	phosphohistidine phosphatase/3‑hydroxyacyl‑CoA dehydrogenase FadJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,456,854	0	A	.	6.8%	235.9 / 10.1	73	intergenic (‑42/+161)	sixA/fadJ	phosphohistidine phosphatase/3‑hydroxyacyl‑CoA dehydrogenase FadJ
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base . (2/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CCATCAGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAG  >  NC_000913/2456710‑2457003                                                                                                                                                 |                                                                                                                                                       ccATCAGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGa                                                                                                                                                   <  1:537519/150‑1 (MQ=255)     cAGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAAcc                                                                                                                                               >  1:911478/1‑150 (MQ=255)       ggCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCag                                                                                                                                             >  2:1317199/1‑150 (MQ=255)        gCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCaga                                                                                                                                            <  1:527288/150‑1 (MQ=255)         cGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCaga‑g                                                                                                                                          <  2:546058/150‑1 (MQ=255)         cGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCaga‑g                                                                                                                                          <  2:173953/150‑1 (MQ=255)             aTTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGa‑gaaag                                                                                                                                      <  1:222674/150‑1 (MQ=255)              ttCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGa‑gaaaga                                                                                                                                     <  1:1145448/150‑1 (MQ=255)                  tCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAAc                                                                                                                                 >  2:1029813/1‑150 (MQ=255)                     caACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTC‑TCTGGTTTCAGATGAAAGAAAACggg                                                                                                                              <  2:1416925/150‑2 (MQ=38)                     caACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTC‑TCTGGTTTCAGATGAAAGAAAACggg                                                                                                                              <  1:864707/150‑2 (MQ=38)                       aCCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTgg                                                                                                                            <  2:188356/150‑1 (MQ=255)                        ccATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGc                                                                                                                           >  1:559876/1‑150 (MQ=255)                          aTTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCg                                                                                                                         >  1:782156/1‑150 (MQ=255)                             aGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTC‑TCTGGTTTCAGATGAAAGAAAACggg                                                                                                                              >  2:1392481/1‑141 (MQ=38)                             aGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTC‑TCTGGTTTCAGATGAAAGAAAACggg                                                                                                                              <  1:1392481/142‑2 (MQ=38)                                ggTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAAATTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCAttt                                                                                                                   >  2:95860/1‑149 (MQ=255)                                  tCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATtgtg                                                                                                                 <  2:1268193/150‑1 (MQ=255)                                   cAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGc                                                                                                                >  2:344874/1‑150 (MQ=255)                                     gagGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCt                                                                                                              >  2:159612/1‑150 (MQ=255)                                     gagGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCt                                                                                                              <  1:1365126/150‑1 (MQ=255)                                      agGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCtt                                                                                                             >  1:805595/1‑150 (MQ=255)                                       ggACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTa                                                                                                            >  1:164939/1‑150 (MQ=255)                                         aCGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATc                                                                                                          >  2:1246458/1‑150 (MQ=255)                                         aCGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATc                                                                                                          <  1:1398491/150‑1 (MQ=255)                                          cGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATcc                                                                                                         <  1:687248/150‑1 (MQ=255)                                          cGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATcc                                                                                                         >  1:804190/1‑150 (MQ=255)                                            aaCGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCAt                                                                                                       <  2:316140/150‑1 (MQ=255)                                               ggAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCa                                                                                                    <  1:1151316/150‑1 (MQ=255)                                               ggAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCa                                                                                                    <  2:658467/150‑1 (MQ=255)                                                  aTCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCg                                                                                                 <  1:428772/150‑1 (MQ=255)                                                   tCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGa                                                                                                >  2:23556/1‑150 (MQ=255)                                                     aCTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAAt                                                                                              <  1:778762/150‑1 (MQ=255)                                                      cTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATg                                                                                             <  1:1317199/150‑1 (MQ=255)                                                        ggcggcATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGaa                                                                                           >  2:393618/1‑150 (MQ=255)                                                         gcggcATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGaaa                                                                                          <  1:446996/150‑1 (MQ=255)                                                          cggcATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAAc                                                                                         >  1:1286998/1‑150 (MQ=255)                                                                 gAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGttt                                                                                  >  2:1127226/1‑150 (MQ=255)                                                                  aGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGtttt                                                                                 <  2:911478/150‑1 (MQ=255)                                                                   gggCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGttttt                                                                                >  1:1064783/1‑150 (MQ=255)                                                                     gCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTa                                                                              <  1:1029813/150‑1 (MQ=255)                                                                       tGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTAcc                                                                            <  1:393618/150‑1 (MQ=255)                                                                         cGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTc                                                                          <  2:1210610/150‑1 (MQ=255)                                                                             gCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATg                                                                      >  2:163427/1‑150 (MQ=255)                                                                             gCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATg                                                                      >  2:848547/1‑150 (MQ=255)                                                                              ccGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATgg                                                                     <  1:621408/150‑1 (MQ=255)                                                                                 tGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTCATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGt                                                                  >  2:1065976/1‑150 (MQ=255)                                                                                     gCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGt                                                              >  1:970022/1‑150 (MQ=255)                                                                                      cATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATg                                                                                             >  2:359459/1‑118 (MQ=255)                                                                                      cATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATg                                                                                             <  1:359459/118‑1 (MQ=255)                                                                                           tAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTc                                                        <  2:340302/150‑1 (MQ=255)                                                                                            aaaaaCTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCa                                                       >  2:340597/1‑150 (MQ=255)                                                                                             aaaaCTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCaa                                                      <  1:1127226/150‑1 (MQ=255)                                                                                                 cTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATg                                                  <  2:782156/150‑1 (MQ=255)                                                                                                 cTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATg                                                  >  1:1316691/1‑150 (MQ=255)                                                                                                   tGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCg                                                >  2:221856/1‑150 (MQ=255)                                                                                                      atatTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGttt                                             <  2:119895/150‑1 (MQ=255)                                                                                                            cACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACATCAATGGCGTAAGTATAGTCAATATGCGTTTACAttt                                       >  2:716699/1‑150 (MQ=255)                                                                                                             aCCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACAtttt                                      <  2:164939/150‑1 (MQ=255)                                                                                                             aCCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACAtttt                                      <  1:944690/150‑1 (MQ=255)                                                                                                              ccGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACAttttt                                     >  1:949769/1‑149 (MQ=255)                                                                                                              ccGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTg                                     <  2:1064783/150‑1 (MQ=255)                                                                                                              ccGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTg                                     <  1:734811/150‑1 (MQ=255)                                                                                                               cGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGa                                    <  2:194521/150‑1 (MQ=255)                                                                                                                gCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGACACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATAGGCGTTTACATTTTGAc                                   >  2:636169/1‑150 (MQ=255)                                                                                                                gCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGAc                                   >  2:136687/1‑150 (MQ=255)                                                                                                                             aGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGGGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCa                      >  2:952799/1‑150 (MQ=255)                                                                                                                               aTTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCAtt                    >  1:108560/1‑150 (MQ=255)                                                                                                                                   caccGTTTTCTTTC‑TCTGGTTTCAGATGAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAgc                >  1:219918/3‑150 (MQ=38)                                                                                                                                     ccGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAgcgc              >  1:64007/1‑150 (MQ=255)                                                                                                                                     ccGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAgcgc              <  2:294511/150‑1 (MQ=255)                                                                                                                                     ccGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAgcgc              >  2:1005613/1‑150 (MQ=255)                                                                                                                                                cATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCAtata   <  1:340597/150‑1 (MQ=255)                                                                                                                                                 aTCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGGTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAg  >  2:29443/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                       CCATCAGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATATTGCACCGCTTTTGTTAACCAGATTCGCCCGTTTTCTTTCATCTGAAACCAGA‑GAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAG  >  NC_000913/2456710‑2457003

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