Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,009,770	T→C	6.3%	W201R (TGG→CGG) ‡	pldB →	lysophospholipase L2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,770	0	T	C	6.3%	238.5 / 4.7	80	W201R (TGG→CGG) ‡	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base C (3/2);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

CCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGA  >  NC_000913/4009622‑4009915                                                                                                                                                     |                                                                                                                                                  ccATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCt                                                                                                                                                   <  1:1484507/149‑1 (MQ=255)   aTCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTgg                                                                                                                                                 >  2:14344/1‑149 (MQ=255)       aGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGgcgcg                                                                                                                                             >  1:1581956/1‑149 (MQ=255)          tgtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGt                                                                                                                                          >  1:157354/1‑149 (MQ=255)          ggtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGt                                                                                                                                          <  2:231218/148‑1 (MQ=255)              tgtgACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGcc                                                                                                                                      >  2:1012252/1‑149 (MQ=255)              tgtgACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGAAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCCCCGGGCGCTGGCGCGCGGTGcc                                                                                                                                      >  1:112357/1‑149 (MQ=255)                tgACGCCATTGCTCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAAACGGCAACGCGCGCCAGCGCCCGGTGccta                                                                                                                                    <  1:956004/149‑3 (MQ=255)                tgACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGt                                                                                                                                    >  1:1387768/1‑149 (MQ=255)                  aCGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGttt                                                                                                                                  <  2:826633/149‑1 (MQ=255)                      cATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTa                                                                                                                              >  2:1459385/1‑149 (MQ=255)                      cATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTa                                                                                                                              <  2:1159263/149‑1 (MQ=255)                       atggcgcTAACGGCGCCAATGGTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTAt                                                                                                                             <  1:1011953/146‑1 (MQ=255)                       aTTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCAACGCGCGCCAGCGCCCGGTGCCTATTGcata                                                                                                                             <  2:1040294/149‑4 (MQ=38)                       aTTGCGCGAACTGCGCCAATGTTGGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCGGGCGCGCGTTGCCGTTTGCTAt                                                                                                                             <  2:159758/149‑1 (MQ=255)                        ttGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTAta                                                                                                                            >  2:618369/1‑148 (MQ=255)                           cgcTAACTGCTCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACTTTTCCGTGATGGATATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAc                                                                                                                         <  1:951908/149‑1 (MQ=255)                           cgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAc                                                                                                                         >  1:364259/1‑149 (MQ=255)                           cgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAc                                                                                                                         >  1:1526207/1‑149 (MQ=255)                           cgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAc                                                                                                                         <  2:1278013/149‑1 (MQ=255)                              tAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTa                                                                                                                      >  2:391543/1‑149 (MQ=255)                                aCTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACt                                                                                                                    >  2:36645/1‑149 (MQ=255)                                  tGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGa                                                                                                                  <  1:187770/149‑1 (MQ=255)                                  tGCGCCAATGTGTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGa                                                                                                                  <  2:1231012/149‑1 (MQ=255)                                    cgcCAATGTTTAGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGAcc                                                                                                                >  2:601112/1‑149 (MQ=255)                                     gcCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGAccc                                                                                                               >  2:1047765/1‑149 (MQ=255)                                     gcCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGAccc                                                                                                               <  1:1317022/149‑1 (MQ=255)                                       cAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCcac                                                                                                             <  1:145090/149‑1 (MQ=255)                                            tttGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCAcagcag                                                                                                        <  1:1335664/149‑1 (MQ=255)                                             ttGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGa                                                                                                       <  1:1268544/149‑1 (MQ=255)                                                   tCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGa                                                                                                 <  1:1012252/149‑1 (MQ=255)                                                     gTGATTCGTATGCCGTCATTTATGGCACGGCAGATCATCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACCTACTGACCCACAGCAGACAGCGata                                                                                               >  2:441255/1‑149 (MQ=255)                                                      ggATTCGTATGCCGTCATTTATGGCTCGGCAGATCCTCTACTGGGGCGTAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGatat                                                                                              <  1:257293/148‑1 (MQ=255)                                                       gATTCGTATGCCGTCATTTATGGCACGGAAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATc                                                                                             <  2:571310/149‑1 (MQ=255)                                                        aTTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCg                                                                                            >  2:59731/1‑149 (MQ=255)                                                        aTTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCg                                                                                            >  1:1278696/1‑149 (MQ=255)                                                                gCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACt                                                                                    <  1:618369/149‑1 (MQ=255)                                                                   gTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTAc                                                                                 >  1:419508/1‑149 (MQ=255)                                                                   gTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTAc                                                                                 <  1:620116/149‑1 (MQ=255)                                                                       tttATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCtt                                                                             >  1:1186364/1‑149 (MQ=255)                                                                         tATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCCCAGCAGCCAGCGATATCGCCGTAACTTACGCTTCt                                                                           >  1:295909/1‑149 (MQ=255)                                                                          aTGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTa                                                                          <  1:1545678/149‑1 (MQ=255)                                                                           tGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTAt                                                                         <  2:1360386/149‑1 (MQ=255)                                                                            ggCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATg                                                                        <  1:538381/149‑1 (MQ=255)                                                                            ggCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATg                                                                        <  1:465983/149‑1 (MQ=255)                                                                                cGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtga                                                                    <  1:1273525/149‑1 (MQ=255)                                                                                   cAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgatga                                                                 <  1:391543/149‑1 (MQ=255)                                                                                    aGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAc                                                                <  2:1186364/149‑1 (MQ=255)                                                                                     gATCCTCTACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAcc                                                               <  1:594458/149‑1 (MQ=255)                                                                                      aTCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGAGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAccc                                                              >  1:911854/1‑149 (MQ=255)                                                                                        ccTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCaa                                                            <  2:157354/149‑1 (MQ=255)                                                                                              cTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTc                                                      <  2:1127814/149‑1 (MQ=255)                                                                                              cTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTc                                                      <  2:117418/149‑1 (MQ=255)                                                                                              cTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTc                                                      <  1:14344/149‑1 (MQ=255)                                                                                                   ccGAAGCGCATCCACGTTTCCGTGAAGGCTATGCAATAGGCACCCCGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGCCAGCGATATCTACGTAACTAACGCTTCTATGCTGATGAACCAACGATTCGCGTc                                                 >  2:886220/1‑149 (MQ=255)                                                                                                      aaGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGt                                              >  1:75456/1‑149 (MQ=255)                                                                                                        gcgcATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTgg                                            <  2:1526207/149‑1 (MQ=255)                                                                                                        gcgcATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTgg                                            <  2:364259/149‑1 (MQ=255)                                                                                                         cgcATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACCGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTggg                                           >  1:523295/1‑149 (MQ=255)                                                                                                           cATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGcc                                         <  2:938846/149‑1 (MQ=255)                                                                                                             tCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGa                                       >  2:607086/1‑149 (MQ=255)                                                                                                              ccACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGAc                                      <  2:930881/149‑1 (MQ=255)                                                                                                               cACGTTTCCGTGATGGATATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCACCGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTAGGCTGATGACCCAACGATTCGCGTCGGTGGGCCGAcc                                     <  1:301003/149‑1 (MQ=255)                                                                                                                aCGTTTCCTTGATGGCTATGCAATAGGCAGCGGGCGCTGGCTCGCGTTGCCGTTTGCTATCTACGTACTGACCCACAGCAGAGAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCt                                    <  2:1382838/149‑1 (MQ=255)                                                                                                                aCGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACAGACCCACAGCAGACAGCGATAACGACGTAACTTACGCTTCTATGCTGATGCCCCAACGATTCGCGTCGGAGCGCCGccca                                    >  1:330052/1‑148 (MQ=255)                                                                                                                  gTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTAc                                  <  2:951934/149‑1 (MQ=255)                                                                                                                        gTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACATCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTgg                            <  1:1141896/149‑1 (MQ=255)                                                                                                                         tGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTAACGCTTCTATTCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTggg                           <  1:441255/149‑1 (MQ=255)                                                                                                                               cTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCCCCGCTTCGCGTCGGGGGGCCGACCTACCATTGGGTAcgcg                     >  2:1550944/1‑149 (MQ=255)                                                                                                                                  tGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGaaa                  <  1:1459385/149‑1 (MQ=255)                                                                                                                                    cAATAGGCACCGGGCGCTGGTGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACGGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGt                >  2:1165059/1‑149 (MQ=255)                                                                                                                                     aaacggCAACGCGCGCCAGCGCCCGGTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACATCGATATCGAAGTAACTAACGCGTCTATGCGGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATAGGGTACGCGAAAGTa               >  2:86232/5‑149 (MQ=255)                                                                                                                                     aaacggCAACGCGCGCCAGCGCCCGGTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTa               >  2:596252/5‑149 (MQ=255)                                                                                                                                     aaTAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCCACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTa               >  1:1454955/1‑149 (MQ=255)                                                                                                                                     aaTAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTa               >  1:1490241/1‑149 (MQ=255)                                                                                                                                      aacggCAACGCGCGCCAGCGCCCGGTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGCTTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTAt              >  1:296142/4‑149 (MQ=255)                                                                                                                                          gCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctg          >  2:382906/1‑149 (MQ=255)                                                                                                                                           cACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGCCCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctgg         >  1:1059073/1‑149 (MQ=255)                                                                                                                                            aCCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggc        <  1:1232540/149‑1 (MQ=255)                                                                                                                                            aCCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggc        >  2:900627/1‑149 (MQ=255)                                                                                                                                                  cgcTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGa  >  1:268268/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGA  >  NC_000913/4009622‑4009915

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