Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	42.4%	145.3 / 20.4	86	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (38/11);  new base C (0/36);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 2.06e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTG  >  NC_000913/907453‑907746                                                                                                                                                         |                                                                                                                                                 cGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑Ac                                                                                                                                                 >  4:377067/1‑151 (MQ=255)  gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑Acc                                                                                                                                                >  2:181166/1‑151 (MQ=255)  gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑Acc                                                                                                                                                <  2:139598/151‑1 (MQ=255)    aGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCg                                                                                                                                               <  4:1161655/150‑1 (MQ=255)     gcgcGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCt                                                                                                                                             <  1:862008/151‑1 (MQ=255)      cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTc                                                                                                                                            <  1:159620/151‑1 (MQ=255)      cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTc                                                                                                                                            <  1:773119/151‑1 (MQ=255)       gcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTc                                                                                                                                            >  4:665197/1‑150 (MQ=255)        cgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCAc                                                                                                                                          >  1:1374345/1‑151 (MQ=255)         ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCAcc                                                                                                                                         <  1:203680/151‑1 (MQ=255)         ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGAGGGCGGTACGGAAGCAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGGGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGTTCAcc                                                                                                                                         >  4:983145/1‑151 (MQ=255)           tATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGc                                                                                                                                       >  3:24838/1‑151 (MQ=255)            atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGcc                                                                                                                                      <  3:31647/151‑1 (MQ=255)             tcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGcc                                                                                                                                      >  4:749184/1‑150 (MQ=255)                  cGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATaaa                                                                                                                                 >  1:530623/1‑150 (MQ=255)                   gCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGCGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGCT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAAt                                                                                                                               >  1:114737/1‑151 (MQ=255)                      agcagATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGtt                                                                                                                            >  2:178709/1‑151 (MQ=255)                      agaagttggtttttCACCCCCTGAATCACCCCCCTGGGGGGGGGGGGCGGTACGGAAGTAATCACATAGAAAACCCCCCCAGAAACCCAGCGGGTGGGGAGGGGGTG‑GGGGATCAACCGGCGG‑CGCTT‑CCCCCTCACCGCCATAAAATGtt                                                                                                                            <  4:842737/140‑1 (MQ=255)                       gcagATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGGTc                                                                                                                           >  3:186458/1‑151 (MQ=255)                          gatgtttttccAGCCCCTAAACCACCGCCGTGGGGGGGGGGGGCGGTAGGAAAATAATCAAATTGAGACCCCACCCCAAACCCCGGCGGGGGGGGGGGGGGGG‑GGGGATCAACCGGGGC‑CTTTT‑CCCGCTCACCGCCATAAAATGtt                                                                                                                            <  4:539198/137‑1 (MQ=255)                                 ttccacccccTGAACCACCCCCGGGGGGGGGGGGGGGGGTAGGGAAGTAATCACATGGAAAACCCACCCCAAAACCCAGGGGGGGGGGGGGGGGGG‑GGGGATCAACCGGCGT‑CGCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                  <  4:795351/144‑1 (MQ=255)                                    cAGCCCGTAAAACAACCCCCGGGCGGTGGTGGGCCTCACGGAATCACTCACATTAAAAACCCACACCAAAGCCCCGCGGGGGGGGAGTGGGTG‑GGCAAACACCCAGCGG‑CTCTT‑CCGGCTAGCCGCCATAAAATGTTCAAACTGtt                                                                                                                   <  4:286868/146‑1 (MQ=255)                                       cccttgaacaaccGCCCGGGGGGGGGGGGGGGGTAGGGAAGTAACCACATTGAAAACCCCCCCCAAACCCCGGGGGGGGGGGGGGGGGTT‑GGGAATCACCCGGGGG‑TGTTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                                            <  4:301978/138‑1 (MQ=255)                                       cccctaatccaccccccTGGCGGTGGTGGGCGGTACGAAAGAAATCACATTGAGAACCCACCCCAAAACCCAGCGGGTGGGAGGGGGTAT‑GGGTACCACCCTCCGC‑TGCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATaa                                                                                                           <  1:474652/137‑1 (MQ=255)                                       ccccTGAATCACCGCCCGGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAAAACCCACACCAAAACCCGGGGGGGGGGGGGGGGGAG‑GGGGATCACCCGGCGG‑CTTTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                                             <  1:1008706/149‑1 (MQ=255)                                       ccccTAAATCACCGCCCGGGGGGGGGGGGGCGGGGCGGAAAGAATCAGATTGGAACCCCACGCCAAAACCCGGGGGGGGGGGGGGGGGAT‑GGGGAAAACCCAGGGG‑GCCTT‑ACCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATaa                                                                                                           <  1:705773/151‑1 (MQ=255)                                       ccccTAAATCACCCCCGGGGGGGTGGGGGGGGGTAGGGAATTAACCACATTGAAAACCCACACCAAAACCCAGGGGGTGGGGGGGGGGAG‑GGGGATCAACCAGCGG‑CGCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATaa                                                                                                           <  1:118146/151‑1 (MQ=255)                                       cccaTGAATCACCGCCGGGGGGGGGGGGGGGGGTACGGAAGAAATCAAATTGAAAACCCCCCCCGAAAGCCGGGGGGGGGGGGGGGGGGG‑GGGGACCACCGGCCGC‑CGTTC‑CCCGCTCGCCGCCATAAAATGTTCAAACTGTTGCTGa                                                                                                              <  4:38117/148‑1 (MQ=255)                                           tGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                        >  1:587619/1‑150 (MQ=255)                                              aTCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                     >  2:17637/1‑150 (MQ=255)                                                caccGCCCGGGGGGTGGGGGGGGGTAGGGAAGTAATCACATTAAGAACCCCCCCCGAAACCCGGGGGGGGGGGGGGGGGAG‑GGGAACCACCCAGCGG‑CCCTT‑CCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAc                                                                                                   <  3:171490/147‑1 (MQ=255)                                                cAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCCGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                  >  1:932850/1‑151 (MQ=255)                                                    gCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTg                                                                                              >  1:237606/1‑151 (MQ=255)                                                        ggggggggggggggggTAGGAAAGAAATCACATGAAAAACCCCCCCCAAAACCCGGGGGGGGGGGAGGGGGTG‑GGGGATCAACCGGCGG‑CGCTT‑CCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTg                                                                                              <  4:603231/137‑1 (MQ=255)                                                         gggggtgggGGGGGTTACGGAAGAAACCACATTGAAAACCCACCCCAAACCCCGGGGGGGGGGGGGGGGGAG‑GGGGATAACCCAGGGG‑CGCTG‑CCCGCTAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCa                                                                                                    <  4:1405635/137‑1 (MQ=255)                                                         gggggggggggggggTAGGGAATTAATCACATTGAAAACCCCCCCCAAAACCCGGGGGGGGGGGAGGGGGAG‑GGGAATCAACCAGCGG‑CGCTT‑CCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                          <  3:136956/141‑1 (MQ=255)                                                         ggCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                          >  3:174848/1‑150 (MQ=255)                                                            ggtggtGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                      >  1:733627/1‑151 (MQ=255)                                                            ggtggtGGGCGGTACGGAAGTAATCACATTAAAAAACCACCCCAAAACCCAGGGGGGGGGGGGGGGGAG‑GGGGATCAACCAGCGG‑CGCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                       <  1:807198/150‑1 (MQ=255)                                                              tggtggGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCggcggc                                                                                    >  3:172333/1‑151 (MQ=255)                                                               ggggggCGGTACGAAATTAATCACATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGA‑‑GGGGCACCACCAGGGGCCGCTC‑CCCCCTCGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCt                                                                                   <  4:1342327/148‑1 (MQ=255)                                                                gggggCGGTACGAAAGAAACCAAATTAAGAACCCCCCAGAAAACCCGGGGGGGGGGGAGGGGGTG‑GGGAATCACCCGGCGG‑CCCTT‑CCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCtt                                                                                  <  4:1266819/149‑1 (MQ=255)                                                                 tggGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                 >  4:1236394/1‑151 (MQ=255)                                                                      gggcgggAGGTATTCGCATTGAAACACCCCACCGAAACCCGGGGGGTGGGGAGGGGGTT‑GGGAATCAACCAGGGG‑CGCTT‑CCCCCCCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGaa                                                                              <  4:63524/144‑1 (MQ=255)                                                                          cGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                        >  4:676163/1‑151 (MQ=255)                                                                             aaGTAATCACATGGAGAAACCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAG‑GGGGATCAACCAGCGG‑CGTTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgt                                                                      <  1:1287673/150‑1 (MQ=255)                                                                              aGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtgt                                                                    >  4:439181/1‑151 (MQ=255)                                                                                       aTTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                           >  1:1473447/1‑151 (MQ=255)                                                                                        ttGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                          >  4:33222/1‑151 (MQ=255)                                                                                        tagaaaacccaaccagAATTCCGGGGGGGGGGGGGGGGGGG‑GGGAATC‑ACCCGCCG‑GGCTTCCCCTCTCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                           <  3:148322/135‑1 (MQ=255)                                                                                          gagaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                         >  4:749688/1‑150 (MQ=255)                                                                                           agaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                       >  2:119736/1‑151 (MQ=255)                                                                                           agaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                       >  2:65763/1‑151 (MQ=255)                                                                                           aaaacccACCCCAAAACCCGGGGGGGGGGGGGGGGGTG‑GGGGACCACCCGGGGC‑CTCTT‑CCCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                          <  4:426458/143‑1 (MQ=255)                                                                                           aaaaaCCACACAAAAACCCGGGGGGGGGGGGGGGGGAG‑GGGGACCACCGGGCGG‑CGCTT‑CCCGCACACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGga                                                        <  1:707027/148‑1 (MQ=255)                                                                                               cccACCCCGAAACCCGGGGGGTGGGGGGGGGGGG‑GGGGATAAACCGGGGG‑CGCTT‑CCCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                    <  1:623639/149‑1 (MQ=255)                                                                                                   cacCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaac                                               >  4:518595/1‑151 (MQ=255)                                                                                                     ccGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                              >  1:1301273/1‑150 (MQ=255)                                                                                                        aaacccggggggggggggggggggg‑ggGAATCACCCAGCGG‑CGTTT‑CCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                           <  4:1075642/141‑1 (MQ=255)                                                                                                        aaacccgggggggggggggggggAG‑GGGAACCACCCGGGGG‑CGCTA‑CCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                          <  4:449243/142‑1 (MQ=255)                                                                                                           tCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                       >  1:1342317/1‑151 (MQ=255)                                                                                                            ccAGGGGGGGGGGGGGGGGAG‑GGGAATCACCCAGCGG‑CGTTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                       <  1:882092/150‑1 (MQ=255)                                                                                                               gggggggggggggggggg‑ggggATCACCCGGCGG‑CGCTC‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCaa                                        <  4:637644/144‑1 (MQ=255)                                                                                                               ggggggggggggggggAG‑GGGGATCAACCGGCGG‑CGCTT‑CCCGCCAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCa                                     <  1:1445475/147‑1 (MQ=255)                                                                                                                 ggggggggggggggAG‑GGGGATCAACCGGCGG‑CGCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                 <  4:281738/151‑1 (MQ=255)                                                                                                                 ggggggggAGGGGGAG‑GGGGACAAACCGGCGC‑CGTTT‑ACCGCCAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                  <  1:665237/150‑1 (MQ=255)                                                                                                                   gTGGGGGGGGGGTGGGGGGATCAACAGGCGG‑CTGTC‑CCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                      <  1:298166/145‑1 (MQ=255)                                                                                                                    tGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTg                              >  4:764316/1‑151 (MQ=255)                                                                                                                     ggggAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgc                             >  1:1206168/1‑151 (MQ=255)                                                                                                                            gggat‑ggggatCACCCAGCGG‑CCCTC‑ACCGCTATCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGt                      <  1:721660/151‑1 (MQ=255)                                                                                                                             ggag‑gggaatCACCCGGGGG‑CGCTT‑CCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                     <  1:951172/151‑1 (MQ=255)                                                                                                                              gat‑ggggatCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                    >  4:1451189/1‑151 (MQ=255)                                                                                                                                t‑ggggatCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTc                  >  4:1165438/1‑151 (MQ=255)                                                                                                                                 tggtattCAACCGGCCG‑CCCTT‑CCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCTGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt                   <  1:618714/144‑1 (MQ=255)                                                                                                                                  ggggatCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                 >  1:549664/1‑151 (MQ=255)                                                                                                                                       tCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc            >  1:553177/1‑151 (MQ=255)                                                                                                                                        cacccAGCGG‑CGCTT‑CCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc            <  1:309924/147‑1 (MQ=255)                                                                                                                                         acccgccgg‑ggCTT‑GCCCCCCGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa           <  4:382788/144‑1 (MQ=255)                                                                                                                                         acccAGCGG‑CGCTT‑CCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcg             <  4:1340455/146‑1 (MQ=255)                                                                                                                                         aaCCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa           >  4:813174/1‑150 (MQ=255)                                                                                                                                          acccggggccgcTT‑CCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcg             <  4:932784/139‑1 (MQ=255)                                                                                                                                              gcgg‑cgCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAAt         <  1:1301311/147‑1 (MQ=255)                                                                                                                                                gg‑cgCTT‑CCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg      <  1:714126/148‑1 (MQ=255)                                                                                                                                                gg‑cgCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCg    >  4:1195180/1‑150 (MQ=255)                                                                                                                                                gg‑ccCTT‑CCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATgg       <  3:17635/147‑1 (MQ=255)                                                                                                                                                 g‑cgcTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg  >  3:94471/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                 CGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGAT‑GGGGATCAACCAGCGG‑CGCTT‑ACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTG  >  NC_000913/907453‑907746

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