Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,941	0	A	C	50.0%	79.4 / 24.4	66	E392D (GAA→GAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (25/8);  new base C (0/33);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 2.65e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026797‑3027089                                                                                                                                                 |                                                                                                                                                     cAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCAcc                                                                                                                                                >  2:753163/1‑151 (MQ=255)  aaaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCg                                                                                                                                               >  2:679797/1‑151 (MQ=255)  aaaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCg                                                                                                                                               >  2:583401/1‑151 (MQ=255)    aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGcc                                                                                                                                             <  1:559121/151‑1 (MQ=255)      tATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCAc                                                                                                                                           >  2:997710/1‑151 (MQ=255)         tGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTcc                                                                                                                                        >  2:663878/1‑151 (MQ=255)           cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGc                                                                                                                                      >  2:934376/1‑151 (MQ=255)                      cTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTg                                                                                                                           >  1:1020306/1‑151 (MQ=255)                          cGCCTCTCGGCATATGAAGCTTTTTCCCTGGCCACGCTCGGGGGAGCGAAATCTCGGGGCTTTGACGTTTTGATGGGCAATTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCt                                                                                                                            <  2:251126/146‑1 (MQ=255)                          cGCCTCTCGGCATATGAAGCTTTTTACCTGGCCACGCTCGGCGGAGCAAAATCTCTGGGCTTTGACGTTTTGATTGGCAATTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                         <  2:623343/149‑1 (MQ=255)                          cGCCTCTCGGCATAGGAAGCGTTTTACCTGGCCACGCTGGGCGGAGCAAAATCTCGGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                       <  2:943773/151‑1 (MQ=255)                          cGCCTCCCGTCATATGAATCGTTTTCCCTGCCCCCGCCGGGCGGAGCGAATTCTCTGGCCCTTGACGATTTGATTGGCAGCTTTTTACCTGGAAAAGAGCCTGTTTTCGGGTTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                       <  2:1028159/151‑1 (MQ=255)                            ccTCTCGGCATATGAGGCTTTTTCCCTGCCCCCGCTCGGCGGGGCGAAATCTCGGGGCTTTGACTTTTTGATGGGCAATTTTTTCCCTGGCAAAGAGGCTGTTTTGGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTgc                                                                                                                          <  1:723863/146‑1 (MQ=255)                            ccTCTCGACATATGAAGTTTTTTCCCGGGCCCCCCTCGGGGGAGCAAAATCTCGGGGCCTTGACGTTTTGTTTGCCACTTTTTTCCCTGGAAAAGGGTCTGTTTCCGTGGTATTGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                       <  1:747477/149‑1 (MQ=255)                               cgcGGCATAGGAAGCTTTTTCCCGGGCCACGTCCGGCGGGGCAAATTCTCGGGGCTTTGACAATTTGTTAGCCAACTTTTTCCCGGACAAAGAGGCATATTTGGTGGTGATGGACCCCGCCGCCACTCCGCTACAGCAGCt                                                                                                                            <  1:238125/139‑1 (MQ=255)                                  ggtatatGAAGTGTTTTCCCGGGCCGCGCTCGGCGGAGCAAAATCTCTGGGCCTTGGCGATTTGATTGGAAATTTTTTACCGGGAAAAGGGGGTTATTGGGTGTTGATGGAACCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                                  <  2:270413/145‑1 (MQ=255)                                  ggCATATAAAGCGTTTTCCCGGGCCACGCTGGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCACCTTTTTCCCGGGAAAAGAGGCTGATTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa                                                                                                                 <  1:1094622/149‑1 (MQ=255)                                          aaGCTTTTTCCCTGGCCACGCTCGGCGGAGCAAAATCTCGGGGCTTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGttt                                                                                                        <  2:63790/150‑1 (MQ=255)                                          aaGCTTTTTCCCGGGCCCCGCTCGGCGGGGCAAAATCTCGGGGCCTTGACGTTTTGATTGCCAACTTTTTCCCGGGAAAAGAGGTTGATTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                         <  1:239891/149‑1 (MQ=255)                                              tttttCCCTGCCCCCGCTCGGCGGAGCAAAATCTCTGGCCTTTGACTTTTTGATTGCCACTTTTTTCCCTGGCAAAGGGGCTGTTTTCGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                      <  1:702890/147‑1 (MQ=255)                                               ttttCCCGGGCCACGCTGGGGGGAGCAAAATCTCGGGCCCTTGACGTTTTGATTGACATTTTTTTCCCTGGCAAAGAGGTTTTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTa                                                                                                   <  1:808486/150‑1 (MQ=255)                                               ttttACCTGCCCACGCTCGGCGGAGCGAAATCTCTGGGCTTTGACGTTTTGATTGGCATTTTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                    <  2:763983/149‑1 (MQ=255)                                                         ccgcgcgcgGCGGGGCGAAATCTGTGGGCTTTGACGTTTTGATTGACAACTTTTTACTTGCAAAAGAGGTTGATTGCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                                                                                            <  2:573518/144‑1 (MQ=255)                                                         ccccGCTCGGGGGAGCAAAATCTCGGGCCTTTGACGTTTTGATGGACACTTTTTTACCTGGAAAAGGGGCTTTTTTCGTGGTGATGGACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt                                                                                         <  1:40853/147‑1 (MQ=255)                                                          cACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCATTTTTTTACCGGGAAAAGAGGTTGATTTCGTGGTGAGGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTg                                                                                       <  2:226801/151‑1 (MQ=255)                                                               tCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGt                                                                                  >  2:509724/1‑151 (MQ=255)                                                                    ggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatg                                                                              >  2:57853/1‑150 (MQ=255)                                                                         gAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACg                                                                        >  1:221774/1‑151 (MQ=255)                                                                         gAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACg                                                                        >  2:49619/1‑151 (MQ=255)                                                                          aaaTCTCTGGGCCTTGACGTTTTGATTGGCACTTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                           <  2:521011/147‑1 (MQ=255)                                                                          aaaTCTCTGGGCCTTGACGATTTGTTTGGCAATTTTTTCCCGGGCAAAGAGGCTGATTTCGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACg                                                                        <  2:217264/150‑1 (MQ=255)                                                                          aaaTCTCGGGGCCTTGACGATTTGATGGGCAACTTTTTCCCGGGAAAAGAGGCTGATTTGGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGt                                                                       <  2:822810/151‑1 (MQ=255)                                                                                  gggCTTTGACGTTTTGTTTGGCAATTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGc                                                                  <  2:1070676/148‑1 (MQ=255)                                                                                  gggCCTTGACGATTTGTTTGGCAATTTTTTCCCTGGAAAAGGGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGa                                                                <  2:230452/150‑1 (MQ=255)                                                                                  gggCCTTGACGATTTGATTGGCAACTTTTTACCTGGAAAAGAGGCTGTTTTCGGGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                                    <  1:239826/146‑1 (MQ=255)                                                                                       ttGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCg                                                          >  2:8922/1‑151 (MQ=255)                                                                                       ttGACAATTTGATTGGCATCTTTTTCCCTGACAAAGAGGTTGTTTGCGGGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGa                                                             <  1:1066353/148‑1 (MQ=255)                                                                                        tGACGTTTTGATTGGCAATTTTTTACCGGGAAAAGGGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCg                                                          <  1:733091/150‑1 (MQ=255)                                                                                          aCGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTc                                                       >  1:1036239/1‑151 (MQ=255)                                                                                            gATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCg                                                      >  2:1123623/1‑150 (MQ=255)                                                                                            gATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCg                                                      >  2:963827/1‑150 (MQ=255)                                                                                            gATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCg                                                      >  2:727437/1‑150 (MQ=255)                                                                                             aTTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                    >  2:531086/1‑151 (MQ=255)                                                                                              ttggATTGGCATTTTTTCCCCTGGAAAAGAGGCTGTTTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                    <  2:556338/147‑1 (MQ=255)                                                                                                tGATTGGCATTTTTTTCCCTGGCAAAGAGGCTGATTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCt                                                  <  2:297326/150‑1 (MQ=255)                                                                                                       cAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAcc                                          >  2:639938/1‑151 (MQ=255)                                                                                                           tttttCCCTGGAAAAGAGGTTGTTTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTa                                        <  1:680443/149‑1 (MQ=255)                                                                                                                  cTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                               >  1:798423/1‑151 (MQ=255)                                                                                                                   tGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                               >  1:222948/1‑150 (MQ=255)                                                                                                                   tGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgt                              >  1:91420/1‑151 (MQ=255)                                                                                                                    ggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                             >  2:962044/1‑151 (MQ=255)                                                                                                                      cAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                           >  1:628122/1‑151 (MQ=255)                                                                                                                       aagggggCTGTTTGCGTGGGGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                             <  2:276842/143‑1 (MQ=255)                                                                                                                       aaaGAGGTTGTTTTGGTGGTGAGGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                          <  1:1002362/151‑1 (MQ=255)                                                                                                                       aaaGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                          >  2:509493/1‑151 (MQ=255)                                                                                                                        aaGAGGCTGTTTTGGTGGTGATGAAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCt                         <  2:4305/151‑1 (MQ=255)                                                                                                                        aaGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCt                         >  2:172019/1‑151 (MQ=255)                                                                                                                          gagGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTg                        <  2:776544/150‑1 (MQ=255)                                                                                                                            ggCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtg                     <  2:1020688/151‑1 (MQ=255)                                                                                                                                  tttCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGCCGTCTGGTGTACg                 <  1:106851/149‑1 (MQ=255)                                                                                                                                        ggtgATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCaa          <  2:451955/150‑1 (MQ=255)                                                                                                                                        ggtgAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGc            <  1:244171/148‑1 (MQ=255)                                                                                                                                         gtgATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt        <  2:384840/151‑1 (MQ=255)                                                                                                                                            agggACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTTTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaat     <  1:184609/149‑1 (MQ=255)                                                                                                                                            aaagAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACta       <  2:464479/146‑1 (MQ=255)                                                                                                                                              ggACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataa   <  2:1132164/151‑1 (MQ=255)                                                                                                                                                aaCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataat  <  2:249727/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAAT  >  NC_000913/3026797‑3027089

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