Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,009,776	G→C	6.4%	A203P (GCG→CCG)	pldB →	lysophospholipase L2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,776	0	G	C	6.4%	245.3 / 8.0	78	A203P (GCG→CCG)	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand):  ref base G (32/41);  new base C (3/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

AGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGT  >  NC_000913/4009628‑4009923                                                                                                                                                     |                                                                                                                                                    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:1159263/149‑1 (MQ=255)                cATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTa                                                                                                                                      >  2:1459385/1‑149 (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                                                                                                                                 <  2:1278013/149‑1 (MQ=255)                     cgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAc                                                                                                                                 >  1:364259/1‑149 (MQ=255)                     cgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAc                                                                                                                                 >  1:1526207/1‑149 (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)                                                                                                                                                  cgcgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgt    <  1:900627/149‑1 (MQ=255)                                                                                                                                                   gcgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgtg   >  2:1306032/1‑149 (MQ=255)                                                                                                                                                   gcgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgtg   >  2:1269680/1‑149 (MQ=255)                                                                                                                                                    cgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgtgt  <  2:1278696/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGT  >  NC_000913/4009628‑4009923

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