Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	50.0%	47.8 / 50.5	72	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (24/12);  new base A (0/36);  total (24/48)
Fisher's exact test for biased strand distribution p-value = 3.15e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTC  >  NC_000913/3177599‑3177898                                                                                                                                                      |                                                                                                                                                       acTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTc                                                                                                                                                        >  2:298330/1‑150 (MQ=255)  cTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTc                                                                                                                                                        <  2:151598/149‑1 (MQ=255)         ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTc                                                                                                                                               <  1:956528/151‑1 (MQ=255)         ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTc                                                                                                                                               <  2:756791/151‑1 (MQ=255)         ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGATAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTc                                                                                                                                               <  1:666334/151‑1 (MQ=255)             ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                           >  2:53176/1‑151 (MQ=255)              tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTccc                                                                                                                                             <  1:904664/148‑1 (MQ=255)              tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                          >  2:681817/1‑151 (MQ=255)               cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                          >  2:53098/1‑150 (MQ=255)                aaTCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCcgccgc                                                                                                                                        >  2:404344/1‑151 (MQ=255)                      cccGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTc                                                                                                                                   >  2:869530/1‑150 (MQ=255)                           cAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACt                                                                                                                             >  1:228315/1‑151 (MQ=255)                                  ctcCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGc                                                                                                                      >  2:741080/1‑151 (MQ=255)                                      aGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                                                                                  >  2:1042937/1‑151 (MQ=255)                                       gTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                                                                                  <  2:698736/150‑1 (MQ=255)                                        tAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGAAATCCGAATCTGCTAAATCAGCACAACTTCATCACGCCCGGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGtt                                                                                                                 <  2:755029/150‑1 (MQ=255)                                         aGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                                >  2:448773/1‑150 (MQ=255)                                              aGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                             >  1:315731/1‑148 (MQ=255)                                               gggtTTTCGCTGGTGTCGTCCGCGGCAACCCGAATCTGCTAAATCACCACACCTTCACCACGCACGGGGTCAAAGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGa                                                                                                          <  1:741057/150‑1 (MQ=255)                                               gggtTTTCGCGGGTGTCGTACGCGGAAACCCGAATCTGCTCAACCAGCACCCCTTCATCGCGCTCTGGGCCAAAGGGTAGCAAGACTGCGGGGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATa                                                                                                            <  2:912946/148‑1 (MQ=255)                                                gggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                        <  2:533160/151‑1 (MQ=255)                                                 ggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                        >  1:279997/1‑150 (MQ=255)                                                  gTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGAt                                                                                                      >  1:775713/1‑151 (MQ=255)                                                    ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAGGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                     <  2:119650/150‑1 (MQ=255)                                                     tCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGAc                                                                                                   >  1:823014/1‑151 (MQ=255)                                                       gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                  >  2:461096/1‑150 (MQ=255)                                                       gCTGGTGTCGTACGCGGCAATCCGAATCTGCAAAACCAGCCCACCTTCACCACGCCCTGGGTAAAAGGGAAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                                    <  2:755651/148‑1 (MQ=255)                                                         tGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                                                                                                >  2:181340/1‑150 (MQ=255)                                                          ggTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                                                                                                <  2:614306/149‑1 (MQ=255)                                                          ggTGTCGAACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTAAAAGGGTAGCAAGCCTGCGGCGGGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt                                                                                              <  1:623436/151‑1 (MQ=255)                                                          ggTCTCGTACGCGGCAATCCAAATCTGCTAAACCAGCACAACTTCACCCCGCCCTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGTTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                                                                                               <  2:552185/150‑1 (MQ=255)                                                             gtCGTACGCGGAAATCCGAATCTGCTCAACCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCCCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                            <  2:290725/150‑1 (MQ=255)                                                               cGTCCGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                         <  2:26505/151‑1 (MQ=255)                                                                      ggaaaTCCGAATCTGCTCAATCAGCACACCTTCACCCCGCACGGGGACAAAGGGTAGCAAGCCTGCGGCGTGCCCGCGCCAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt                                                                                    <  2:988067/146‑1 (MQ=255)                                                                         aaTCCGAATCTGCTAAATCAGCACACCTTCACCCCGCCCGGGGTAAAAGGGTAGCAAGACTGCGGCGGGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt                                                                                  <  2:769821/148‑1 (MQ=255)                                                                            ccGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                             >  1:630290/1‑150 (MQ=255)                                                                             cGAATCTGCTCAATCAGCACAACTTCATCCCGCCCTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaa                                                                           <  2:836796/151‑1 (MQ=255)                                                                             cGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaa                                                                            >  2:514052/1‑150 (MQ=255)                                                                               aaTCTGCTCAATCAGCCCAACTCCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                          <  1:298325/150‑1 (MQ=255)                                                                               aaTCTGCTCAACCAGCACACCTTCATCACGCACTGGTTCAAAGGATACAATGACGCCGGCGTGACCGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                         <  2:568725/151‑1 (MQ=255)                                                                               aaTCTCCTCAACCAGCACAACTTCATCACCCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                          <  1:555455/150‑1 (MQ=255)                                                                                     ctaaaTCAGCACACCTCCATCACGCACTGGGTAAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                        <  2:804170/143‑1 (MQ=255)                                                                                        aaTCAGCACAACTTCATCCCCCACGGGGCCAAGGGGTAGCAAAACTGCGGCGTGCCCGCCCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                        <  2:70350/143‑1 (MQ=255)                                                                                        aaTCAGCACAACTTCACCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                                 <  2:1020452/150‑1 (MQ=255)                                                                                         accaccacaACTCCATCACGCACTGGGTCAAAGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                                 <  1:1043799/144‑1 (MQ=255)                                                                                            aGCACAACTTCACCACGCACTGGGTAAAAGGGTACCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                               <  1:277818/148‑1 (MQ=255)                                                                                              cacaACTTCATCACGCACTGGGTAAAAGGGTAGCAAGACGGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAg                                                           <  1:669519/150‑1 (MQ=255)                                                                                               acaACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgt                                                          <  1:335042/150‑1 (MQ=255)                                                                                                 aaCTTCATCCCGCCCGGGTTAAAGGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtg                                                         <  2:936696/149‑1 (MQ=255)                                                                                                   cTTCACCCCGCGCTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                        <  2:689612/148‑1 (MQ=255)                                                                                                    ttCATCCCGCACTGGTTCAAAGGGTAACAAGACGGGGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                      <  1:315328/149‑1 (MQ=255)                                                                                                    tccctccccctcTGGGCAAAAGGTTAGCAAGTCGGCGGGGTGCCCGCCTAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtg                                                         <  1:231517/135‑1 (MQ=255)                                                                                                     tcatcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                   >  1:640216/1‑151 (MQ=255)                                                                                                      catcaCGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                   <  2:578126/150‑1 (MQ=255)                                                                                                       atcaCGCACGGGGTCAAAGGGTAGCAAGACTGGGGCGTGCCCGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                      <  2:678196/146‑1 (MQ=255)                                                                                                         cccGCACTGGGTAAAAGGTTAGCAGGGCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                                     <  2:629327/143‑1 (MQ=255)                                                                                                          aCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                               >  1:488957/1‑150 (MQ=255)                                                                                                           cGCACTGGGTCAAAGGGTAACAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGt                                                 <  1:962751/147‑1 (MQ=255)                                                                                                             cACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                            >  1:978091/1‑150 (MQ=255)                                                                                                               cTGGGCCAAAGGGTAGCAAGACTGCGGCGGGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                            <  2:700747/148‑1 (MQ=255)                                                                                                                 gggTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataat                                         <  2:715799/149‑1 (MQ=255)                                                                                                                 gggTAAAAGGGTAGCAAGCCGGCGGCGTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAtt                                        <  2:626629/150‑1 (MQ=255)                                                                                                                 gggTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataa                                          <  1:318231/148‑1 (MQ=255)                                                                                                                          ggTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGt                              <  2:489013/151‑1 (MQ=255)                                                                                                                             aGCAAGACTGCGGCTTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtt                             <  2:864374/149‑1 (MQ=255)                                                                                                                             aGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttt                            <  1:553836/150‑1 (MQ=255)                                                                                                                                aaGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttt                            <  1:448762/147‑1 (MQ=255)                                                                                                                                   aCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa                      >  1:572222/1‑150 (MQ=255)                                                                                                                                     tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAACCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaaa                    >  2:282152/1‑150 (MQ=255)                                                                                                                                        ggcgTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa                <  1:336502/151‑1 (MQ=255)                                                                                                                                           gtccccGCGCCCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg             <  2:299763/147‑1 (MQ=255)                                                                                                                                              aCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc          >  1:339725/1‑151 (MQ=255)                                                                                                                                                  cgcTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAgg       <  1:521667/150‑1 (MQ=255)                                                                                                                                                    cTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTg    >  2:439478/1‑151 (MQ=255)                                                                                                                                                      cAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTc  >  2:517843/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       ACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTC  >  NC_000913/3177599‑3177898

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