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

CTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCA  >  NC_000913/3177608‑3177891                                                                                                                                              |                                                                                                                                                ctctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc                                                                                                                                       >  1:890529/1‑151 (MQ=255)   ctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccg                                                                                                                                     >  1:679519/1‑151 (MQ=255)    ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                    >  1:209382/1‑151 (MQ=255)    ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                    >  2:970854/1‑151 (MQ=255)    ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                    >  1:727665/1‑151 (MQ=255)    ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                    >  2:758954/1‑151 (MQ=255)     tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                   >  1:900813/1‑151 (MQ=255)          catcCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCAcc                                                                                                                              >  1:741065/1‑151 (MQ=255)                  cAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACt                                                                                                                      >  1:475270/1‑151 (MQ=255)                   acccATCTCAAGTAGCAAAGGGTTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCCCAATCAGCACACCTCCATCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTc                                                                                                                     <  1:325794/149‑1 (MQ=255)                   aaCCATCTCCGGTAGCCAAGGGGTTTCGCTGGTGTCTTCCGCGCCAACCCGAAT‑CTCCTCAATCAGCACAACTTCATCCCGCACGGGGTAAAAGGGGAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACt                                                                                                                      <  1:906863/150‑1 (MQ=255)                         ctcCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCAAAT‑CTGCTCAATCAGCACAACTTCAGCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCAtt                                                                                                                  <  2:366656/148‑1 (MQ=255)                           ccgttAGCCAAGGGGTTTCGCGGGTGTCGTAC‑CCCCCATCCCAATCCTCCCAAATCACCACACCTCCACCACTCCCTAGGTCAAAGGTATGAATAACTGCGGCGTGCCCCCGCTCAAAATTTCCCCGCCGCACCTCATGACTCATTTGccc                                                                                                             <  1:710757/147‑1 (MQ=255)                           ccAGTAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGccc                                                                                                             <  2:913117/151‑1 (MQ=255)                            cAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCg                                                                                                            >  1:180394/1‑151 (MQ=255)                              gTAGCCAGGGGGTTTCGTTGGTGCCGTACGGGGCAACCAGAAT‑CTACTCATTCACCACATCTTAACACCGCTCTGGGTAAAAGGGTGGAAAGACTGCGCCTTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                                                                           <  1:1010692/150‑1 (MQ=255)                                aGCCAGGGGGTTGCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCCCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                         <  1:804458/150‑1 (MQ=255)                                 gCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGa                                                                                                        >  1:99416/1‑150 (MQ=255)                                   cagggggTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGATCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATa                                                                                                     <  2:109142/148‑1 (MQ=255)                                    aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATa                                                                                                     <  2:127703/150‑1 (MQ=255)                                      ggggTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc                                                                                                  <  2:429296/151‑1 (MQ=255)                                           ttCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                              >  2:603749/1‑150 (MQ=255)                                            tCGCTTGTCTCGTGCGGGGCAACCCAAAT‑CTGCAAAATCACCAAAACTCCACCACCCAGTGGCTCAAAGGGTAGCAAGGCTGCGGCGAAACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                             <  1:112813/150‑1 (MQ=255)                                             cGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                             <  2:32635/149‑1 (MQ=255)                                             cGCTGGTGTCGTACGCGGCAACCCGAAT‑CTGCTCAATCAGCACACCTTCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGAc                                                                                            <  2:720218/150‑1 (MQ=255)                                                 ggtgcctaaCGCGGCAACCCAAAT‑CGGCTCAACCAGCACACCTCCACCACCCACGGGGTAAAAGGGTAGCAAGACTGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                                                                                          <  1:700095/140‑1 (MQ=255)                                                       gTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                  >  2:246889/1‑150 (MQ=255)                                                           gcgcaaaTCCAAAT‑CTGCTCAAGCAGCACACCTCCATCCCCCCCGGGAAAAAAGGGCAGCCAGCCGGGGGCGTAACCGCGTTAAAAAATCTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                              <  1:98759/145‑1 (MQ=255)                                                           gcgcCAATCCGAAT‑CTGCCCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCCCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAg                                                                               <  2:483731/149‑1 (MQ=255)                                                           gcgACAATCCGATA‑CTGCCCAAACAGCACATCTTCACCGCGCTGGGGGTAAACGGGGAGAAAGCCGGCGGCGAGACCGCGAGCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                              <  2:876853/150‑1 (MQ=255)                                                             ggaaaTCCGAAT‑CTGCTCAACCAGCACACCTTCACCCCGCACTGGGTCAAGGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt                                                                             <  1:575442/146‑1 (MQ=255)                                                              gCAACCCGAAT‑CTGCTCATCCACCAAAACTTCAACCCGCACTGGGTCAAAGGGTAGAAAGCCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa                                                                          <  2:538104/151‑1 (MQ=255)                                                               cAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa                                                                          >  1:319402/1‑150 (MQ=255)                                                               cAACCCGAAT‑CTGCTCAATCAGCACACCTTCATCCCGCACTGGGTCAAAGGGTAGAAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa                                                                          <  2:243518/150‑1 (MQ=255)                                                                aaTCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa                                                                          <  2:336790/149‑1 (MQ=255)                                                                 aTCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                       >  1:899058/1‑151 (MQ=255)                                                                   ccGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                      >  1:642847/1‑150 (MQ=255)                                                                        t‑ctcCTAAATCGGCACATCTTCATCCCGCCCTGGGTCAAAGGGTAACAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                  <  2:945088/149‑1 (MQ=255)                                                                          cTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                 >  1:459716/1‑149 (MQ=255)                                                                            gCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc                                                             >  2:194093/1‑151 (MQ=255)                                                                             cTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc                                                             <  2:23227/150‑1 (MQ=255)                                                                             cTCAATCAGCACAACTTCACCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAg                                                               <  2:988490/148‑1 (MQ=255)                                                                               cAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                          >  2:206468/1‑151 (MQ=255)                                                                                aaTCAGCACAACTTCACCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                          <  1:707205/150‑1 (MQ=255)                                                                                      cacacctcccccccGCACGGGGTAAAAGGGTAGCAAGACTGCGGGGTCACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAg                                                    <  1:352264/137‑1 (MQ=255)                                                                                         aaCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGtt                                                >  2:613341/1‑151 (MQ=255)                                                                                         aaCTTCACCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtg                                                  <  2:308500/149‑1 (MQ=255)                                                                                            ttCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                               <  2:986413/149‑1 (MQ=255)                                                                                            ttCAACACGCGTTGGACCAAGGGGAAGCAGGTCTGCAGCGGGCCCGGGTAAAAAAATATTCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                                      <  1:179538/142‑1 (MQ=255)                                                                                              catcaCGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtct                                             <  2:727787/149‑1 (MQ=255)                                                                                              catcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCg                                           >  2:1037588/1‑151 (MQ=255)                                                                                              caccccgcccggggTAAAAGGGTAGAAAGCCTGCGGCGTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                            <  2:604797/139‑1 (MQ=255)                                                                                                 caCGCACTGGGTAAAAGGATACCAAGCCTGGGGCGTCCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                            <  1:745698/147‑1 (MQ=255)                                                                                                  aCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                        >  2:428545/1‑150 (MQ=255)                                                                                                   cGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCa                                       >  1:343965/1‑150 (MQ=255)                                                                                                    gCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                     >  2:1035515/1‑151 (MQ=255)                                                                                                         gggTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAtt                                 <  1:876993/150‑1 (MQ=255)                                                                                                                 gggTAGCAAGACTGGGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                         <  2:798221/150‑1 (MQ=255)                                                                                                                       cAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAATTCTATAAAGATCTAATGAAAACAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTg                  >  1:241171/1‑151 (MQ=255)                                                                                                                        aaGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTg                  <  2:376631/150‑1 (MQ=255)                                                                                                                                       ccGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCa  <  2:393204/151‑1 (MQ=255)                                                                                                                                        cgcgCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc   >  1:282729/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                CTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAAT‑CTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCA  >  NC_000913/3177608‑3177891

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