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

ATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGG  >  NC_000913/4332643‑4332938                                                                                                                                                      |                                                                                                                                                    aTGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcc                                                                                                                                                   >  2:1563246/1‑151 (MQ=255)   gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                                 >  1:332267/1‑151 (MQ=255)   gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                                 <  2:1521866/151‑1 (MQ=255)      gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAccccc                                                                                                                                                <  2:1157069/149‑1 (MQ=255)      gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGc                                                                                                                                              <  2:416381/151‑1 (MQ=255)          aTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGt                                                                                                                                          >  2:734071/1‑151 (MQ=255)           tttccGGAAAAAAATGACTGTCCGGCCCGCGCCGGTTGCAGCAGCTGGGAGCCTCTCTCCTTCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCGGATGCAAACGCCCCCCCGCCa                                                                                                                                            <  1:1458274/145‑1 (MQ=255)              cggaaaaaaaaTGCCTGCCCGCCACGCCCCGGTTGCAGCGGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTCCCACTACCGGTG‑CTACATAAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAACCCCCCCCCCGCCAGTg                                                                                                                                         <  1:669924/141‑1 (MQ=255)                ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGtt                                                                                                                                    >  1:433422/1‑151 (MQ=255)                  aagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCg                                                                                                                                  >  2:1641884/1‑151 (MQ=255)                    gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                 >  1:82176/1‑150 (MQ=255)                     aaaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                 >  1:1488583/1‑149 (MQ=255)                      aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAg                                                                                                                               >  1:1499267/1‑150 (MQ=255)                      aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGt                                                                                                                              >  1:1299275/1‑151 (MQ=255)                       aaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                             >  2:1632934/1‑151 (MQ=255)                        aTGTCTGCCCGGCACGCGCCAGTGGCAGCAGGGGGGAGCCGCTCTCCATCATCTGATCAAGCCGTGCCACAAAGGGTC‑CTACATCAATGTGATGTGTTTTCGCCAGCATTTCCAAATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                            <  1:1658925/151‑1 (MQ=255)                         tGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                           >  2:253670/1‑151 (MQ=255)                         tGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                           >  1:41715/1‑151 (MQ=255)                          gacggccgggCCCGCGCCAGTGGCAGCAGGGGGGAGACGCTCCCCATCTTCTGAACAAGCCGTCCCCCTAACGGTG‑CTACATCAATGTGATGCTTTTCCGCCAGCAGTCCAAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                           <  2:798647/142‑1 (MQ=255)                           aCTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTg                                                                                                                         >  2:1400179/1‑151 (MQ=255)                                ccGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGcgcg                                                                                                                     >  2:1027240/1‑150 (MQ=255)                                  ggCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACg                                                                                                                  >  1:516217/1‑151 (MQ=255)                                    ccggcgcgttttGAAGCAGTGGGGGGCCCCTCCCCATCATCTCATCAGGCCGTGCCACTACGGGTC‑CTACATCATGGTGGTGCGTTTTCCCCAGCTTTTCAAAATCCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGa                                                                                                                    <  2:1699019/137‑1 (MQ=255)                                    cccgccCCTGTTGCAGCAGCGGGGAGCCGCTCCCCATCTTCTGACCAACCCGTGCCACTACGGGTC‑CTACATCATTGTGATGCTTTTCCCCCAGCATTCCCAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                <  2:962493/149‑1 (MQ=255)                                    cACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                >  1:1090811/1‑151 (MQ=255)                                    cACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                >  1:1426/1‑151 (MQ=255)                                      cgccccgtttGCAGCAGCGGGGAGCCCCTCCCCTTCATCTAATAACGCCGTCCCACAAACGGTC‑CTACATCAATGTGATGCTTTTCCGCCAGCAGTCCCAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                <  2:56107/141‑1 (MQ=255)                                        cgcCAGTTGCAGCAGGGGGGAGCCGCTCTCCATCATCTGATCAAGCCGGGCCACTAAGGGTG‑CTACATAAATGTGATGTTTTTCCGCCAGCAGTTCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                              <  2:955247/149‑1 (MQ=255)                                         gcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                                           >  1:244253/1‑151 (MQ=255)                                          ccggtggcagcgGCTGGGAGCCGCTCTCCATCATCTGACCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCTTTTTCGCCAGCATTTCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGt                                                                                                                          <  1:1605711/129‑1 (MQ=255)                                          ccggTTGCAGCGGCGGGGGGCCCCTCCCCATCACTTAAAAAACCCGTGCCACTAACGGTG‑CTAAACAAATGTGATGTTTTTCCGCCGGCGGTCCCAGAGCCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                              <  1:1089041/144‑1 (MQ=255)                                                   gcagcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCt                                                                                                 >  1:1258934/1‑151 (MQ=255)                                                     ggtggggAGCCCCTCCCCATCATCTGATCAAGCGGTGCCA‑TAAGGGGGCCAACATAAATGGGTTGTTTTTCCGCCAGCATTTCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                               <  2:421773/147‑1 (MQ=255)                                                      ggggggAGCCCCCCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACACCAATGTGATGCGTTTCCGCCGGCATTCCCAGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                     <  1:172173/141‑1 (MQ=255)                                                        tGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                            >  2:971965/1‑151 (MQ=255)                                                        gggggGACCCTCTCCACCATCTGATCAAGCCGTGCCACTACGGGTG‑CTACATCAATGTGATGCGTTTCCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                            <  2:1159353/150‑1 (MQ=255)                                                        ggggAGCCCCTCCCCATCATCTGAACAGCCCGGGCCACTAACGGTG‑CAACATAAATGTGATGCTTTTCCGCCAGCATTCCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCg                                                                                             <  1:1446424/149‑1 (MQ=255)                                                         gggAGACCCTCCCCATCATCTGACCAAGCCGTGCCAAAAACGGTG‑CTACATCAATGTGATGCTTTTCCGCCAGCAGTCCCAGAGGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                               <  2:592254/147‑1 (MQ=255)                                                               ccccctcctccatcTGACCAGCCCGTCCCACTACCGGTC‑CTACATAAATTTGATGCTTTTCCGCCACCAGTTCCAGATCAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCa                                                                                      <  1:647792/140‑1 (MQ=255)                                                                  tctcCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGt                                                                                  >  2:1694791/1‑151 (MQ=255)                                                                      catcatCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                              >  1:724340/1‑151 (MQ=255)                                                                            cTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                        >  1:1201989/1‑151 (MQ=255)                                                                            cTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                        >  2:722322/1‑151 (MQ=255)                                                                             taaTCAGGCCGTGCCACTAACGGTC‑CTACATCAATGTGATTCTTTTTCCCCAGCATTTCAAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                              <  2:1657036/142‑1 (MQ=255)                                                                               aTCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGGGGTGTTTTCGCCAGCGTTCCCAAATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGa                                                                     <  2:1666899/151‑1 (MQ=255)                                                                                  aggCCGGGCCACTAACGGTG‑CTACATCAATGTGATGTTTTTTCGCCAGCGGTTCCAGATCAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAg                                                                         <  1:983017/142‑1 (MQ=255)                                                                                  aacccGTGCCACTAACGGTG‑CTACATCAATGTGATGCTTTTCCGCCAGCGGTTCCAGACGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                  <  1:1147227/148‑1 (MQ=255)                                                                                  aaGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                  >  2:866362/1‑151 (MQ=255)                                                                                     ccGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                               >  2:757788/1‑151 (MQ=255)                                                                                             cTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGtt                                                       >  2:1587325/1‑151 (MQ=255)                                                                                             cTAACGGTG‑CTACATAAATGTGAGGCTTTTCCGCCAGCATTCCCAGAGAAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                        <  2:516268/150‑1 (MQ=255)                                                                                              tAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                      >  2:181493/1‑151 (MQ=255)                                                                                                 cGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAAc                                                   >  1:77178/1‑151 (MQ=255)                                                                                                  ggTG‑CTACATCAATGTGATGTGTTTCCGCCAGCAGTCCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCg                                                         <  2:1210852/144‑1 (MQ=255)                                                                                                       cTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                               >  1:1299806/1‑150 (MQ=255)                                                                                                        tACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                              >  1:1480539/1‑150 (MQ=255)                                                                                                         aCATCAATGTGATGCGTTTTCGCCAGCGTTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                                <  2:405364/147‑1 (MQ=255)                                                                                                         aCATAAATGTGATGCGTTTTCCCCAGCAGTTCCAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaa                                             <  2:332293/150‑1 (MQ=255)                                                                                                          cATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                           >  1:660899/1‑151 (MQ=255)                                                                                                           ataaaTGTGTTGCTTTTTCGCCAGCATTTCCAGAGGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                               <  1:1551841/143‑1 (MQ=255)                                                                                                           ataaaTGTGATGCTTTTTCGCCAGCAGTCCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                           <  1:678676/147‑1 (MQ=255)                                                                                                           ataaaTGTGATGCTTTTCCGCCAGCAGTCCAAGAGGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                           <  2:779113/147‑1 (MQ=255)                                                                                                             cAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                        >  1:919694/1‑151 (MQ=255)                                                                                                              aaggtgATGCTTTTCCGCCAGCATTTCCAGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                        <  2:1488682/147‑1 (MQ=255)                                                                                                              aaTGGGTTGCTTTTTCCCCGGCAGTACAAGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                           <  1:1448844/147‑1 (MQ=255)                                                                                                                      gCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAt                               >  2:1649314/1‑151 (MQ=255)                                                                                                                      gCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAt                               >  1:365288/1‑151 (MQ=255)                                                                                                                       cGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                              >  2:1571281/1‑151 (MQ=255)                                                                                                                        ttttcCGCCAGCAGTCCCAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGa                                <  2:441331/147‑1 (MQ=255)                                                                                                                          tttCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                           >  1:498862/1‑151 (MQ=255)                                                                                                                          tttCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACACCGAAACCACCGCTTCGATTTCg                           >  1:91354/1‑151 (MQ=255)                                                                                                                          tttCGCCAGCAGTCCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                            <  2:1244151/150‑1 (MQ=255)                                                                                                                              gCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                       >  1:1645616/1‑151 (MQ=255)                                                                                                                                cagcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                     >  1:1598936/1‑151 (MQ=255)                                                                                                                                 agcagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTgg                    >  2:1249549/1‑151 (MQ=255)                                                                                                                                   cagTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcg                  >  2:536604/1‑151 (MQ=255)                                                                                                                                     tttCCAGAGGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTgg                    <  2:179097/146‑1 (MQ=255)                                                                                                                                      ttCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg               >  1:627976/1‑151 (MQ=255)                                                                                                                                         cAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg             >  1:1555944/1‑150 (MQ=255)                                                                                                                                            atcaaaccGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg         <  2:1311320/144‑1 (MQ=255)                                                                                                                                            aTGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAAt          <  2:1299369/150‑1 (MQ=255)                                                                                                                                               aaaaCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCa      <  2:1075194/150‑1 (MQ=255)                                                                                                                                                aaccGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCaa     <  2:82179/148‑1 (MQ=255)                                                                                                                                                   cGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATgg  >  2:950336/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGG  >  NC_000913/4332643‑4332938

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