Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	289,251	0	T	C	44.1%	84.9 / 20.2	59	intergenic (‑89/+50)	yagI/argF	CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (28/5);  new base C (0/26);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.86e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCC  >  NC_000913/289105‑289400                                                                                                                                                     |                                                                                                                                                           gCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTtctcccgct‑ccc                                                                                                                                                        <  1:977802/150‑1 (MQ=255) gCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTctccttct‑ccc‑t                                                                                                                                                      <  2:409777/151‑1 (MQ=255)    cGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑t‑gtg                                                                                                                                                  >  2:501870/1‑151 (MQ=255)     gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTgg                                                                                                                                                 <  1:266354/151‑1 (MQ=255)       tGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTggg                                                                                                                                                >  2:37890/1‑150 (MQ=255)              gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAgg                                                                                                                                           <  1:146033/148‑1 (MQ=255)                cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑g                                                                                                                                      <  2:853495/150‑1 (MQ=255)                   ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGcc                                                                                                                                        >  2:980463/1‑146 (MQ=255)                   ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGcc                                                                                                                                        <  1:980245/146‑1 (MQ=255)                       cTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAgg                                                                                                                              >  1:965706/1‑151 (MQ=255)                         gCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑c                                                                                                                           >  1:95553/1‑151 (MQ=255)                            aCGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CAcc                                                                                                                        >  1:650558/1‑151 (MQ=255)                                    tCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACgcgc                                                                                                                 >  1:262088/1‑150 (MQ=255)                                     cAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACgcgc                                                                                                                 >  1:823626/1‑149 (MQ=255)                                      aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGc                                                                                                              >  1:585003/1‑151 (MQ=255)                                               aataatCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAAt                                                                                                     >  1:754540/1‑151 (MQ=255)                                                         aTAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCa                                                                                            >  1:924911/1‑150 (MQ=255)                                                          tAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCa                                                                                            >  2:796961/1‑149 (MQ=255)                                                              gCTGGCATGTCCACGCTGGAATTAGGGCGAAAA‑CGTCCCTTTTTCGGG‑ATGAAAAATCAAGGAAATGGAGCCGGTGTTCTCCCTCC‑CCC‑T‑GGGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCg                                                                                        <  1:71429/149‑1 (MQ=255)                                                                tGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGc                                                                                    >  1:649465/1‑151 (MQ=255)                                                                 ggCATGTCCACGCTGGAATTAAGGCAAAAA‑CTTCCCTTTTTCGGG‑ATAAAAACTCAATGAAATGGAGCCGGTGTTCCCCCTCC‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGc                                                                                    <  1:558255/150‑1 (MQ=38)                                                                      tccccccgttgAATTAAAGGAAAAA‑CCTTCCTTTTTCGGG‑GAAAAAAAAAAAAAAAAGGGGGGCGGTTTTCCCCCCCC‑CCCGG‑GGGGGGGGGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                                                                                  <  2:751815/138‑1 (MQ=255)                                                                      gTCCCCGCGGGAATTACGGCGAAAG‑CTTCCCTTTTTGGGG‑AGGAAAAATCAATGAAAGGGGGCGGTTTTTCTCCCCCT‑CCC‑G‑GAGGGAAAGGCCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                                                                                       <  1:715073/142‑1 (MQ=255)                                                                        ccccGCTGAATTAAAGGCAAAAC‑CGTCCTTTTTTCGGG‑ATAAAAAACAAATAAAAGGGGGCGGGTTTTCCCCCCCC‑CCT‑T‑GGGGGGGGGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                                                                               <  2:700286/145‑1 (MQ=255)                                                                        cccacggtgAAATAAGGGAAAAA‑CCTTCCTTTTTTGGG‑AAAAAAAAATAAAAAAAGGGGGGCGGTTTTCCCCCCCC‑CCC‑TGGGGGGAGGGGGCG‑GGGGGGGGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatc                                                                                <  2:942164/140‑1 (MQ=255)                                                                        ccACGC‑GGAAATAAGGGAAAAACCGTCCTTTTTTGGGA‑AAAAAAAAAAAATGAAAGGGGGCGGGTTTTCCCCCTCC‑CCT‑G‑GGGGGGGGGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCAtcatca                                                                            <  2:920924/151‑1 (MQ=255)                                                                         ccccgtGGAAATAAAGGGAAAA‑CCTTCCTTTTTGGGGAATAAAAAAAAAATAAAAGGGGCCGGGTTTTCCCCCCCC‑CCC‑T‑GGGGGGGGGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                                                                             <  2:775096/145‑1 (MQ=255)                                                                           cgcggaaaTTAAGGCGAAAA‑CGTCCCTTTTTCGGG‑ATAAAAAATCAATGAAAGGGAGCCGGTGTTCCCCCCCC‑CCC‑T‑GGGGGAGAGGGCC‑GGGGTGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc                                                                          <  1:1026620/144‑1 (MQ=255)                                                                            gCTGAAATAAGGGCAAAAC‑CGTCCTTTTTTGGGG‑ATAAAAAACAAATGAAGGGAGGCCGGTTTCCCCCCCCC‑CCC‑G‑GGGGGAGGGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCg                                                                         <  2:43215/150‑1 (MQ=255)                                                                                 aaTTAGGGCGAAAG‑CGTACCTTTTCCGGG‑TTGAAAAATCAATAAAGTGGAGCCGTTGTTCTCCCCCC‑CCC‑T‑GTGGGGGAGGGCG‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTa                                                                    <  1:718329/150‑1 (MQ=255)                                                                                  aTTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCtt                                                                     >  1:107627/1‑148 (MQ=255)                                                                                    tAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCg                                                                >  1:346492/1‑151 (MQ=255)                                                                                        ggtaaaactgtcccTTTTGCGGA‑ATAAAAAATCAATGAAGAGGAGCCGTTTTCCCCCCTCC‑CCC‑T‑GTGGGAGAGGCCC‑GGGGTGGGTG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                              <  1:553249/138‑1 (MQ=255)                                                                                         cGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGc                                                           >  1:695897/1‑151 (MQ=255)                                                                                          gAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGc                                                           >  2:537815/1‑150 (MQ=255)                                                                                            aaa‑CCTTCCTTTTTGGGG‑AAAAAAAAAAAAAGAAAGGGGGGCGGGTTTCCCCCCCC‑CCC‑TGGGGGGGGGGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATGACCGCCTTAATCGTATGc                                                           <  2:985498/149‑1 (MQ=255)                                                                                               accttccTTTTTGGGG‑AAAAAAAAATAAAGAAAGGGGGGCGGTTTTTCCCCCCC‑CCC‑T‑GGGGGGGGGGGCCGGGGGTGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATcc                                                       <  2:837082/146‑1 (MQ=255)                                                                                                    ccTTTTTCGGG‑AGGAAAAACGAATGAATGGGACCCGGTTTTCCCCCTCC‑CCC‑T‑TGGGGAGAGGCCC‑GGGGGGAGGC‑CCCCAACGCGCGGCGGACCCAACCACTCCCCAAGCGTTACCATCATCACCGCCTTAATCGTATGCATCCGGtttt                                                 <  1:666319/151‑1 (MQ=255)                                                                                                     cTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTc                                                >  2:968711/1‑151 (MQ=255)                                                                                                      tttttGGGG‑AAAAAAAAACAAATAAAAGGGGGCGGTTTTCCCCCCCC‑CCCGT‑TGGGGGGGGGGCG‑GGGGTGGGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGt                                                    <  2:605271/147‑1 (MQ=255)                                                                                                      tttttCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTcc                                               >  2:55372/1‑151 (MQ=255)                                                                                                        tttCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGc                                             >  2:204501/1‑151 (MQ=255)                                                                                                        ttccGGG‑ATGAAAAATAAAGGAAAGGGAGCCGGGTTCCCCCCCCC‑CCC‑G‑GTGGGGGAGGGCC‑GGGGTGGGGG‑CACCAACGCGCAGCGGACCCAATCGCTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTc                                                <  1:53364/145‑1 (MQ=255)                                                                                                            ggg‑ATAAAAAACAAGTGAAGGGGGGGCGGTTTCCCCCCTCC‑CCC‑T‑GGGGAGGAGGGCG‑GGGGTGAGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTg                                          <  1:468773/150‑1 (MQ=255)                                                                                                              g‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                       >  1:504109/1‑151 (MQ=255)                                                                                                              g‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                       >  1:851276/1‑151 (MQ=255)                                                                                                                aTGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCg                                      >  1:348414/1‑151 (MQ=255)                                                                                                                 aaaaaaaaaaaaGAAAGGGGGGCGGTTTTCCCCCCCC‑CCC‑TGGGGGGGGGGGGCG‑GGGGTGGGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                       <  2:944128/148‑1 (MQ=255)                                                                                                                   aaaagTCAATAAAGGGGAGCCGGTTTCCCCCCTCC‑CCC‑T‑GGGGGGGAGGGCG‑GGGGTGGGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCg                                      <  1:945702/148‑1 (MQ=255)                                                                                                                    aaaaTCAAGGAAAGGGGGCCGGTTTTCTCCCTCC‑CCC‑T‑GGGGGAGAGGGCG‑GGGGTGAGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAaca                                  <  1:73462/151‑1 (MQ=255)                                                                                                                     aaaTCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAacac                                 >  2:580171/1‑151 (MQ=255)                                                                                                                      aaTCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACg                                >  2:237005/1‑151 (MQ=255)                                                                                                                        tCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                              >  1:835690/1‑151 (MQ=255)                                                                                                                         aaaaaaaggggcgCCGGTTTTCCCCCCCC‑CCT‑T‑GGGGGAGAGGGCG‑GGGGGGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                              <  2:754707/138‑1 (MQ=255)                                                                                                                          aaggAAATGGAGCGGGTGTTCTCCCTCC‑CCC‑G‑GTGGGGGGGGGCC‑GGGGTGAGGC‑CCCCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATg                             <  1:585700/147‑1 (MQ=255)                                                                                                                          aaTGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGc                            >  2:816884/1‑151 (MQ=255)                                                                                                                              aaaGGGGGCGGGTGTCCCCCCCCC‑CCC‑G‑GGGGGGGGGGGCG‑GGGGGGGGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgc                        <  2:178563/151‑1 (MQ=255)                                                                                                                               aagggggCCGGTGTTCCCCCCCC‑CCC‑TGGGGGGAGGGGGGG‑GGGGGGGGGC‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgc                        <  2:403106/145‑1 (MQ=255)                                                                                                                                  gggaccGGTTTTCCCCCCCC‑CCC‑T‑GGGGGGGGGGGGC‑GGGGGGGGGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgccg                      <  2:690544/146‑1 (MQ=255)                                                                                                                                   gggCCGGTTTTCCCCCCCCCCCC‑T‑GGGGGGGGGGCCG‑GGGGTGGGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                       <  2:180521/146‑1 (MQ=255)                                                                                                                                   gAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGccgccg                   >  1:126259/1‑151 (MQ=255)                                                                                                                                                 cctcc‑ccc‑t‑TTGGGAGGGGCCC‑GGGGTGAGGC‑CCCCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAaca          <  1:370047/146‑1 (MQ=255)                                                                                                                                                   tct‑ccc‑t‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtc   >  2:961153/1‑151 (MQ=255)                                                                                                                                                    ct‑ccc‑t‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtcc  >  2:291672/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                           GCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAA‑CGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCT‑CCC‑T‑GTGGGAGAGGGCC‑GGGGTGAGGG‑CACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCC  >  NC_000913/289105‑289400

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