Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,850,007	0	G	A	43.4%	88.9 / 42.8	76	P44S (CCG→TCG)	uhpB	sensory histidine kinase in two‑component regulatory sytem with UhpA
Reads supporting (aligned to +/- strand):  ref base G (39/4);  new base A (0/33);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.84e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACAT  >  NC_000913/3849858‑3850153                                                                                                                                                       |                                                                                                                                                    cGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACgg                                                                                                                                                     >  2:1448798/1‑150 (MQ=255)    cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAAt                                                                                                                                                 >  1:714333/1‑151 (MQ=255)       tAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCCGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAAc                                                                                                                                              >  1:1222976/1‑151 (MQ=255)        aaaCTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAAc                                                                                                                                              >  2:1152117/1‑150 (MQ=255)          aCTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑Gc                                                                                                                                          <  2:1443568/151‑1 (MQ=255)             aCCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCAcc                                                                                                                                       >  1:804070/1‑151 (MQ=255)               cGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGc                                                                                                                                     >  2:1441289/1‑151 (MQ=255)                gATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGcc                                                                                                                                    >  1:96502/1‑151 (MQ=255)                gATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGcc                                                                                                                                    >  2:386152/1‑151 (MQ=255)                gATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGcc                                                                                                                                    <  1:1072257/151‑1 (MQ=255)                   catcaATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCat                                                                                                                                  >  2:1085768/1‑150 (MQ=255)                      caATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATc                                                                                                                               >  1:560375/1‑150 (MQ=255)                            cGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGgcgc                                                                                                                        >  2:1374521/1‑151 (MQ=255)                               naaaTGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTc                                                                                                                      >  2:1377594/2‑150 (MQ=255)                               aaaaTGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTc                                                                                                                      >  1:272174/1‑150 (MQ=255)                                 aaTGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAc                                                                                                                   >  1:1071286/1‑151 (MQ=255)                                  agtggTTAAACCCGCCCGCCTGCTTTAGCAATAAAACCAGCGGCCCCCCCGCCCCCAGCAATCCGGGCCAGTTTCCGCGCGGGCATTGCAGCTTTAGCCCCAGCCGCAGACAAAGCGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAAc                                                                                                                   <  1:1212839/148‑1 (MQ=255)                                   tGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAcc                                                                                                                  >  1:756401/1‑150 (MQ=255)                                    gggTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAcc                                                                                                                  >  1:1001019/1‑149 (MQ=255)                                      gTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAcc                                                                                                                  >  1:1333736/1‑147 (MQ=255)                                        tAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            >  2:136404/1‑151 (MQ=255)                                        tAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            >  1:892956/1‑151 (MQ=255)                                             ccgcccgcCGGCGTTAGCAAGAAAATCAGCGCCCCCCCCGCGCCCAGCAATACGGGCCAGTCCCCGCGGGGGCATTGCAGCATTAGCCCCAGACGCAGACAAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGc                                                                                                        <  1:1127072/151‑1 (MQ=255)                                               gaccgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTg                                                                                                      >  2:734071/1‑151 (MQ=255)                                               gaccgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTg                                                                                                      >  2:832785/1‑151 (MQ=255)                                               ccccgCCTGCTTTGCCAAGAAAATCACCAGCCCCCCCGCCCCCAGAAATCCGGGCCAGTACCCGCGCGGGCATTGCAGCTTTAGCCCCAGCCGCAGACCAAACGAAAATAAAA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTg                                                                                                      <  1:105495/149‑1 (MQ=255)                                                       gctttggCCTGTAATTCACCTGCCACCCCGCGCCAAGCAAAACGGCCCGATACCCGCGCGGCGTTTGCAGCATTAGCCCCAGACGCAGACCAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTcc                                                                                               <  1:820491/144‑1 (MQ=255)                                                            aGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGc                                                                                         >  2:1233115/1‑151 (MQ=255)                                                              ccataaaaTCAGCACCCACCCCGCGCCCAGAAATACGGGCCAGTACCCGCGCGGGCATTGCAGCATTACCCCCGGACGCAGACAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCaa                                                                                       <  1:1423669/146‑1 (MQ=255)                                                                      acagcccccACCCCCCCCCCAAAAATACGGGCCAGTACCCGCGCGGGCATTGCAGTATAGGCCCCAGCCACAAACCAACCGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGc                                                                               <  1:63220/144‑1 (MQ=255)                                                                       cagcagCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCg                                                                              >  2:703667/1‑151 (MQ=255)                                                                          cagcCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                            >  2:909736/1‑150 (MQ=255)                                                                           gcccACCCCGCCCCCAGCAACGCGGGCCATTTTCCGCGCGGGAATTCCAGCATTACCCCCAGACCCAGACCAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCag                                                                          <  1:1334379/149‑1 (MQ=255)                                                                               aCTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                       >  2:1155915/1‑150 (MQ=255)                                                                                ccccGCGCCCAGAAATACGGGCCAGTACCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                                                                     <  1:1068593/149‑1 (MQ=255)                                                                                ccccGCCCCCACAAATCCGGGCCAGTCCCCGGGGGGGCTTTGCCACTTTGCCCCCAAACGCAGCCAAAAGGAAAATACCA‑CAACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                                                                     <  1:41273/149‑1 (MQ=255)                                                                                cTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                       >  2:880868/1‑149 (MQ=255)                                                                                 cccGCCCCCAGCAATGCGGGCCAGTATCCGCGGGGGCATTGCAGCATTAGCCCCAGCCGAAGACAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaaga                                                                    <  1:1458248/150‑1 (MQ=255)                                                                                     cgcCCAGAAATCCGGCCCATTCTCCGCGGGGGCATGGCAGCATTAGCCCCAGACGCAGACAAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaa                                                                <  1:1413509/151‑1 (MQ=255)                                                                                      ccccaaaacagcGGCCCAGAACCCGCGGGGCCTTGGCAGCATTACCCCCAGACGCGGACAAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaa                                                               <  2:487007/140‑1 (MQ=255)                                                                                            cAATCCGGGCCAGTACCCGCGGGGGCATGCCAGCATTCCCCCCAGCCGCAGACAAAACGAAAAAAACA‑CCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAg                                                         <  2:264847/151‑1 (MQ=255)                                                                                             aaTACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAgg                                                        >  2:808720/1‑151 (MQ=255)                                                                                              agaCGGGCCAGTACCGGCGCGGGAATTGCAGCATAAGCCCCAGACCCAGACAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGc                                                       <  2:316994/149‑1 (MQ=255)                                                                                              aTACGGGCAAGTACCGGGGGGGGCATGGCAGCATTAGCCCCAGACGCAGCCCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGc                                                       <  2:1351605/151‑1 (MQ=255)                                                                                                aCGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCa                                                      >  2:1309864/1‑150 (MQ=255)                                                                                                  gggCCAGTACCCGCGCGGGCATTGCAGCATTACCCCCAGACGCAGACAAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATa                                                   <  1:158605/151‑1 (MQ=255)                                                                                                       aGTACCCGCGGGGGCTTTGCAGCATTACCCCCAGACGCAAACAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGt                                              <  1:278297/151‑1 (MQ=255)                                                                                                          aTCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCCGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAAt                                           >  1:950519/1‑151 (MQ=255)                                                                                                          aTCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAAt                                           >  1:1109467/1‑151 (MQ=255)                                                                                                             ggcgcgGGCATTGCACCATTCCCCCCAGACGCAAACCAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTaa                                        <  2:113402/150‑1 (MQ=255)                                                                                                             cgcgggGGCTGTGCAGCGTAAGCCCCAGACGAAGCCAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTaa                                        <  1:1226167/151‑1 (MQ=255)                                                                                                              ccggggccATTGCACCTTTCCCCCCACCCCCAGCCCAAGGGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAg                                       <  1:832789/144‑1 (MQ=255)                                                                                                               cgcgGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGc                                      >  1:261736/1‑151 (MQ=255)                                                                                                                gcgGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCg                                     >  2:171844/1‑151 (MQ=255)                                                                                                                  gggCATTGCAGCATTACCCCCAGACGCAGACCAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                   <  1:296188/151‑1 (MQ=255)                                                                                                                     cATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                   >  2:1407872/1‑148 (MQ=255)                                                                                                                     cATTGCAGCATTACCCCCAGCCGCAGACCAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgaga                                <  1:378125/151‑1 (MQ=255)                                                                                                                     cATTGCACCATTACCCCCAGCCGCAGACAAACCGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgaga                                <  1:157883/151‑1 (MQ=255)                                                                                                                        tGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                              >  2:226037/1‑150 (MQ=255)                                                                                                                         ggaacaTCAGCCCCTGCCCCAAACCAAAGAAAAATACCG‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                            <  2:135529/147‑1 (MQ=255)                                                                                                                         gaagaaTTACCCCCAGCCGCAGACCAAACGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaaca                             <  1:998210/145‑1 (MQ=255)                                                                                                                          cagcaTTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                           >  1:174663/1‑151 (MQ=255)                                                                                                                            gaatgcgCCCCAGACGCAGACAAAGCGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGt                         <  2:1222983/145‑1 (MQ=255)                                                                                                                             caTTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTc                        >  1:1157954/1‑151 (MQ=255)                                                                                                                                cacccccAGACGAAAACAAAACGAAAAAAACC‑CCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtc                     <  2:1290211/148‑1 (MQ=255)                                                                                                                                aacccccAGACGCAGCCAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtc                     <  2:1117556/148‑1 (MQ=255)                                                                                                                                  cccccAGACGCAGACCAAAGGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtcat                   <  1:1441286/150‑1 (MQ=255)                                                                                                                                    cccAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtcatc                  >  2:571289/1‑150 (MQ=255)                                                                                                                                     ccaaacacaacccAAACGAAAAAAACA‑CCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAc                <  2:104906/140‑1 (MQ=255)                                                                                                                                      cAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAcc               >  2:993188/1‑151 (MQ=255)                                                                                                                                       agccGCAAACAAAAGGAAAAAACCA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCa              <  1:1448795/148‑1 (MQ=255)                                                                                                                                       aGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAcc               >  1:1016448/1‑150 (MQ=255)                                                                                                                                       aGACGCAGACAAAACGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCa              <  1:54147/151‑1 (MQ=255)                                                                                                                                                caaaaCGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCaaa     <  2:714330/149‑1 (MQ=255)                                                                                                                                                aacaaCCGAAATAACAGGCACCGCCATAACAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGACCAACGTCTTCATCACCAGCCATCat      <  2:96504/148‑2 (MQ=255)                                                                                                                                                  aaaCGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACa   <  1:1152113/151‑1 (MQ=255)                                                                                                                                                   aaCGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACAt  >  2:597201/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    CGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACAT  >  NC_000913/3849858‑3850153

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