Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,089	0	A	C	42.4%	97.2 / 11.2	61	Y299D (TAT→GAT)	yagA	CP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (27/7);  new base C (0/25);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.14e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACG  >  NC_000913/280942‑281239                                                                                                                                                    |                                                                                                                                                        cccGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcg                                                                                                                                                      <  2:564600/151‑1 (MQ=255) cccGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcg                                                                                                                                                      >  2:523772/1‑151 (MQ=255)             tGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTAc                                                                                                                                         >  2:423581/1‑151 (MQ=255)               aaCGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTAcc                                                                                                                                        >  1:699024/1‑150 (MQ=255)                 cGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCATTCATCGTATTCCGGGGGCGTTGTTTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTAcc                                                                                                                                        <  2:1112468/148‑1 (MQ=255)                   ccTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCat                                                                                                                                   <  1:712386/151‑1 (MQ=255)                   ccTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCat                                                                                                                                   >  1:650532/1‑151 (MQ=255)                   ccTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCat                                                                                                                                   >  1:587388/1‑151 (MQ=255)                        ccTGCGCTCAGACTTCCCCCTTTCCCGCTCAGCTTTCCGCTGATATCCACTTCCCTGCCCATCCCCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCG‑GTACCGCCata                                                                                                                                  <  1:150300/147‑1 (MQ=255)                            cgcTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAgcgc                                                                                                                          >  2:558411/1‑151 (MQ=255)                             gcTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAg                                                                                                                             >  1:936265/1‑147 (MQ=255)                                catccTTCCCCCTTTCACGCTCAGCTTTCCGCGTATATCCACTTCCCTGACCATCCCCCCTCCATCGTATCCCGGGGGCGTTGTGTCGCCGCGGTACTGCCGCGCTGCCGGCTGATCCCGCGAGCCGG‑GTACCGCCATATCCAGCGCCTc                                                                                                                       <  2:904382/146‑1 (MQ=255)                                   aCTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGc                                                                                                                   >  1:129747/1‑151 (MQ=255)                                   aCTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGc                                                                                                                   >  2:533992/1‑151 (MQ=255)                                     ttcccccTTTCACGCCCAGCTTCCCGCGGATACCCCCTTTCCTGCCCCCCCCCCCTTCATCGTATCCGGGGGGCTTTGTGTTCCCGCTGTACTGCCGCTCGGCCGTCTGATCCCGCGAGCCGG‑GTACCGCCATATCCAGCGCCTCGTg                                                                                                                    <  2:964340/145‑1 (MQ=255)                                       cccccTTTCCCGCTCAGCTTTCCGCTGATATCCCCTTTCCTGACCATCCCCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCggg                                                                                                                <  1:948527/149‑1 (MQ=255)                                       aCCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGc                                                                                                               >  2:47087/1‑151 (MQ=255)                                        ccccTTTCACGCTCAGCTTTCCGCTGATACCCCCTTCCCTGACCATCCCCCCTTCATCGTATTCCGGGGGCTTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑TTACCGCCATATCCAGCGCCTCGTGCGGGc                                                                                                               <  1:281412/150‑1 (MQ=255)                                        ccccTTCCACGCTCAGTTTTCCGCTGATATCCACTTTCCTGACCATCCCCCCTTCATCGTATTCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCGGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGc                                                                                                               <  1:209227/150‑1 (MQ=255)                                         cccTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGt                                                                                                             >  1:762664/1‑151 (MQ=255)                                           cttacccccTCAGCTTCCCGCTGATATCCTCTTTCCTGACCACCCCCCCTTCATCGTCTCCGGGGGGCGTTGTGTTCCCGCTTTCCGGCCGCGCTGGCGGCTGATCCCGCGAGCCCG‑ATACCGCCATATCCAGCGCCTCGTGCGGGCGtt                                                                                                            <  2:1057865/142‑1 (MQ=255)                                           cTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTc                                                                                                           >  1:1039850/1‑151 (MQ=255)                                             ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                                                                                          >  2:466765/1‑150 (MQ=255)                                              tCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                                                                                          >  1:624223/1‑149 (MQ=255)                                              ccccccccAGCTTCCCGCTGATCTCCCTTTCCCCGACCACCCCCCCTTCATCGTATCCCGGGGGCGTTGTGTTGCCGCTGTTCGGCCGCGCTGCCGGCTGATCCCGCGAGCCAG‑TTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAg                                                                                                        <  2:213466/144‑1 (MQ=255)                                                 cccTCAGTTTCCCGCTGATGCCCCCTTCCCTCCCCATCCCCCCTTCACCGTTTCCGGGGGGCGTTTTGTTCCCGCTGTCCCGCCGCGCTGACGGCGGATCCCGCGACCCCG‑TTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGt                                                                                                      <  2:469474/148‑1 (MQ=255)                                                 cGCTCAGCTTCCCGCTGATCCCCCCTTCCCTGCCCACCCCCCCTTCATCTTATCCGGGGGGCTTTGTGTTGCCGCTGTCCTGCCGCGCTGACGGCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGtt                                                                                                     <  2:750184/151‑1 (MQ=255)                                                       gCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGAc                                                                                               >  2:123098/1‑151 (MQ=255)                                                         ttcccGCTGATATCCTTTTTCCTGACCATCCCCCCTCCGGCGTATCGGGGGGGCGTTGTGTGCCCCCGGTACCGCCGCGTTGCTGGGGGAT‑CCCCGGGCCCGTGTTCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAg                                                                                                 <  2:937699/144‑1 (MQ=255)                                                                 gATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGTGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGcc                                                                                       >  2:728919/1‑149 (MQ=255)                                                                  atatCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGt                                                                                    >  1:1077354/1‑151 (MQ=255)                                                                       cACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCg                                                                               >  1:532716/1‑151 (MQ=255)                                                                         cTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGa                                                                              >  2:273010/1‑150 (MQ=255)                                                                              cgcaccaccCCCCCTTCATGGTATTCGGGGGGTGTTGTTTTGCCGTTCAACGGCCGCGCTGACGGCTGACCCCGGCACCCCT‑GAACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAgg                                                                           <  1:647165/145‑1 (MQ=255)                                                                                  ccccccccccTCCTTCGTATCCCGGGGGCGTTGTGTTCCCGCTGTCCGCCCGCGCTGACGGCTGACCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcg                                                                       <  2:50414/142‑1 (MQ=255)                                                                                  ccacccccccTTCATCGTATCCCGGGGGCTTTGTGTTGCCGTTGTACTGCCGCGCTGACGGCTGATGCCGCGAGCCGG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCt                                                                     <  1:258853/144‑1 (MQ=255)                                                                                   cATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCAGCCAGCGGTCGAAGGCGCGCTg                                                                    >  1:573045/1‑150 (MQ=255)                                                                                        ccccTTCCCCGTATCCCGGGGGCTTTGTGTGGCCGCGGTACTGCCGCGCGGGCGGCTTATCCCGCGCGCCCG‑TTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTc                                                              <  2:91308/151‑1 (MQ=255)                                                                                        ccccTTCACCGTATCCCGGGGGCGTTGTGTTGCCGCTGTCCTGCCGCGCTGCGGGCTGATCCCGCGAGCCGG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTc                                                              <  2:563133/151‑1 (MQ=255)                                                                                               aTCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTg                                                       >  2:160353/1‑151 (MQ=255)                                                                                                tCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGt                                                      >  2:366098/1‑151 (MQ=255)                                                                                                 cGTATCCGGGGGGCGTTTTTTTGCCGCTTTACGGCCGCGCTGCCGCCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCt                                                        <  1:607909/148‑1 (MQ=255)                                                                                                   tttgcggggggTTTTGTTTTGCCGCTTTCCTCCCGCCCGGGGGGCTGACCCCGCGAGCCGG‑TTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCt                                                    <  2:906710/144‑1 (MQ=255)                                                                                                   tttccggggggCGTTGTTTTGCCGCTGTCCTCCCGCGCTGCCGGCTGTTCCCGCGAGCCGG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCg                                                          <  2:1001799/138‑1 (MQ=255)                                                                                                   tATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTg                                                   >  2:903323/1‑151 (MQ=255)                                                                                                   tATCCCGGGGGCTTTTTGTCCCCGCTGTCCTGCCGCCCTGCCGGCTGATACCGGGCCCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTg                                                       <  2:961833/147‑1 (MQ=255)                                                                                                   tATCCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGt                                                      <  1:369283/148‑1 (MQ=255)                                                                                                    atcccGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTg                                                   <  1:209542/147‑1 (MQ=255)                                                                                                       ccGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGa                                                >  1:12154/1‑150 (MQ=255)                                                                                                         gggggCGTTGTGTCCCCGGTGTCCTCCCGCGCTGCCGGCTGATACCGGGACCCCG‑GTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGaa                                               <  2:82909/149‑1 (MQ=255)                                                                                                           gggttttGTTTTCCCGCTGTCCCGCCGCGCGGCGGGATGATCCCGCGAGCCGG‑CTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCg                                                 <  2:493822/140‑1 (MQ=255)                                                                                                           gggCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCAt                                           >  1:309215/1‑151 (MQ=255)                                                                                                                    gtTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgc                                  >  2:603897/1‑151 (MQ=255)                                                                                                                     ttcccGCTGTACCGCCGCGCGGACGGCTGAGCCCGCGCGCCGG‑TTGCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgc                                  <  2:811067/147‑1 (MQ=255)                                                                                                                     ttGCCGCTGTCCGGCCGCGCTGCCGGCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCt                                    <  1:24737/148‑1 (MQ=255)                                                                                                                        ccgttttccTGCCGCGCGGGGGGCTGATCCCGCGAGCCGG‑TAACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgcagca                              <  2:819796/143‑1 (MQ=255)                                                                                                                        ccGCTTTCCTGCCGCGCTGCCGGCTGACCCCGCGACCCCG‑TTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgcag                                <  2:677291/149‑1 (MQ=255)                                                                                                                        ccGCTGTCCTGCCGCGCGGGGGGCTGATCCCGGGAGCCGG‑TTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgcagca                              <  2:95711/151‑1 (MQ=255)                                                                                                                              tcctcccGCGCTGCCGGCTGATCCCGGGGCCCGG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCAc                             <  2:55453/141‑1 (MQ=255)                                                                                                                                      cgcTGCCGGCTGATCCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAg                 <  1:782765/150‑1 (MQ=255)                                                                                                                                               cTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTg        >  1:707396/1‑150 (MQ=255)                                                                                                                                                    aCCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACg  >  2:889240/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCG‑GTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACG  >  NC_000913/280942‑281239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑