Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,791	0	A	C	43.9%	94.2 / 17.8	57	V161G (GTG→GGG)	basS	sensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (32/0);  new base C (0/25);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.01e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGC  >  NC_000913/4332641‑4332933                                                                                                                                                        |                                                                                                                                               aCATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCa                                                                                                                                                <  1:152215/151‑1 (MQ=255)     gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                            >  1:996210/1‑151 (MQ=255)               aCCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGc                                                                                                                                  >  1:946501/1‑151 (MQ=255)                  ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCCGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGtt                                                                                                                               >  2:95821/1‑151 (MQ=255)                     agaaAATGCCTGCCCGGCGCGCGCCAGTTGAAGCAGCGGGGAGACGCTCT‑CCATCACCTAATCAACCGGTGCCACTAACGGTGCTACATCAAGGTGATGCTTTTTCGCCAGCAGTTCCAGAGGCAACCCCCCCCCCGCCAGTGGCGTTCGc                                                                                                                            <  2:273574/151‑1 (MQ=255)                       aaaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                            >  1:194928/1‑149 (MQ=255)                        aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGt                                                                                                                         >  2:1025573/1‑151 (MQ=255)                           tGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                      >  2:307643/1‑151 (MQ=255)                            gACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGt                                                                                                                     >  2:1045017/1‑151 (MQ=255)                              cTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                                   >  1:728504/1‑151 (MQ=255)                               tGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgcg                                                                                                                  >  1:8193/1‑151 (MQ=255)                                gTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                 >  1:489078/1‑151 (MQ=255)                                       aCGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                           >  2:962694/1‑150 (MQ=255)                                         gcacgaGTAGCAGCAGGTGGGAGACGCTCC‑CCTTCTTCGAATCAGCCCGTGCCATTAATGTTGTTACACAAATGTGTTGCTTTTCCGCCAGCAGTCCAAGACGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGAc                                                                                                              <  1:981926/140‑1 (MQ=255)                                          cgcCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCAAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCg                                                                                                       >  2:793496/1‑151 (MQ=255)                                            cgaggcggggcagcGGGGAGCCG‑TCTCCCACCATTGGACCAAGCCGTCCCACTAACGGTGCTCCACCAATGGGTTGTGTTTCCGCCGGCGGTCCGAAGTGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                           <  2:124930/136‑1 (MQ=255)                                            ccAGTGGCAGCAGCGGGGAGCCCCTCC‑CCATAATCTAACCAAGCCGTGCCACTAAGGGTGCTACATCAATGTGATGTTTTTTCCCAAGCAGTTCCGGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                     <  1:582894/151‑1 (MQ=255)                                                ttGCAGCGGGGGGGAGCCCCTCC‑CCATCATCTAATCAAGCCGGCCCACAAACGGTGCTACATAAATGTGGTGCTTTTCCGCCAGCAGTCCCAGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                         <  1:672376/143‑1 (MQ=255)                                                tggcagcaGCGGGGAGCCCCTCT‑CCATCATCTGATCAAGCCGTGCCACTAAGGGTGCTACATAAATGTGATGCTTTTCCGCCAGCATTTCCAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaa                                                                                                   <  1:21760/147‑1 (MQ=255)                                                   ctgcctcTGGAAGTCACTCT‑CCATCTTTGGCTCAAGCCGCCCCACTAACTGTGCTCCCTCTATGTCCTTCTTTTTCCCCAGCAGTCCAATCCGAAAGCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCAACGTCAGCGGTAaacaa                                                                                               <  2:259452/144‑1 (MQ=255)                                                       agcTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                          >  1:448526/1‑151 (MQ=255)                                                                 cGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                <  1:622287/151‑1 (MQ=255)                                                                 cGCTCC‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACAACAATGTGATGCGTTTCCGCCAGCAGTTCCAGAGACAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                                    <  1:388305/147‑1 (MQ=255)                                                                   ctct‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                              >  1:314608/1‑151 (MQ=255)                                                                   ctct‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                              >  1:195649/1‑151 (MQ=255)                                                                     ct‑cCATCTTCTGAACAACCGGTGCCACAAACGTGTTTACATCAGTGTGTTGGTTTTTCGCCAGCAGTCCCAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                <  2:835382/147‑1 (MQ=255)                                                                     ct‑cCATCATCTGATCAAGCCGGGCCACTAACGGTGCTACATAAATGTGATGCTTTTCCGCCAGCAGTCCCAGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGt                                                                             <  1:463284/150‑1 (MQ=255)                                                                        ccATCATCTGATCAAGCCGTGCCACTAACGGTACTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGc                                                                           >  1:956814/1‑150 (MQ=255)                                                                            cgtcTGATCAGCCCTGGCCACTACGGGGCCTACACAAATGTGAGCCGTTTCGGCCGGCATTCCCGGGGGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTgg                                                                        <  2:427077/147‑1 (MQ=255)                                                                                tGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                   >  2:964650/1‑150 (MQ=255)                                                                                 gATCAGGCCTTGCCACTACCGGTGCTACTTCAATTTGTGGCGTTTTCGCCAGCTGTTCCAGACCCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                  <  2:1008214/150‑1 (MQ=255)                                                                                    cAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                              >  1:861496/1‑151 (MQ=255)                                                                                     aaGCCGTGCCACTACCGGTGCTACATCAAGGTGATGCTTTTTCGCCAGCAGTTCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                              <  1:727637/150‑1 (MQ=255)                                                                                     aaGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                              <  1:350105/150‑1 (MQ=255)                                                                                     aaGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                             >  1:783215/1‑151 (MQ=255)                                                                                             ccACTAACGGTGCAACATCAAGGGGATGCGTTTTCGCCAGCAGTCCGAGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                     <  2:612269/151‑1 (MQ=255)                                                                                                cTAACGGTTCTACAACAATGTGAGGCGTTTTCCCCACCAGTACCAAATACAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                   <  1:615091/150‑1 (MQ=255)                                                                                                  aaCGTTGCTACATAAATGTGATGCTTTTCCGCCAGCAGTCCAAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTa                                                <  1:434380/151‑1 (MQ=255)                                                                                                    cGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAAc                                              >  2:297394/1‑151 (MQ=255)                                                                                                     ggTGCTACATCAATGTGATGCTTTTTCGCCAGCAGTTCCAGAGCCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                             <  1:1045673/151‑1 (MQ=255)                                                                                                      gTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                            >  2:1031408/1‑151 (MQ=255)                                                                                                          tACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                         >  1:921254/1‑150 (MQ=255)                                                                                                                  tgtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                                >  1:659001/1‑151 (MQ=255)                                                                                                                       tGCGTTTCCGCCAGCAGTTCCAGAGGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                            <  1:73041/150‑1 (MQ=255)                                                                                                                                 ccAGCATTCCCAGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGa                     <  1:1038138/147‑1 (MQ=255)                                                                                                                                 ccAGCAGTTCCAGATCCAAACGCACCCCCGCCAGTGGCGTTCGCAGTACGTGCGCGACGTCAGCGGTAAACAACCTCTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCATATCTTTCACCACCCCAACCACCGCTTCGactcacaac                   >  2:855888/1‑141 (MQ=255)                                                                                                                                   agaagTTCAAGAGGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                  <  2:734101/145‑1 (MQ=255)                                                                                                                                      agTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc            >  2:357496/1‑151 (MQ=255)                                                                                                                                       gTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCt           >  1:845621/1‑151 (MQ=255)                                                                                                                                        ttCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg          >  2:667482/1‑151 (MQ=255)                                                                                                                                          ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg        >  1:423140/1‑151 (MQ=255)                                                                                                                                           cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg        >  1:49628/1‑150 (MQ=255)                                                                                                                                            aGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGa       >  2:194592/1‑150 (MQ=255)                                                                                                                                              ctgcagccGCCCCGCCGCCAGTGGCGTTCGAAGTTCGTGCGCGACGTCAGAGGTAAACAGCCTTTCGTGAGCCAGCGTGCTGGTTGGGCGTCTGACCAGATCGTTTAACGCCGAATCCACCGCTTCGATTTCGAGTGTTGCGCtgtg        <  2:28706/140‑1 (MQ=255)                                                                                                                                               tGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATgg   >  1:425079/1‑151 (MQ=255)                                                                                                                                                gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc  >  2:773571/1‑151 (MQ=255)                                                                                                                                                  aaccGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGa       <  1:918716/141‑1 (MQ=255)                                                                                                                                                  aaccGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc  <  1:37674/146‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               ACATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCT‑CCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGC  >  NC_000913/4332641‑4332933

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