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

GTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCT  >  NC_000913/4332645‑4332924                                                                                                                                                    |                                                                                                                                       gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCAccc‑c                                                                                                                                   >  1:217602/1‑151 (MQ=255) gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCAccc‑c                                                                                                                                   >  2:239317/1‑151 (MQ=255)   tGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCg                                                                                                                                 >  1:9104/1‑151 (MQ=255)    gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGc                                                                                                                                >  2:157529/1‑151 (MQ=255)          tACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTgg                                                                                                                          >  1:398257/1‑151 (MQ=255)           aCCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGc                                                                                                                         >  2:750149/1‑151 (MQ=255)            ccGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCg                                                                                                                        >  2:748430/1‑151 (MQ=255)              ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGtt                                                                                                                      >  1:1024376/1‑151 (MQ=255)                aagaaAACGACTGCCCGGCACGCGCCGGTTGCAGCGGTGGGGAGGCCCTCTCCATCATCTGTTAAGGCCGTGCCACTAACGGTGCTACATCAATGTGATGGGTTTTCGCCCGCATTT‑CCAAATGCAAACGCCCCC‑CCGCCAGTgg                                                                                                                          <  2:238159/145‑1 (MQ=255)                 agaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGc                                                                                                                   >  2:396366/1‑151 (MQ=255)                    aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGt                                                                                                                >  2:49859/1‑151 (MQ=255)                     aaTGCCTGCCCGGCCCGCCCCGTTTGCAGCGGCGGGGAGCCGCTCCCCATCATCTGATCAGCCCGCGCCACTAAGGGTGCTACATAATTGTTATGCTTTTTCCCCAGCAGTT‑CCAGATCAAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGtt                                                                                                               <  1:51773/151‑1 (MQ=255)                          cTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                          >  1:3167/1‑151 (MQ=255)                              ccggccgccgcgggTTGCAGCGGGGGGGACACCCTCTCCATCATCACACAAAGCCGTCCCCCTAACGGTGCACCATCAATGTGGGGTGTTTCCGCCAGCTTTC‑CCGGGTGAAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGa                                                                                                      <  2:276343/138‑1 (MQ=255)                              ccgccccgcgCCGGTTGCAGCGGCGGGGAGCCCCTCTCCATCATCTAATCAAGCCGTCCCACTACCGGTCCTACATCAATGTGAGGCTTTTTCGCCGGCAGTC‑CAAGATGCAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTc                                                                                                              <  2:432378/137‑1 (MQ=255)                              ccgccccgcgCAGGTTGCGGCGGGGGGGAGCCCCTCCCCATCTTTTGACCAAGCGGTCCCACAAACGGTCCTACACCAATGTGATGCTTTTCCCCCAGCAGTT‑CAAGATGAAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTg                                                                                                           <  2:893441/140‑1 (MQ=255)                              ccGGCACGCCCCGTTGGCAGCAGCTGGGAGCCCCTCTCCATCATCTGACCAAGCCGTGCCACTAACGGTGCTACATCAAGGTGATGCTTTTTCGCCGGCAGTC‑CAAGATGAAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTc                                                                                                              <  2:506468/143‑1 (MQ=255)                                ggCACGCGCCAGTTGCAGCAGGGGGGAGCCCCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAGGTGATGCTTTTTCGCCAGCAGTT‑CCAGATGCAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                        <  2:307705/147‑1 (MQ=255)                                  ccggcgcCATTTGCGGCAGCGGGGGGACGCTCCCCATCATCTGATCAGGCCGTGCCACTAGCGGTGCTACATCATTGTGATGCTTTTGCCCCAGGAGTC‑CCAGATGCAACCCCTCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGAc                                                                                                     <  1:698921/145‑1 (MQ=255)                                    cgcgCCGGTTGCAGCGGCGGGGAGACGCTCCCCATCATCTGATCAAGCCGTGCCACTAAGGGTGCTACATCAATGTGATGCTTTTCCGCCAGCAGTT‑CCAGATGCAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACg                                                                                                    <  1:27915/147‑1 (MQ=255)                                       gcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                             >  2:136292/1‑151 (MQ=255)                                        ccggTTGCGGCGGCGGGGAGCCCCTCCCCTTCTTTTGACAAGGCCGTCCCACAAGGGGTGCTACATCAATGTGATGCTTTTTCGCCAGCAGTC‑CCAGATGAAACCGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                            <  2:159430/148‑1 (MQ=255)                                        ccggTTGCAGCGGGGGGGAGCCCCTCTCCATCTTCTGATCAGGCCGGCCCACAACCGGTCCTACATAAATGTGATGCGTTTCCGCCAGCAGTC‑CCAGATGAAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                            <  1:66629/148‑1 (MQ=255)                                        catgtggcggcgcgtGGGATACTCTCTCCATATTCGGAACAGTCCGTGGCCGTAACGGAGGTACATAAATGTTATGGGTTTTCGCCAGAAGTT‑CAAGATCCAAACCAACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                            <  1:206772/137‑1 (MQ=255)                                         cAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                            >  2:977335/1‑150 (MQ=255)                                          aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                           >  2:297705/1‑150 (MQ=255)                                                agcagcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAAcc                                                                                    >  2:259216/1‑151 (MQ=255)                                                  cagcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCtt                                                                                  >  2:425158/1‑151 (MQ=255)                                                  cagcGGGGGGTCCCTCCCCTTCATCGCATCAGGCCTTCCCATAACCGGTGTTACATAAGTGGGTTGCGTTTCCCCCACCATTT‑CCAGAGGCAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaacaac                                                                                     <  1:774353/148‑1 (MQ=255)                                                   agggggaagccGCTCTCCATCATCTGACCAAGCCGTGCCACTAGCGGTGCAACATCAATGTGAGGCGTTTCCGCCAGCAGTT‑CCAGTTGAAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCtt                                                                                  <  1:639861/140‑1 (MQ=255)                                                    ggggggAGCCCCCCTCCATCATCTGACCAAGCCGTCCCACTAACGGTCCTCCATCAATGTGAGGCTTTTCCCCCGGCATTT‑CCAGAGGCAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaaa                                                                                         <  2:690686/139‑1 (MQ=255)                                                    gcggggAGCCGCTCTCCACCTTCTGATCAAGCCGTGCCACTAAGGGTGCTAAATAATGGTGATGCTTTTTCGCCAGAAGTC‑CCAGATGCAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                 <  2:543327/147‑1 (MQ=255)                                                      tGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCg                                                                               >  1:332134/1‑150 (MQ=255)                                                      tGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                              >  2:474499/1‑151 (MQ=255)                                                       gggagcccccccccATCATCTGATCAAGCCGTCCCACAACCGGTCCTACACAAATGGGAGCTTTTTTCGCCAGCGTTC‑CCAGAGGCAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGtt                                                                             <  2:90628/139‑1 (MQ=255)                                                         gagcccctctCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATAAATGTGATGCTTTTCCGCCAGCAGTC‑CAAGATGAAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                           <  2:144272/145‑1 (MQ=255)                                                         gagaCGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                           >  2:157798/1‑151 (MQ=255)                                                             ccctctCCATCATCTGATAAAGCCGTGCCACTAACGGTGCTACATAATTGTGATGCGTTTTCGCCAGCAGTC‑CCAGATGCAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATc                                                                          <  1:738433/146‑1 (MQ=255)                                                             ccctccCCATCATCTGACCAAGCCGTGCCACTAGCGGTGATACAAAAATGTGATGCGTTTCCGCCAGCAGTT‑CGAGATGCAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                       <  2:85522/149‑1 (MQ=255)                                                             ccctccCCATCATCGAATCAAGCCGCCCCCCTACCGGTCCTACATCAATGTGTTGCTTTTTCGCCAGCAGTC‑CGAGATGCAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                           <  1:348735/145‑1 (MQ=255)                                                                      tcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGt                                                              >  1:731988/1‑151 (MQ=255)                                                                       cttcTGATCAGGCCGGCCCACTAAGGGTGCTACATAAATGTGTTGCTTTTTCGCCGGCATTT‑CCAGATGCAACCGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                <  1:893213/146‑1 (MQ=255)                                                                            gataaggCCGTGCCACGAGGGGTGCTACATCAATGTGTTGCTTTTCCGCCAGCAGTC‑CCAGATGCAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                          <  2:84840/143‑1 (MQ=255)                                                                            gATAAACCCGTCCCACTAACGGTGATACAAAAATGTGATGTTTTTTCGCCAGCAGTC‑CAAGATGAAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                         <  2:775575/150‑1 (MQ=255)                                                                             aTCAACCCGTCCCACCAACGGTGCTACACCAATGTGATGTTTTTCCGCCAGCAGTC‑CAAGATGAAACCGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                         <  2:988136/149‑1 (MQ=255)                                                                                   ccGTCCCACAAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCATTT‑CCAGATGCAACCGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGa                                                       <  1:329120/145‑1 (MQ=255)                                                                                     gTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                 >  2:405594/1‑149 (MQ=255)                                                                                     gTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                 >  1:472624/1‑149 (MQ=255)                                                                                      tccctcTAACGGTGCAACATAAAGTTGATCCTTTTTCGCCAGCATTC‑CCAGATGAAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                      <  1:73405/138‑1 (MQ=255)                                                                                        ccACTAGCGTTTTTACATCATTGTGTTGCTTTTTCGCCAGCTGTC‑CCAGGCGCAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGAt                                             <  2:960773/150‑1 (MQ=255)                                                                                         cACTAGCGGTGTTACATAAATTTGGTCTGTTTTCCCCAG‑AGTTCCCAGACGCAATCGCCCCG‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                         <  1:961021/137‑1 (MQ=255)                                                                                                ggTGCTACATAAGTGTGTTGGTTTTCCGCCAGCATTT‑CCAGAGAAAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACg                                    <  2:855853/151‑1 (MQ=255)                                                                                                      aCATCAATGTGATGCTTTTTCCCCAGCATTT‑CCAGATGCAAACGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                 <  1:157534/148‑1 (MQ=255)                                                                                                       cATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                             >  2:131773/1‑151 (MQ=255)                                                                                                        aTCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCTCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGCAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGCTTAACGCCGAAacc                            >  1:804201/1‑151 (MQ=255)                                                                                                             tgtgATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                       >  1:250787/1‑151 (MQ=255)                                                                                                                gATGCGTTTCCGCCGGCAGTC‑CGAGATGCAACCGCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                      <  1:543757/149‑1 (MQ=255)                                                                                                                     gTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt               >  1:35373/1‑151 (MQ=255)                                                                                                                     gTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt               >  2:850401/1‑151 (MQ=255)                                                                                                                      ttttCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc              >  2:261609/1‑151 (MQ=255)                                                                                                                      tttcccccAGCTGTC‑CCAGATGAAACCCCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTc                    <  2:680540/139‑1 (MQ=255)                                                                                                                      tttcccccAGCAGTT‑CCAAATGCAACCCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAt                 <  1:793355/142‑1 (MQ=255)                                                                                                                      tttcCGCCGGCATTC‑CCAGATGAAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAt                 <  1:23626/148‑1 (MQ=255)                                                                                                                      tttcCGCCAGCAGTC‑CAAGATGAAAACCCCCCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt               <  1:1005722/150‑1 (MQ=255)                                                                                                                       tttCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg             >  1:141698/1‑151 (MQ=255)                                                                                                                            caagacgTT‑CCAGCT‑CAAACCCCCCCACCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGa            <  1:168580/141‑1 (MQ=255)                                                                                                                              agcagTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTgg      >  1:943701/1‑151 (MQ=255)                                                                                                                                 agTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc   >  2:621139/1‑151 (MQ=255)                                                                                                                                  gTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGGAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCt  >  1:389532/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                       GTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTT‑CCAGATGCAAACGCACCC‑CCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCT  >  NC_000913/4332645‑4332924

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