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

ACATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332641‑4332941                                                                                                                                                       |                                                                                                                                                       aCATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCa                                                                                                                                                        <  2:354874/151‑1 (MQ=255)  cATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAc                                                                                                                                                       >  2:836559/1‑151 (MQ=255)     gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                                    <  2:458961/151‑1 (MQ=255)     gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                                    >  1:198831/1‑151 (MQ=255)        gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCg                                                                                                                                                  <  2:246629/150‑1 (MQ=255)        gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGc                                                                                                                                                 >  1:222004/1‑151 (MQ=255)           aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGcc                                                                                                                                                >  1:405660/1‑149 (MQ=255)           aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGcc                                                                                                                                                >  1:718354/1‑149 (MQ=255)           aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCa                                                                                                                                               <  1:11890/150‑1 (MQ=255)             ttACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTg                                                                                                                                            >  2:166380/1‑151 (MQ=255)                  ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGtt                                                                                                                                       >  2:3231/1‑151 (MQ=255)                    aagaaAATGCCTGCCCGGCCCGCGCCTGTGGCAGCGGCGGGGAGCCGCTCTCCACCATCTGACCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCCCCAGCAGTTCCAGATGCAACCCCCCCCCCGCCAGTGGCGtt                                                                                                                                       <  2:280098/149‑1 (MQ=255)                    aagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCg                                                                                                                                     >  2:970061/1‑151 (MQ=255)                    aagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCg                                                                                                                                     >  2:564310/1‑151 (MQ=255)                    aaaaaaaTGACTGTCCGGCACGCGCCGGTTGCAGCAGCGGGGAGACCCTCTCCATCATCTGACCAAGCCGTGCCACTAACGGTGCTCCATAAATGTGATGCTTTTTCGCCAGCATTTCAAGATGCAACCGCCCCCCCGCCAGTGGCg                                                                                                                                         <  2:394410/144‑1 (MQ=255)                      gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                    >  2:410885/1‑150 (MQ=255)                        aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGGGGGGAGCCCCTCCCCATCATCTGATCAAGCCGTGCCACTACCGGTGCTACATCAATGTGAGGCGTTTTCCCCAGCAGTTCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGt                                                                                                                                 <  2:25884/151‑1 (MQ=255)                        aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGt                                                                                                                                 >  2:258028/1‑151 (MQ=255)                         aaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                                >  2:522587/1‑151 (MQ=255)                            gACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGt                                                                                                                             >  2:739401/1‑151 (MQ=255)                                      cccgccCCAGTTGCAGCAGCGGGGAGACCCTCTCCATCTTCTGATCAAGCCTTGCCACTACGGGGGCTACATCATTGTGATGCGTTTTCCCCGGCAGTCCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGt                                                                                                                    <  1:694926/148‑1 (MQ=255)                                       aCGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                   >  1:817731/1‑150 (MQ=255)                                            ccAGTTGCAGCAGCGGGGAGACGCTCCCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCATTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                             <  1:436796/151‑1 (MQ=255)                                                   cagcggCGGGGAGCCGCTCCCCATCATCTGACCAAGCCGTCCCACAACCGGTGCAACAACAATGTGATGCTTTTTCCCCAGCATTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaacaa                                                                                                       <  2:708559/150‑1 (MQ=255)                                                       agcTGGGAGCCCCTCTCCATCATCTAATCAACCCGTCCCACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTCCCAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCtt                                                                                                   <  1:854315/150‑1 (MQ=255)                                                           gggATACTCTCCCCATCTTCGTATCAAGCCGTGCCACTACCGGAGTTACATAAATGGTTTGTTTTTTCCACAGCAGTCCCAGAAAAAAACCCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                               <  2:595215/150‑1 (MQ=255)                                                                  gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                       >  1:133090/1‑151 (MQ=255)                                                                   ctctCCACCACCTGATCAAGCCGTGCCACTAACGGTGCTACATAAATGTGATGCTTTTCCGCCAGCAGTCCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                       <  1:232398/150‑1 (MQ=255)                                                                    tctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                                      >  2:952271/1‑150 (MQ=255)                                                                       ccaccaTGTGACCAAGCCGTCCCACTAACGGTGCTACATAAATGTGTTGCGTTTTCCCCAGCAGTCCCAATTGCAATCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                       <  1:404064/146‑1 (MQ=255)                                                                           catcTGATCAAGCCGTGCCACTACCGGTGCTACATAAATGTGATGGTTTTTCCCCAGCATTTCAAGAGGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTg                                                                                    <  1:543813/145‑1 (MQ=255)                                                                            atcCGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                              >  2:586351/1‑150 (MQ=255)                                                                                    aaGCCGTGCCACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                          <  1:519667/146‑1 (MQ=255)                                                                                       ccGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                  >  1:102401/1‑151 (MQ=255)                                                                                       ccGTCCCCCTAACGGTGCTACATCAATGTGAGGCTTTTCCGCCAGCAGTTCAAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGa                                                                    <  1:490230/149‑1 (MQ=255)                                                                                       ccGTCCCACTAACGGTCCTACATCAAGGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAc                                                                   <  1:792994/150‑1 (MQ=255)                                                                                       ccGCGCCCCTACCGGTGCTACATCAAGGTGTGGCTTTTCCCCCAGCAGTCCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                       <  2:505440/146‑1 (MQ=255)                                                                                            ccACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                  >  1:358630/1‑146 (MQ=255)                                                                                                  aCGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTa                                                        >  1:417348/1‑150 (MQ=255)                                                                                                   cGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAAc                                                      >  1:42276/1‑151 (MQ=255)                                                                                                    ggtccaacacCATTGTGATGGTTTTCCGCCAGCAGTCCAAGAGGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                                     <  2:30841/145‑1 (MQ=255)                                                                                                    ggTGCTACATCAATGTGTTGCGTTTTCGCCAGCATTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                                     <  2:568127/151‑1 (MQ=255)                                                                                                               attgtgATGCTTTTCCCCCAGCAGTTCCAAATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                             <  2:103731/146‑1 (MQ=255)                                                                                                               aaTGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                          >  1:605430/1‑151 (MQ=255)                                                                                                                aTGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCg                                         >  1:135964/1‑151 (MQ=255)                                                                                                                    gATGCTTTTCCCCCAGCAGTCCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                      <  2:827148/150‑1 (MQ=255)                                                                                                                      tGCGTTTTCGCCAGATTTTCCAGAAGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCg                                         <  2:202259/145‑1 (MQ=255)                                                                                                                       gctttttCCCCAGCAGTTCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGa                                   <  1:696100/147‑1 (MQ=255)                                                                                                                          tttcCGCCAGCGTTTCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                                       <  1:999208/143‑1 (MQ=255)                                                                                                                           tttCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                              >  2:383052/1‑151 (MQ=255)                                                                                                                               ccatgcttttCCAAACGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                          <  1:281531/143‑1 (MQ=255)                                                                                                                                ccggcattcccAGTTGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                          <  2:499054/141‑1 (MQ=255)                                                                                                                                         ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                >  1:193961/1‑151 (MQ=255)                                                                                                                                             aTGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg            <  2:691484/151‑1 (MQ=255)                                                                                                                                              ggaaaccGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGa               <  1:621312/141‑1 (MQ=255)                                                                                                                                               gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc          >  2:621550/1‑151 (MQ=255)                                                                                                                                                cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc          >  1:563369/1‑150 (MQ=255)                                                                                                                                                       aCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCg  >  1:962154/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       ACATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332641‑4332941

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: ‑