Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,840,646	0	T	C	44.8%	67.2 / 30.1	59	D365G (GAT→GGT)	nepI	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (29/3);  new base C (0/26);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGCCCC  >  NC_000913/3840496‑3840795                                                                                                                                                          |                                                                                                                                                      cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAAtattat                                                                                                                                                       >  2:944478/1‑151 (MQ=255)    ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCa                                                                                                                                                    <  1:525667/151‑1 (MQ=255)             aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCTCCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAAc                                                                                                                                           <  2:1085194/151‑1 (MQ=255)                     aacaacTACTGACCACCGCTTCTTTCCCTGATACCAGGTTTTCTTCATTTTCCCCTTTACAGTAACCAACAATGCAGTCAGCAACATCAATGTCCCGC‑ACAA‑CATC‑ACCGGCGAAGTCAGCCCAATATTTCCCAGCGCATACCCGCCGAtt                                                                                                                                   <  2:1067348/151‑1 (MQ=255)                        aacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGcc                                                                                                                                >  1:379690/1‑151 (MQ=255)                            actGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCgcgc                                                                                                                            >  1:657315/1‑151 (MQ=255)                             ctGCCCAACGCTTTTTTACCTAATATCAGGATTTCTTCATTTTCCCCTTTGCAGAACCAAACAATGCAGTCAACAACATCAATGTGCCGG‑ACAA‑CACA‑AACGGCGAAGCCAGACAAATATTACCCAGCGCAAAACCGCCGATTGCcgc                                                                                                                              <  2:411055/148‑1 (MQ=255)                             ctGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGcc                                                                                                                           >  1:1114609/1‑151 (MQ=255)                                  cAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAgg                                                                                                                      >  1:1181436/1‑151 (MQ=255)                                  cAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAgg                                                                                                                      >  1:130312/1‑151 (MQ=255)                                    aCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGt                                                                                                                     >  1:352615/1‑150 (MQ=255)                                        ttttttCCCTAATATCAGGTTTTTTTCTTTTCCCCTTTTGAAGAAACAAAAAATGCAGCAACAAACATAAATGGGCCGA‑ACAA‑CACC‑ACCGGGAAAGCCAGACAAATATTCCCCGGCCAAAAACCGCCGATTGCCGCGCCACAGGTATTa                                                                                                                 <  2:1084408/147‑1 (MQ=255)                                        ttcttTACCTGATATCAGTATTTCTTCATTTCCCCTTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGCCAGACCAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTa                                                                                                                 <  1:1129856/150‑1 (MQ=255)                                         tttttACCTGATATCAGGATTTCTTCATTTTCCCCTTGGAAGAAACAAACAATGCAGTCAGCAACATCAATGTGCCGA‑ACAA‑CACC‑AACGGCGAAGTCAGACAAATTTTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAg                                                                                                                <  2:816130/148‑1 (MQ=255)                                            ttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCa                                                                                                              >  1:726385/1‑149 (MQ=255)                                            ttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCa                                                                                                              >  1:772323/1‑149 (MQ=255)                                             tACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGc                                                                                                           >  1:554159/1‑151 (MQ=255)                                             tACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGc                                                                                                           >  1:1165577/1‑151 (MQ=255)                                                cTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTg                                                                                                         >  1:1100801/1‑150 (MQ=255)                                                   atatCAGATTTTCTTTAATTTCACCTTCGCATTAACCAACAGAGAGGTCAGCAACAGCATTGTGCGGA‑ACAA‑CATC‑AACGGCGAATTCAGCCAAATATAACCCACCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCt                                                                                                          <  1:304124/146‑1 (MQ=255)                                                       cAGGTTTTCTTCTTTTTCACTTTTGCAGAAACCAACAATGCAGTCAGCAACATAAATGGGCCGG‑AAAA‑CATC‑AACGGCGAAGTCAGCCCAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAAcc                                                                                                  <  1:557208/150‑1 (MQ=255)                                                           aTTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCAc                                                                                             >  2:863298/1‑151 (MQ=255)                                                               cttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTg                                                                                          >  1:498364/1‑150 (MQ=255)                                                                    ttttCCCCTTTGGAGTAACCAAAAATGCGGCCGGCAACATAAATGTCCGGG‑ACAA‑CACA‑ACCGGCGAGGCCAGCCCAATATCAGCCAGCCCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAAt                                                                                       <  1:109400/148‑1 (MQ=255)                                                                    tttccccttttGCAGTAACAAACAATGCAGTCAGCAAAATCAATGTGCGGA‑ACAA‑CCTC‑ACCGGCGAAGTCAGACAAATATTACCCGCCGAAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTg                                                                                          <  1:628083/137‑1 (MQ=255)                                                                      tccacctgggaattAACCAACAACGCAATCAGAAAAATCAAT‑TGGCGG‑AAAACCATC‑AACGGCGAAGCCAGCCCAATATTATCCAGAGGATAGCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATa                                                                                      <  2:143526/134‑1 (MQ=255)                                                                       tCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGAccc                                                                                 >  1:314161/1‑151 (MQ=255)                                                                        cccttttgaataaacaaacaATGCGGTCAACACCAT‑AATTGGCCGGAACAA‑CATC‑AGCGGCGAAGTCAGCCAAATATTACCCCGCCCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCg                                                                                <  2:943335/139‑1 (MQ=255)                                                                        ccctttggaagaaaCAAAAAATGCGGCCACAAACATCAAGGTGCCGA‑ACAA‑CACC‑ACGGGCGAAGCCAAACAAATATTACCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGAccc                                                                                 <  2:1119114/138‑1 (MQ=255)                                                                         cccTTTGCAAAAACCAAAAACGCAAAAAGCAAAATCAATTTGCGGA‑AAAA‑CTCC‑ACGGGCGAATCCAAACAATTATTACCCAGCCAATACCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGAccc                                                                                 <  1:364715/148‑1 (MQ=255)                                                                            tttGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCtt                                                                            >  1:588938/1‑151 (MQ=255)                                                                             ttGCAGAAACCAAAAAGGCAGTCAGCAACATAAATGTGCCGA‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTACCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCtt                                                                            <  1:243103/150‑1 (MQ=255)                                                                               gCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGTGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCttttt                                                                         >  1:1012829/1‑151 (MQ=255)                                                                                 agaaaCCAACAAGGCAGTCACCAACATAATTGTGCCGG‑AAAA‑CATC‑AACGGCAAAGCCAGCCAAATATTCCCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTc                                                                        <  1:1157417/147‑1 (MQ=255)                                                                                 aGTAACCAACAATGCAGGCAGCAACACAAATGGCCCGG‑AAAA‑CACC‑AACGGCGAACTCAAACCAATATTACCCACCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCt                                                                       <  1:598849/151‑1 (MQ=255)                                                                                   tAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGc                                                                     >  2:523883/1‑151 (MQ=255)                                                                                    aaCCAAAAAGCCGGTCACCAACATCAATGGGCCGG‑AAAA‑CATC‑ACGGGCGAGGTCAAACAAATATTACCCAGCCCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCttttt                                                                         <  1:901150/146‑1 (MQ=255)                                                                                         aaaatgaagccaaaaaCATCAAGGTGCCGG‑A‑AA‑CAACAAACGCCGAGGGAACACAAATATTACCCCGCAAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTg                                                                      <  2:866944/130‑1 (MQ=255)                                                                                           aaTGCAGTCAGCAACATCATTGTGCGGA‑AAAA‑CACC‑AGCGGCGAAGCCAACCAAATTTAACCCAGCCCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCg                                                              <  1:240849/150‑1 (MQ=255)                                                                                                aGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCagc                                                        >  2:733497/1‑151 (MQ=255)                                                                                                       aaCATCAATGTCCCGG‑AAAA‑CATC‑AACGGCAAAGTCAGACAAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCggg                                                 <  1:675763/151‑1 (MQ=255)                                                                                                       aaCATAAAGGGGCCGA‑ACAA‑CACC‑AAGGGCAAAGCCAGCCAAATATTACCCAGCCCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCag                                                         <  2:1078953/143‑1 (MQ=255)                                                                                                        acacCAATGCGCCGA‑ACAA‑CAAC‑AGCGGCGAAGTAAGACAAATATTCCCCAGCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCgg                                                  <  2:337666/149‑1 (MQ=255)                                                                                                        aCATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGt                                                >  1:221160/1‑151 (MQ=255)                                                                                                          aacactgtgCCGG‑ACAA‑AATA‑AGCGCAGAATTCATCCAAATATTCCCCCCCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTg                                               <  2:517814/145‑1 (MQ=255)                                                                                                          aTCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGACTTTTTCTGCCTGATCGGCCAGCGAGCGGGt                                                >  1:337901/1‑149 (MQ=255)                                                                                                             acggcgCCGG‑ACAA‑CAAA‑ACCGACGAAGTCAGACCAATATTCCCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATc                                            <  2:341201/145‑1 (MQ=255)                                                                                                                     gg‑ACAA‑CATC‑ACCGGCGAAGCCAAACAAATATTACCCAGCGCATACCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGtt                                      <  2:286011/148‑1 (MQ=255)                                                                                                                     gg‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGAc                                   >  2:1074513/1‑151 (MQ=255)                                                                                                                     gg‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGAc                                   >  1:1011987/1‑151 (MQ=255)                                                                                                                      gaacaa‑caCC‑ACCGGGGAAGTAAGCCCAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGtt                                      <  1:780875/146‑1 (MQ=255)                                                                                                                        acaa‑caTC‑AACGGCGAAGTCAGACCAATATAATCAAGGGCATAACCTCCCATTGCCGCGCCACAGGTACTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCCGATCGGCCACCGAGCGGGTGATCCACGTTCAc                                   >  2:1894/1‑149 (MQ=255)                                                                                                                         caa‑caTC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGAcc                                  >  1:1195343/1‑149 (MQ=255)                                                                                                                           a‑caTC‑ACCGGCGAAGTCAGACAAATATTACCCAGCCCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGcc                              <  1:453768/151‑1 (MQ=255)                                                                                                                             caTC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCg                             >  2:133221/1‑151 (MQ=255)                                                                                                                              aTC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCg                             >  1:219728/1‑150 (MQ=255)                                                                                                                               tC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGAc                           >  2:67110/1‑151 (MQ=255)                                                                                                                                  aaCGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACgg                         >  2:695589/1‑151 (MQ=255)                                                                                                                                     ggCGAACCCAGACAAATATTACCCAGCCCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACggg                        <  2:232457/149‑1 (MQ=255)                                                                                                                                               gACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCaa             >  2:645040/1‑150 (MQ=255)                                                                                                                                                       ttatCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGcc    >  2:1069197/1‑151 (MQ=255)                                                                                                                                                         atCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGcccc  >  2:897408/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                      CCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGG‑ACAA‑CATC‑AACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGCCCC  >  NC_000913/3840496‑3840795

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