Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,173	0	T	C	40.0%	105.7 / 19.3	65	I310V (ATT→GTT)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (31/8);  new base C (0/26);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.83e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTA  >  NC_000913/2921023‑2921317                                                                                                                                                       |                                                                                                                                                 agGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAaataat                                                                                                                                                  <  1:424855/151‑1 (MQ=255)  ggCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATc                                                                                                                                                 >  1:348663/1‑151 (MQ=255)     cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATccca                                                                                                                                              <  1:88515/151‑1 (MQ=255)     cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATccc                                                                                                                                               >  2:807632/1‑150 (MQ=255)     cAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATccc                                                                                                                                               >  2:125428/1‑150 (MQ=255)        cgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                           >  2:331052/1‑151 (MQ=255)         gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                           >  1:93583/1‑150 (MQ=255)         gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCACAACTCCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccca                                                                                                                                          <  1:755686/151‑1 (MQ=255)            aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                                                                                        >  1:255872/1‑150 (MQ=255)              aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                      <  1:265678/150‑1 (MQ=255)                    gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgccg                                                                                                                               <  1:766451/151‑1 (MQ=255)                       gatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGAt                                                                                                                             >  2:42971/1‑150 (MQ=255)                          catcCACTCAACGTTGACGTAGTTGCAGAACCCCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCCCCCAGCACACCGCCGATaaaa                                                                                                                         <  2:885476/151‑1 (MQ=255)                              cACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                     >  2:500469/1‑151 (MQ=255)                                 tCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaa                                                                                                                   >  1:502525/1‑150 (MQ=255)                                    aCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                               >  1:223887/1‑151 (MQ=255)                                         gACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGAGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcgg                                                                                                           >  1:481800/1‑150 (MQ=255)                                                ttGCAGAACCCCACCCCCATCGACAGCAACATGCCCATTACGATCGGTGTTTTCCGCGAAATGTTCAGCGACCCCGTGGGGCGCATCAGCCAGCCGAAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                       <  2:743329/147‑1 (MQ=255)                                                    aGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                               >  1:567965/1‑151 (MQ=255)                                                      acccccttcacctccGACAGCAACAGCCCCATTCCGATCGGTTTTTACCGCGAAATGTCCAGCAATCCGGGGCGGCGCATCAGCCAGTCGGAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaa                                                                                             <  1:3311/144‑1 (MQ=255)                                                      acccccacccccaccGACAGCAACAGCCCCATACCGAGCGGTGTTTTACGCGAAATGTCCAGCGACCCCGTGGGCCGCACCGCCCAGCCGAAAAAAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGc                                                                                                  <  1:117218/136‑1 (MQ=255)                                                      aaCACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcaca                                                                                              >  1:872302/1‑150 (MQ=255)                                                        caccatcaccatCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaa                                                                                             >  2:165397/1‑149 (MQ=255)                                                          ccacccccaacaaCAGCAACACGCCCGTTACGATGGGGGTTTTACGCGCAGTGTGGAGCAACCCCGGGCGGCGCTTAACCCTGTCAAAAATACCCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                                                                                          <  1:279513/138‑1 (MQ=255)                                                           catccccaccGACGGCAACAGCCCCATTCCGATCGGTTTTTTACGCGCAATTTTCAGCGACCCCGTGCGGCGCATCAGCCAGCCGGAAAAAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                        <  1:273432/151‑1 (MQ=255)                                                           catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                        >  1:38027/1‑151 (MQ=255)                                                                ccatcGACACCAACATGCCCTTTACGATCGGTTTTTTCCGCGAAATGTTCAGCGACCCCGTGCGGCGCATCACCCAGTCGAAAAAAATCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCtt                                                                                    <  2:1077395/150‑1 (MQ=255)                                                                  atcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                  >  2:639867/1‑150 (MQ=255)                                                                     gACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                               >  1:919366/1‑150 (MQ=255)                                                                      aCAGCAACATGCCAATTACGATCGGTTTTTTACGCGCAATGTCCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATACCCCCACCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                               <  1:836341/149‑1 (MQ=255)                                                                      aCAGAAACATGCCCATTACGATCGGTTTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa                                                                                 <  2:580760/147‑1 (MQ=255)                                                                           aaCATGCCCATTACGATGGGTGTTTTACGCGCAATGTTCAGCGACCCCGTGCGGCGCATCAGCCAGTGGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                            <  2:880085/147‑1 (MQ=255)                                                                            aCATCCCCATTACGAGCGGTTTTTTCCGCGAATGGTCCAGCGCCCCCGGGCGGCGCACCAGCCATGCGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGa                                                                          <  1:541800/148‑1 (MQ=255)                                                                                 cccATTACGATCGGTTTTTTCCGCGCAATGTCCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATcc                                                                     <  2:608479/148‑1 (MQ=255)                                                                                 cccATTACGACCGGTTTTTTCCGCCCAATGTCCAGCGACCCCGGGCGGCGCACCACCCAGCCGAAAAAACCCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                  <  1:366693/151‑1 (MQ=255)                                                                                  ccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                 >  2:470714/1‑151 (MQ=255)                                                                                       aCGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                            >  2:1026838/1‑151 (MQ=255)                                                                                             ggTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                       >  1:540647/1‑150 (MQ=255)                                                                                                 ttttcggcgcATTGTCCAGCGACCCCGTGGGGCCCATCACCCAGCCGGAAAAAACCCCCCCCAGAACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGAt                                                   <  1:573438/144‑1 (MQ=255)                                                                                                    tACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                               >  1:894337/1‑151 (MQ=255)                                                                                                      cgcgAAATGTCCAGCGATCCCGGGCGGCGCAACAGCCAGTCGAAAATAACCCCACCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                               <  1:145277/149‑1 (MQ=255)                                                                                                          cAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                          <  2:223891/150‑1 (MQ=255)                                                                                                           aaTGTCCAGCGCCCCCGGGCGGCGCACCGGCCGGTCGGAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                              <  1:807461/145‑1 (MQ=255)                                                                                                               tccaccgcccccGTGCGCCGCAACAGCCAGTCGGAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                               <  1:332596/131‑1 (MQ=255)                                                                                                                    cGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                >  2:1064218/1‑150 (MQ=255)                                                                                                                    cGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                >  2:873342/1‑150 (MQ=255)                                                                                                                    cGACCCCGTGCGGCGCATCACCCAGTCGAAAATAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaa                                 <  2:173319/149‑1 (MQ=255)                                                                                                                    cGACCCCGTGCGGCCCATCACCCAGCCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                <  1:846572/150‑1 (MQ=255)                                                                                                                      aTCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                              >  1:8063/1‑150 (MQ=255)                                                                                                                           ggggggcgCACCACCCAGCGGAAAATAACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaag                         <  1:597550/146‑1 (MQ=255)                                                                                                                             gcggcgCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                       >  1:730794/1‑150 (MQ=255)                                                                                                                              cggcgCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                       >  1:1025123/1‑149 (MQ=255)                                                                                                                               ggcgCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                    >  2:856987/1‑151 (MQ=255)                                                                                                                                gcgcATCAGCCAGTCGGAAATAATCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                        <  2:919388/146‑1 (MQ=255)                                                                                                                                   cATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAgg                >  1:119052/1‑151 (MQ=255)                                                                                                                                      cAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg             >  2:92088/1‑151 (MQ=255)                                                                                                                                         ccAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                   <  2:354083/142‑1 (MQ=255)                                                                                                                                           ggccgaaaaTAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTg          <  1:499893/143‑1 (MQ=255)                                                                                                                                            gTCGGAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATa       <  2:172361/151‑1 (MQ=255)                                                                                                                                            gTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt        <  1:495411/150‑1 (MQ=255)                                                                                                                                              cGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc      >  2:1113621/1‑150 (MQ=255)                                                                                                                                               gaaaaTAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa     <  2:1038412/148‑1 (MQ=255)                                                                                                                                                 aaaTAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt   >  1:283259/1‑150 (MQ=255)                                                                                                                                                 aaaTAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt   <  2:960280/150‑1 (MQ=255)                                                                                                                                                 aaaTAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt   <  2:508338/150‑1 (MQ=255)                                                                                                                                                 aaaTAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTa  <  2:817870/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 AGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTA  >  NC_000913/2921023‑2921317

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