Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,0070GA50.5% 79.5 / 163.7 107P44S (CCG→TCG) uhpBsensory histidine kinase in two‑component regulatory sytem with UhpA
Reads supporting (aligned to +/- strand):  ref base G (46/7);  new base A (0/54);  total (46/61)
Fisher's exact test for biased strand distribution p-value = 3.51e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-17

Only 100 of 108 total aligned reads displayed.

CGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACAT  >  NC_000913/3849858‑3850153
                                                                                                                                                        |                                                                                                                                                   
cGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑Cgg                                                                                                                                                     >  4:1448798/1‑150 (MQ=255)
cGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGa                                                                                                                                                    >  3:389571/1‑151 (MQ=255)
 gTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGaa                                                                                                                                                   >  1:483034/1‑151 (MQ=255)
   cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAAt                                                                                                                                                 >  2:714333/1‑151 (MQ=255)
      tAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAAt                                                                                                                                                 <  3:444500/148‑1 (MQ=255)
      tAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAAc                                                                                                                                              <  3:288338/151‑1 (MQ=255)
       aaaCTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAAc                                                                                                                                              >  4:1152117/1‑150 (MQ=255)
         aCTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑Gc                                                                                                                                          <  4:1443568/151‑1 (MQ=255)
            aCCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCAcc                                                                                                                                       >  2:804070/1‑151 (MQ=255)
              cGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGc                                                                                                                                     >  4:1441289/1‑151 (MQ=255)
               gATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGcc                                                                                                                                    >  2:96502/1‑151 (MQ=255)
               gATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGcc                                                                                                                                    >  4:386152/1‑151 (MQ=255)
               gATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGcc                                                                                                                                    <  2:1072257/151‑1 (MQ=255)
                  catcaATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCat                                                                                                                                  >  4:1085768/1‑150 (MQ=255)
                    aaaaaaaacgAAAATGGGGTAAACCCGCCCCCCTGCTTTAGCAAAAAAACCACCACCCCCCCCGCCCCAAGAAATACGGGCCAGTACCCGCGCGGGC‑ATGCCACCATTCCCCCCAGCCGAAGACCAAA‑CGAAAATAACA‑GCACCGCCatat                                                                                                                                <  3:93706/142‑1 (MQ=255)
                     caATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATc                                                                                                                               >  2:560375/1‑150 (MQ=255)
                           cGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGgcgc                                                                                                                        >  4:1374521/1‑151 (MQ=255)
                              naaaTGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTc                                                                                                                      >  4:1377594/2‑150 (MQ=255)
                              aaaaTGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTc                                                                                                                      >  2:272174/1‑150 (MQ=255)
                                aaTGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAc                                                                                                                   >  2:1071286/1‑151 (MQ=255)
                                 agtggTTAAACCCGCCCGCCTGCTTTAGCAATAAAACCAGCGGCCCCCCCGCCCCCAGCAATCCGGGCCAGTTTCCGCGCGGGC‑ATTGCAGCTTTAGCCCCAGCCGCAGACAAAG‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAAc                                                                                                                   <  2:1212839/148‑1 (MQ=255)
                                  tGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAcc                                                                                                                  >  2:756401/1‑150 (MQ=255)
                                   gggTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAcc                                                                                                                  >  2:1001019/1‑149 (MQ=255)
                                   gggTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAg                                                                                                                >  1:451351/1‑151 (MQ=255)
                                     gTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAAcc                                                                                                                  >  2:1333736/1‑147 (MQ=255)
                                       tAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            >  4:136404/1‑151 (MQ=255)
                                       tAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            >  2:892956/1‑151 (MQ=255)
                                        aacccgcccccccTGGGTAACCCAGAAAAACAGCCCCCCCCCCGCCCCCAACAAAACGGGCCAGTACCCGCGCGGGC‑ATTGCAGCTTAACCCCCAGCCGAAAACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            <  3:55195/140‑1 (MQ=255)
                                        aaacccgcccgcCGGCGTTGCCAAGAAAACCAGCACCCCCCCCCCCCCCAAAAATCCGGCCCAGTCCCCGCGGGGGC‑ATGGCAGCATTACCCCCAGACGCAGACAAAA‑GGAAAAAAACA‑CCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            <  3:166185/147‑1 (MQ=255)
                                        aaaccccccccccTGCTTTACCCAGAAAACCACCACCCACCCCCCCCCCACAAAGCCGGGCCAGTACCCGGGCGGCC‑TTTGCACCATACGCCCCAGACGCAGCCAAAA‑GGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            <  3:472889/140‑1 (MQ=255)
                                        aaacccacccgCCTGGGTAGGCAAGAAAACCAGCCCCCCCCCCGCCCCCAACAAGCGGGGCCAGTCCCCGCGCGGGC‑ATGGCAGCATTACCCCCAGACGCAGCCAAAA‑CGAAAATAACA‑GCCCCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                            <  3:414086/143‑1 (MQ=255)
                                            ccgcccgcCGGCGTTAGCAAGAAAATCAGCGCCCCCCCCGCGCCCAGCAATACGGGCCAGTCCCCGCGGGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACAAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGc                                                                                                        <  2:1127072/151‑1 (MQ=255)
                                             cgcccgCCTGCGTTAGCCAGAAAACCGGCACCCCCCCCGCGCCCAACAAAGCGGGCCATTACCCGCGCGGGC‑ATTGCAGCATTACCCCCAGACGCAGACCAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCt                                                                                                       <  3:363338/148‑1 (MQ=255)
                                              gaccgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTg                                                                                                      >  4:734071/1‑151 (MQ=255)
                                              gaccgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTg                                                                                                      >  4:832785/1‑151 (MQ=255)
                                              ccccgCCTGCTTTGCCAAGAAAATCACCAGCCCCCCCGCCCCCAGAAATCCGGGCCAGTACCCGCGCGGGC‑ATTGCAGCTTTAGCCCCAGCCGCAGACCAAA‑CGAAAATAAAA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTg                                                                                                      <  2:105495/149‑1 (MQ=255)
                                                      gctttggCCTGTAATTCACCTGCCACCCCGCGCCAAGCAAAACGGCCCGATACCCGCGCGGCG‑TTTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTcc                                                                                               <  2:820491/144‑1 (MQ=255)
                                                       cGTTAGCCAGAAAACCAGCACCCCCCCCGCCCCCAGAAATACGGGCCAGTACCCGCGGGGGC‑ATTGCAGCATTAGCCCCAGACGAAGACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCcac                                                                                             <  1:389576/151‑1 (MQ=255)
                                                          tAGCAAGAAAACCAGCACCCCCCCCGCCCCCAGAAATACGGGCCAGTACCCGCGGGGGC‑ATTGCAGCATTACCCCCAGACGCAGACCAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAgg                                                                                          <  1:220630/151‑1 (MQ=255)
                                                           aGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGc                                                                                         >  4:1233115/1‑151 (MQ=255)
                                                             ccataaaaTCAGCACCCACCCCGCGCCCAGAAATACGGGCCAGTACCCGCGCGGGC‑ATTGCAGCATTACCCCCGGACGCAGACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCaa                                                                                       <  2:1423669/146‑1 (MQ=255)
                                                                  aaaCCAGCACCCCCCCCGCGCCCAGAAATACGGGCAAGTCCCCGCGGGGGC‑TTGGCAGCATTAGCCCCAGCCGCAACCAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCa                                                                                  <  1:219263/151‑1 (MQ=255)
                                                                     acagcccccACCCCCCCCCCAAAAATACGGGCCAGTACCCGCGCGGGC‑ATTGCAGTATAGGCCCCAGCCACAAACCAAC‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGc                                                                               <  2:63220/144‑1 (MQ=255)
                                                                      cagcagCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCg                                                                              >  4:703667/1‑151 (MQ=255)
                                                                         ccccccccccGCGCCCAGAAATCCGGGCCATTACCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                            <  3:237516/142‑1 (MQ=255)
                                                                         cagcCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                            >  4:909736/1‑150 (MQ=255)
                                                                          gcccACCCCGCCCCCAGCAACGCGGGCCATTTTCCGCGCGGGA‑ATTCCAGCATTACCCCCAGACCCAGACCAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCag                                                                          <  2:1334379/149‑1 (MQ=255)
                                                                              aCTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                       >  4:1155915/1‑150 (MQ=255)
                                                                               ccccGCGCCCAGAAATACGGGCCAGTACCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                                                                     <  2:1068593/149‑1 (MQ=255)
                                                                               ccccGCCCCCACAAATCCGGGCCAGTCCCCGGGGGGGC‑TTTGCCACTTTGCCCCCAAACGCAGCCAAAA‑GGAAAATACCA‑CAACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                                                                     <  2:41273/149‑1 (MQ=255)
                                                                               cTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                       >  4:880868/1‑149 (MQ=255)
                                                                                cccGCCCCCAGCAATGCGGGCCAGTATCCGCGGGGGC‑ATTGCAGCATTAGCCCCAGCCGAAGACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaaga                                                                    <  2:1458248/150‑1 (MQ=255)
                                                                                    cgcCCAGAAATCCGGCCCATTCTCCGCGGGGGC‑ATGGCAGCATTAGCCCCAGACGCAGACAAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaa                                                                <  2:1413509/151‑1 (MQ=255)
                                                                                     ccccaaaacagcGGCCCAGAACCCGCGGGGCC‑TTGGCAGCATTACCCCCAGACGCGGACAAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaa                                                               <  4:487007/140‑1 (MQ=255)
                                                                                        cgacaatgggggCCAGTACCCGCGCGGGC‑ATGGCATCTTTCCCCCCAGCCGAAGACAAAA‑CGAAAAAACCA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAg                                                            <  1:509966/142‑1 (MQ=255)
                                                                                           cAATCCGGGCCAGTACCCGCGGGGGC‑ATGCCAGCATTCCCCCCAGCCGCAGACAAAA‑CGAAAAAAACA‑CCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAg                                                         <  4:264847/151‑1 (MQ=255)
                                                                                            aaTACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAgg                                                        >  4:808720/1‑151 (MQ=255)
                                                                                             aTACGGGCAAGTACCGGGGGGGGC‑ATGGCAGCATTAGCCCCAGACGCAGCCCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGc                                                       <  4:1351605/151‑1 (MQ=255)
                                                                                               aCGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCa                                                      >  4:1309864/1‑150 (MQ=255)
                                                                                                 gggCCAGTACCCGCGCGGGC‑ATTGCAGCATTACCCCCAGACGCAGACAAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATa                                                   <  2:158605/151‑1 (MQ=255)
                                                                                                  ggCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAgg                                                                                          >  1:275706/1‑111 (MQ=255)
                                                                                                  ggCCAGTACCCGCGCGGGC‑ATGGCACCATTACCCCCAGCCGCAGACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAgg                                                                                          <  3:275704/111‑1 (MQ=255)
                                                                                                    ccAGTACCCGCGCGGGC‑ATTGCAGCATTACCCCCAGACGCAGACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAAc                                                 <  1:36827/150‑1 (MQ=255)
                                                                                                    ccAGTACCCGCGCGGGC‑ATGCCAGCATCACCCCCAGCCGAAGACCAAA‑CGAAAAAAACA‑CCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAAc                                                 <  3:9902/150‑1 (MQ=255)
                                                                                                       gTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTa                                             >  1:431510/1‑151 (MQ=255)
                                                                                                         aTCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAAt                                           >  2:1109467/1‑151 (MQ=255)
                                                                                                           ccGCGCGGGC‑ATTGCAGCATTACCCCCAGACGCAGACCAAA‑GGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTa                                         <  1:356621/151‑1 (MQ=255)
                                                                                                            ggcgcgGGC‑ATTGCACCATTCCCCCCAGACGCAAACCAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTaa                                        <  4:113402/150‑1 (MQ=255)
                                                                                                            cgcgggGGC‑TGTGCAGCGTAAGCCCCAGACGAAGCCAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTaa                                        <  2:1226167/151‑1 (MQ=255)
                                                                                                             ccggggcc‑ATTGCACCTTTCCCCCCACCCCCAGCCCAAG‑GGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAg                                       <  2:832789/144‑1 (MQ=255)
                                                                                                              cgcgGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGc                                      >  2:261736/1‑151 (MQ=255)
                                                                                                               gcgGGCAATTGCAGCATGACCCCCAGACGCGGCAAAAA‑CAAAAATAACA‑ACACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTa                                         <  3:530935/148‑1 (MQ=255)
                                                                                                               gcgGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCg                                     >  4:171844/1‑151 (MQ=255)
                                                                                                                 gggC‑ATTGCAGCATTACCCCCAGACGCAGACCAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                   <  2:296188/151‑1 (MQ=255)
                                                                                                                    c‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                   >  4:1407872/1‑148 (MQ=255)
                                                                                                                    c‑ATTGCAGCATTACCCCCAGCCGCAGACCAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgaga                                <  2:378125/151‑1 (MQ=255)
                                                                                                                    c‑ATTGCACCATTACCCCCAGCCGCAGACAAAC‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgaga                                <  2:157883/151‑1 (MQ=255)
                                                                                                                       ttGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaa                               >  1:517471/1‑150 (MQ=255)
                                                                                                                       tggcagcaTTAGCCCCAGACGCA‑ACCAAAGCGGAAAAAACA‑GCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                              <  1:223093/149‑1 (MQ=255)
                                                                                                                        tGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                              >  4:226037/1‑150 (MQ=255)
                                                                                                                         ggaacaTCAGCCCCTGCCCCAAACCAAA‑GAAAAATACCG‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                            <  4:135529/147‑1 (MQ=255)
                                                                                                                         gaagaaTTACCCCCAGCCGCAGACCAAA‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaaca                             <  2:998210/145‑1 (MQ=255)
                                                                                                                         ccagcaTCACCCCCAGCCGAAGACCAAA‑GGAAAAAAACA‑GCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                            <  3:232994/150‑1 (MQ=255)
                                                                                                                          cagcaTTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                           >  2:174663/1‑151 (MQ=255)
                                                                                                                          cagcaTTACCCCCAGACGCAAACAAAA‑CGAAAAAAACA‑CCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                            <  1:105605/150‑1 (MQ=255)
                                                                                                                            gaatgcgCCCCAGACGCAGACAAAG‑CGAAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGt                         <  4:1222983/145‑1 (MQ=255)
                                                                                                                             caTTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTc                        >  2:1157954/1‑151 (MQ=255)
                                                                                                                                aacccccAGACGCAGCCAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtc                     <  4:1117556/148‑1 (MQ=255)
                                                                                                                                  cccccAGACGCAGACCAAA‑GGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtcat                   <  2:1441286/150‑1 (MQ=255)
                                                                                                                                    cccAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtcatc                  >  4:571289/1‑150 (MQ=255)
                                                                                                                                     ccaaacacaacccAAA‑CGAAAAAAACA‑CCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAc                <  4:104906/140‑1 (MQ=255)
                                                                                                                                       agccGCAAACAAAA‑GGAAAAAACCA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCa              <  2:1448795/148‑1 (MQ=255)
                                                                                                                                       aGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAcc               >  2:1016448/1‑150 (MQ=255)
                                                                                                                                       aGACGCAGACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCa              <  2:54147/151‑1 (MQ=255)
                                                                                                                                              gACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGt                         >  3:323358/1‑133 (MQ=255)
                                                                                                                                              gACAAAA‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGt                         <  1:323360/133‑1 (MQ=255)
                                                                                                                                                caaaa‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCaaa     <  4:714330/149‑1 (MQ=255)
                                                                                                                                                aacaa‑CCGAAATAACAGGCACCGCCATAACAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGACCAACGTCTTCATCACCAGCCATCat      <  4:96504/148‑2 (MQ=255)
                                                                                                                                                  aaa‑CGAAAAAAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACa   <  2:1152113/151‑1 (MQ=255)
                                                                                                                                                   aa‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACAt  >  4:597201/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGC‑ATTGCAGCATTAGCCCCAGACGCAGACCAAA‑CGGAAATAACA‑GCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACAT  >  NC_000913/3849858‑3850153

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