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

ACATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCA  >  NC_000913/4332641‑4332934                                                                                                                                                        |                                                                                                                                                aCATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCa                                                                                                                                                 <  1:700063/151‑1 (MQ=255)   aTGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCAcc                                                                                                                                               >  1:573253/1‑151 (MQ=255)   aTGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCAcc                                                                                                                                               <  1:410645/151‑1 (MQ=255)     gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCAACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCAcccc                                                                                                                                             <  1:442136/151‑1 (MQ=255)         aTAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGcc                                                                                                                                         <  2:498390/151‑1 (MQ=255)                cggaaagaaAATGCCTGCCGGGCACGCGCCGGTGGCAGCGGGGGGGAGCCCCTCCCCACCATCTGACCAAGCCGTGCCACTAACGGTGCTACATCAATGGGATGCGTTTCCGCCAGCATTTCC‑AGATGCAAACGCCCCCCCGCCAGTg                                                                                                                                     <  2:981060/144‑1 (MQ=255)                     aaaaaaTGCCTGCCCGCCCCGCGCCAGTTGCAGCAGCGGGGAGCCCCTCCCCATCATCTGATCAAGCCGTCCCACAACCGGTGCTACATCAATGGGATGCGTTTCCGCCAGCAGTCCA‑AGATGAAAACCCCCCCCCGCCAGTGGCGTTCGc                                                                                                                             <  2:983521/149‑1 (MQ=255)                      gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                             >  1:1086003/1‑150 (MQ=255)                       aaaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                             >  1:982628/1‑149 (MQ=255)                       aaaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                             >  2:805570/1‑149 (MQ=255)                                gTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                  >  1:943587/1‑151 (MQ=255)                                  ccGGCCCGCGCCAGTTGCAGCAGCGGGGAGCCGCTCTCCATCATCTGATCAAGCCGTGCCCCTAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCAGTCCC‑AGATGAAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTg                                                                                                                     <  1:837645/146‑1 (MQ=255)                                  ccGGCACGCGCCAGTTGCAGCAGCGGGGAGCCGCTCTCCATCATCTGACCAAGCCGTGCCACTAACGGTGCTACATAAATGTGATGCGTTTTCGCCAGCAGTCCC‑AGAGGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                  <  1:632549/149‑1 (MQ=255)                                    gccccgcgCCGGTGGCAGCGGCGGGGAGCCGCTCCCCATCATCTGACCAAGCCGTGCCCCTAACGGTGCTCCATAATTGTGATGCGTTTTCGCCAGCAGTTCC‑AGAGGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                                    <  2:1187255/141‑1 (MQ=255)                                          cgcCAGTTGCAGCGGGGGGGACCCGCCCTCCATCTTCGCATCAGGCCGTCCCACTAACGGTGCTACATCAATATAATGCTTTTCCCCCGGCATTTCC‑AGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCg                                                                                                        <  2:684001/151‑1 (MQ=255)                                          cgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCg                                                                                                        >  1:691786/1‑151 (MQ=255)                                           gcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                                       >  2:1202816/1‑151 (MQ=255)                                           cccggtggcagcgGCGGGGAGCCCCTCCCAATCATCTGACCAACCCGGCCCACTACCGGTGCTACATCAATGTGAGGCTTTTCCGCCAGCAGTTCC‑AGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                                       <  2:844792/144‑1 (MQ=255)                                            ccAGTTGCAGCAGCTGGGAGACCCTCCCCTTCATCTGATCAAGCCGTGCCACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTCCC‑AGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGt                                                                                                             <  1:263216/144‑1 (MQ=255)                                               ttggcagcgGCGGGGAGCCCCTCTCCATCATTTGATCAACCCGTCCCCCTACCGGGGCTACATCAAGGGGATGTTTTTTCCCCAGCTTTCCC‑AAACGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                      <  1:674369/145‑1 (MQ=255)                                                   cagcggggggggggcGCTCCCCACCATCGGACCAAGCCGCCCCACGAACGGCGCCCCATAATGGGGTTGTGTTTTCGCCAGCAGTACC‑AGACGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                                       <  2:905225/129‑1 (MQ=255)                                                    agcagcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAAcc                                                                                              >  2:93084/1‑151 (MQ=255)                                                    agcagcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAAcc                                                                                              >  1:33120/1‑151 (MQ=255)                                                        gcgggggGACGCTCCCCTTCATCTGATCAAGCCGCCCCACTAAGGGTGCTACAACAGTGTGATGCTTTTCCGCAAGCAGTTCC‑AGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                           <  2:220251/147‑1 (MQ=255)                                                        gcggggacgcGCTCCCCATCACCTGACCAGCCGGGGCCACTACGGGTGTGACATCAATTGGCTGGGTTTTCGCCGGGATTCCA‑AGTTAAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCtt                                                                                            <  2:1194381/140‑1 (MQ=255)                                                          ggggAGCCGCCCCCCATCATTTGACCAAGCCGTCCCAATAACGGTGATACATAAATGTGATGTTTTTCCGCCACCAGTCCC‑AGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                        <  1:1227428/150‑1 (MQ=255)                                                          ggggAGCCCCTCCCCTTCATCTGACCAAGCCGTGCCACTAACGGTGCTACATCAATGTGTTGCGTTTCCGCCAGCAGTCCC‑AGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                        <  2:841933/150‑1 (MQ=255)                                                                   ctccCCATCTTCTGACCAAGCCGTGCCACTAACGGTGCTACACAAATGTGATGCTTTTTCGCCAGCATTTCC‑AGAGGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                <  1:1207720/150‑1 (MQ=255)                                                                    tctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                               >  1:126243/1‑150 (MQ=255)                                                                        catcatCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                          >  1:101358/1‑151 (MQ=255)                                                                          tcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGt                                                                        >  1:1195750/1‑151 (MQ=255)                                                                            atcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                       >  1:889247/1‑150 (MQ=255)                                                                            atcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                       >  1:978049/1‑150 (MQ=255)                                                                              ctgacaaggccggcccACAAAGGGTGCTACATAAATGTGATGCTTTTCCGCCAGCAGTCCC‑AGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGt                                                                              <  2:78674/127‑1 (MQ=255)                                                                              cTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                    >  1:968321/1‑151 (MQ=255)                                                                                 aTCACGCGGTGCCCATAACGGTGCTACATAATTGTGTTGCTTTTCCCCCAGTATTCCC‑ATATCCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAg                                                                     <  1:781360/147‑1 (MQ=255)                                                                                  tCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                >  2:579859/1‑151 (MQ=255)                                                                                   cAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                               >  1:1180812/1‑151 (MQ=255)                                                                                    aaGCCGTGCCACTAACGGTGCTACATAAATGTGATGCGTTTCCCCCAGCAGTTCA‑AGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGa                                                                 <  1:750139/148‑1 (MQ=255)                                                                                    aaGCCGTCCCACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCAGCATTCCC‑AGAGGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                    <  1:1058244/145‑1 (MQ=255)                                                                                       ccgccccACTACCGGTCCTAAATAATTGTGATGCGTTTTCCCCAGCAGTTCC‑AGATCCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                  <  2:819975/139‑1 (MQ=255)                                                                                       ccGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                           >  2:1194057/1‑151 (MQ=255)                                                                                        cGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                           >  2:809348/1‑150 (MQ=255)                                                                                         gTGCCACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCGGTCCC‑AGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                   <  2:1050840/141‑1 (MQ=255)                                                                                         gTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                           >  2:543040/1‑149 (MQ=255)                                                                                            ccACAAAGGGTGCAACAACAGTGGGTTGCGTTTCCGCC‑GGTGTTCCAAGAGGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                      <  2:620848/151‑1 (MQ=255)                                                                                             cACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                           >  1:680674/1‑145 (MQ=255)                                                                                              aCTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                    >  2:621784/1‑151 (MQ=255)                                                                                                      tGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                            >  2:522040/1‑151 (MQ=255)                                                                                                       gCTACATAAATGTGTTGGTTTTTCCCCAGCAGTTCA‑AAATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                           <  1:275968/151‑1 (MQ=255)                                                                                                          acataaatggattGCGTTTTCGCCGGCAGTTCA‑AAAGCAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                           <  2:1251924/136‑1 (MQ=255)                                                                                                                aTGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCg                                  >  2:1096077/1‑151 (MQ=255)                                                                                                                 tgtgATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                                 >  2:209755/1‑151 (MQ=255)                                                                                                                  gtgATGCTTTTTCGCCAGCAGTCCC‑AGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                                <  1:128170/151‑1 (MQ=255)                                                                                                                  gtgATGCGTTTCCGCCAGCAGTTCC‑AGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                   <  1:997101/148‑1 (MQ=255)                                                                                                                      tGCTTTTCCGCCAGCAGTTCC‑AGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                             <  1:426773/150‑1 (MQ=255)                                                                                                                        cGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                          >  1:84511/1‑151 (MQ=255)                                                                                                                          tttcccccAGCATTTCC‑AGATCAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTc                              <  2:549137/139‑1 (MQ=255)                                                                                                                          tttcCGCCAGCAGTCCC‑AGACGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                                <  2:1163268/143‑1 (MQ=255)                                                                                                                           tttCGCCGGCATTTCA‑AGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                        <  2:415111/150‑1 (MQ=255)                                                                                                                              cGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAgg                    >  2:578865/1‑151 (MQ=255)                                                                                                                                ccAGCATTCCC‑AGATCCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAg                     <  1:879733/148‑1 (MQ=255)                                                                                                                                ccAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                  >  2:1153740/1‑151 (MQ=255)                                                                                                                                 caccagTCCA‑AGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                  <  2:968312/147‑1 (MQ=255)                                                                                                                                  agaggTTCC‑AGACTCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGa                      <  2:723409/141‑1 (MQ=255)                                                                                                                                      gTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCt            >  1:196216/1‑151 (MQ=255)                                                                                                                                         cc‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg         >  1:566567/1‑151 (MQ=255)                                                                                                                                                 cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc   >  2:699529/1‑150 (MQ=255)                                                                                                                                                 cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc   >  1:815056/1‑150 (MQ=255)                                                                                                                                                 cAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc   >  1:624017/1‑150 (MQ=255)                                                                                                                                                  aaaCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCa  >  1:902440/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                ACATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCA  >  NC_000913/4332641‑4332934

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