Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,698,874	0	A	G	52.9%	8.2 / 103.9	68	K221E (AAA→GAA)	bcsG	DUF3260 family cellulose production inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (16/16);  new base G (16/20);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCAAGAACT  >  NC_000913/3698729‑3699024                                                                                                                                                  |                                                                                                                                                       aacGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCTGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCg                                                                                                                                                   <  1:280398/151‑1 (MQ=255)  acgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGa                                                                                                                                                  <  2:462648/151‑1 (MQ=255)   cgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGAc                                                                                                                                                 >  1:624194/1‑151 (MQ=255)     acgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCt                                                                                                                                               >  2:370763/1‑151 (MQ=255)          tggtAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCg                                                                                                                                          >  2:234862/1‑151 (MQ=255)              aaCGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCtt                                                                                                                                      >  1:441444/1‑151 (MQ=255)              aaCGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCtt                                                                                                                                      >  1:363751/1‑151 (MQ=255)              aaCGCAGCGGCAACCGTTGCGGAGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCt                                                                                                                                       <  2:516836/150‑1 (MQ=255)                cGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCg                                                                                                                                    <  2:238521/151‑1 (MQ=255)                cGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCg                                                                                                                                    <  2:267034/151‑1 (MQ=255)                  cAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCt                                                                                                                                  <  2:279471/151‑1 (MQ=255)                   aGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTg                                                                                                                                 >  2:352132/1‑151 (MQ=255)                    gCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGc                                                                                                                                >  2:52780/1‑151 (MQ=255)                     cGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGcc                                                                                                                               >  2:313205/1‑151 (MQ=255)                       gCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCg                                                                                                                             <  2:616318/151‑1 (MQ=255)                                 cggcgACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCa                                                                                                                   <  2:570027/151‑1 (MQ=255)                                   gcgACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGc                                                                                                                 >  2:620619/1‑151 (MQ=255)                                            ggCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGaa                                                                                                         <  1:156434/150‑1 (MQ=255)                                                cACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAactact                                                                                                    >  2:157194/1‑151 (MQ=255)                                                  ccGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTg                                                                                                   <  1:142589/150‑1 (MQ=255)                                                   cGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGt                                                                                                 >  1:508934/1‑151 (MQ=255)                                                    ggTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTg                                                                                                <  1:581951/151‑1 (MQ=255)                                                       aGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGAt                                                                                              <  2:439207/150‑1 (MQ=255)                                                        gTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGAtt                                                                                             >  1:220755/1‑150 (MQ=255)                                                          gggTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAAc                                                                                          <  1:141766/151‑1 (MQ=255)                                                              gATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATc                                                                                       <  2:235198/150‑1 (MQ=255)                                                               atatGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCt                                                                                      <  2:128073/150‑1 (MQ=255)                                                               atatGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTg                                                                                     <  1:437777/151‑1 (MQ=255)                                                                 atGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGtt                                                                                   <  1:320646/151‑1 (MQ=255)                                                                   gCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCg                                                                                 <  2:483141/151‑1 (MQ=255)                                                                    cccGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGc                                                                                <  2:153104/151‑1 (MQ=255)                                                                       gCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCttt                                                                             >  2:595639/1‑151 (MQ=255)                                                                        cacaAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTc                                                                            >  1:582966/1‑151 (MQ=255)                                                                             aCTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGt                                                                       >  2:417267/1‑151 (MQ=255)                                                                               tGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCg                                                                     >  1:408146/1‑151 (MQ=255)                                                                               tGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCg                                                                     >  2:123704/1‑151 (MQ=255)                                                                                gCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTc                                                                      <  2:595743/149‑1 (MQ=255)                                                                                gCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCgg                                                                    >  2:617903/1‑151 (MQ=255)                                                                                 cACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCgg                                                                    >  2:628250/1‑150 (MQ=255)                                                                                    cGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGatat                                                                <  2:241439/151‑1 (MQ=255)                                                                                     gCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGAtata                                                               <  1:15308/151‑1 (MQ=255)                                                                                        aacaacGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGaa                                                            <  2:191822/151‑1 (MQ=255)                                                                                        aacaacGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGaa                                                            <  2:627129/151‑1 (MQ=255)                                                                                        aacaacGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGa                                                             >  1:54014/1‑150 (MQ=255)                                                                                         acaacGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAg                                                           <  2:37325/151‑1 (MQ=255)                                                                                           aacGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAacc                                                         >  2:192299/1‑151 (MQ=255)                                                                                            acGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAgccg                                                        <  1:199466/151‑1 (MQ=255)                                                                                              ggCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAgccgcc                                                      <  1:462610/151‑1 (MQ=255)                                                                                                cGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCgg                                                    >  1:593946/1‑151 (MQ=255)                                                                                                cGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCgg                                                    >  2:459386/1‑151 (MQ=255)                                                                                                  aaCCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGt                                                  <  1:154448/151‑1 (MQ=255)                                                                                                  aaCCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGt                                                  >  2:36255/1‑151 (MQ=255)                                                                                                          cGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCg                                          >  1:468863/1‑151 (MQ=255)                                                                                                           gCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGc                                         <  2:195153/151‑1 (MQ=255)                                                                                                            ccTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCa                                        <  2:614060/151‑1 (MQ=255)                                                                                                                   gAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACt                                  <  2:565435/150‑1 (MQ=255)                                                                                                                   gAATAATTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGACGCCGGGTTGATGTCGCATCCACtg                                 >  1:304393/1‑151 (MQ=255)                                                                                                                   gAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCAc                                   <  2:421086/149‑1 (MQ=255)                                                                                                                        atTTCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTc                            <  2:626604/151‑1 (MQ=255)                                                                                                                           tCTATAACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGc                          >  2:503858/1‑150 (MQ=255)                                                                                                                              ataACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCa                         <  2:31819/148‑1 (MQ=255)                                                                                                                              ataACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCAt                        <  1:137187/149‑1 (MQ=255)                                                                                                                               taACGCGGAGGCGAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGAGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTc                     >  2:19438/1‑151 (MQ=255)                                                                                                                                           gAAACGTGAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTc         >  2:205450/1‑151 (MQ=255)                                                                                                                                            aaaCGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTc         >  1:184348/1‑150 (MQ=255)                                                                                                                                              aCGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCAAg      <  1:417274/151‑1 (MQ=255)                                                                                                                                               cGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCa        >  2:296174/1‑148 (MQ=255)                                                                                                                                                gTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCAAGa     >  2:185863/1‑150 (MQ=255)                                                                                                                                                  aaaTCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCAAGAAc   <  1:61963/150‑1 (MQ=255)                                                                                                                                                  aaaTCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCAAGAACt  >  1:297884/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                       AACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTCGCTTTCCTGGTCGGATATAGAAGCCGCCGGGTTGATGTCGCATCCACTGTGGTCGCATTTCGATATTGAGTTCAAGAACT  >  NC_000913/3698729‑3699024

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