Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,0070GA43.2% 85.1 / 39.2 74P44S (CCG→TCG) uhpBsensory histidine kinase in two‑component regulatory sytem with UhpA
Reads supporting (aligned to +/- strand):  ref base G (36/6);  new base A (0/32);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 3.15e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAC  >  NC_000913/3849859‑3850151
                                                                                                                                                       |                                                                                                                                                 
gTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGa                                                                                                                                                 >  1:136200/1‑150 (MQ=255)
  cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGaaa                                                                                                                                               <  1:1452848/150‑1 (MQ=255)
  cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAAt                                                                                                                                              >  2:1453262/1‑151 (MQ=255)
  cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAAt                                                                                                                                              <  1:1329535/151‑1 (MQ=255)
      aaaCTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAAc                                                                                                                                           <  2:874488/150‑1 (MQ=255)
      aaaCTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAAc                                                                                                                                           >  1:248159/1‑150 (MQ=255)
          tACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCAc                                                                                                                                      >  2:1086253/1‑151 (MQ=255)
           aCCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCAcc                                                                                                                                     >  1:369427/1‑151 (MQ=255)
             cGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGc                                                                                                                                   >  1:460495/1‑151 (MQ=255)
               atcatcAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGcc                                                                                                                                  >  1:34888/1‑150 (MQ=255)
                tcatcaATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCat                                                                                                                                >  2:598036/1‑151 (MQ=255)
                 catcaATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCat                                                                                                                                >  1:703125/1‑150 (MQ=255)
                     aaTAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCAACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATc                                                                                                                             >  2:198732/1‑149 (MQ=255)
                       tAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAggg                                                                                                                         >  1:510181/1‑151 (MQ=255)
                              aaaaGGGTTAAACCCGCCCGCCTGCGTTAGCAAGAAAACCAGCAGCCCCCCCGCGCCCAGAAATACGGGCCAGTACCCGCGCGGGCATGGCAGCATTAGCCCCAGACGCAGACCAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCa                                                                                                                   <  2:787267/151‑1 (MQ=255)
                              aaaaGGGTTAAACCCGCCCCCCTGGGTTACCCAGAAAACCAGCAGCCCCCCCGCCCCCAGAAATACGGGCAAGTACCCGCGCGGGCATTGCAGCATTACCCCCAGCCGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCa                                                                                                                   <  2:181978/151‑1 (MQ=255)
                                aTGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAAcc                                                                                                                >  2:1038934/1‑151 (MQ=255)
                                   ggTTAAACCCCCCCCCTTGCTTTAGCCAGAAAATCAGCAGCCACCCCCCCCCCAGCAATACGGGCAGGAACCCGCGGGGCATTTGCAGCATTAGCCCCAGACGCACCCAAAAAGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAg                                                                                                              <  1:645869/151‑1 (MQ=255)
                                    gTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAAcc                                                                                                                >  1:762841/1‑147 (MQ=255)
                                      tAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                          >  1:1188571/1‑151 (MQ=255)
                                       aaaCCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                          >  1:485022/1‑150 (MQ=255)
                                               ccgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGAt                                                                                                  >  1:234563/1‑151 (MQ=255)
                                                 gCCTGCGTTAGCAAAAAAACCAGCACCCCCCCCGCGCCCAGAAATCCGGGCCAGTACCCGCGCGGGCATTGCAGCATCACCCCCAGACGCAGACAAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATAc                                                                                                <  1:862844/151‑1 (MQ=255)
                                                        gTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCcac                                                                                           >  2:1213939/1‑150 (MQ=255)
                                                               aGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCaaa                                                                                    >  2:1419643/1‑150 (MQ=255)
                                                                      cagcacCCCCCCCGCGCCCAGAAATGCGGGCAAGTACCCGCGCGGGCATTGCACCATTAGCCCCAGACGCAGCCAAAACGAAAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCg                                                                            <  1:597981/151‑1 (MQ=255)
                                                                       agccgcccccccccccccAGCAATCCGGCCCAGTACCCGCGGGGGCATGGCAGCATTACCCCCAGCCGAAAACAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCgg                                                                           <  2:879523/136‑1 (MQ=255)
                                                                         cagcCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                          >  1:1055477/1‑150 (MQ=255)
                                                                           gcCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                          >  2:488667/1‑148 (MQ=255)
                                                                             ccccccGCGCCAAGAAATCCGGGCCAGTACCCGCGGGGGCATTGCAGCATTACCCCCAGACGCAGACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                     <  2:732032/147‑1 (MQ=255)
                                                                             ccccccGCCCCCAGCAATGCGGGCCAGTACCCGCGCGGGCATTGCAGCATTAGCCCCGGACGCAGACAAAAGGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                     <  2:1308149/147‑1 (MQ=255)
                                                                               ccccGCCCCCAGCAATACGGGCCAGTACCCGCGCGGGCATTGCAGCATTACCCCCAGACGCAAACAAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                                                                   <  2:665853/149‑1 (MQ=255)
                                                                                 ccGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                     >  1:1311977/1‑147 (MQ=255)
                                                                                     gcCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaa                                                              >  1:715744/1‑150 (MQ=255)
                                                                                     ccccAGAAAACCGGGCCAGTACCGGCGCGGGCATTGCAGCATTACCCCCAGACGCAGACAAAACGAAAAAAACGGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaa                                                             <  1:1453266/150‑1 (MQ=255)
                                                                                      cccagaaaccgggcccAGTACCCGCGGGGGAATTGCAGCATTACCCCCAGCCGAAGACAAAACGAAAAAAACCGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaa                                                            <  1:1204382/137‑1 (MQ=255)
                                                                                      cccAGAAATCGGGCCCAGTACCCGGGGGGGCATTGCACCATGAGCCCAAGCCGCAGCCAAAAGGAAAATAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaa                                                            <  2:136193/151‑1 (MQ=255)
                                                                                       ccgccaaggcGGGCCAGTACCCGCGGGGGCATTGCAGCATAGGCCCGAGACGAAAAACAAGCGAAAAAACCGGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaaa                                                           <  1:938557/142‑1 (MQ=255)
                                                                                         aGCAATACGGCCCATTACCCGCGGGGGCATTGCAGCATTACCCCCAGACGCAGCCAAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGc                                                         <  2:507434/151‑1 (MQ=255)
                                                                                           cAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGc                                                         >  1:520503/1‑149 (MQ=255)
                                                                                            aaTACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAgg                                                      >  1:1432974/1‑151 (MQ=255)
                                                                                                 gggCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAAt                                                  >  1:454017/1‑150 (MQ=255)
                                                                                                 gggCCAGTACCCGCGGGGGCATTGCAGCTTTAGCCCAAGACGCAAACCAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATa                                                 <  2:40927/151‑1 (MQ=255)
                                                                                                    ccattccccgggggggCATTGAATCTTTACCCCCAGACGCAGACCAAACGAAAAAACCAGCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACg                                              <  1:435485/138‑1 (MQ=255)
                                                                                                      agtccccGCGCGGGCATTGCAGAATTACCCCCAGACGCAGACAAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGt                                            <  1:1291505/146‑1 (MQ=255)
                                                                                                       gTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGt                                            >  1:729431/1‑150 (MQ=255)
                                                                                                         ccccGCGCGGGCTTGGCAGAATAGCCCCCAGCCCAAGCCAAACGGAAAAAAAAAACACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAAt                                         <  1:635167/149‑1 (MQ=255)
                                                                                                         acccGCGCGGGCATTGCAGCATTACCCCCAGACGCAGACAAAACGAAAAAAACAGCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTaa                                          <  1:1177260/148‑1 (MQ=255)
                                                                                                           ccGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTa                                       <  2:8464/151‑1 (MQ=255)
                                                                                                            cgcgggGGCTTTGCAGCATTAGCCCCAGACGCAGACCAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTaa                                      <  2:369420/151‑1 (MQ=255)
                                                                                                               gcgGGCATTGAAGCATGAGCCCCAGACGAAAACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCg                                   <  1:486973/151‑1 (MQ=255)
                                                                                                                 gggCATTGCAGCATTACCCCGAGACGCAGACCAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                 <  2:1156143/151‑1 (MQ=255)
                                                                                                                   gCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                 >  2:734480/1‑149 (MQ=255)
                                                                                                                   gCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                 >  1:992763/1‑149 (MQ=255)
                                                                                                                    cATGGCAGCATCAGCCCCAGACGCAGACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgaga                              <  1:1452181/151‑1 (MQ=255)
                                                                                                                      atagcagatttccccccAGCCGCAGACAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaa                             <  2:1220276/138‑1 (MQ=255)
                                                                                                                       gtGCGGCATTACCCCCAGACGCAAACAAAAAGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                            <  2:223089/150‑1 (MQ=255)
                                                                                                                        ggcagcaTTACCCCCAGCCGCAGACCAAACGGAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                            <  2:510175/149‑1 (MQ=255)
                                                                                                                        ggcagcaTTACCCCCAGCCGAAGACCAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                            <  2:1071053/149‑1 (MQ=255)
                                                                                                                         gcaccaTTACCCCCAGACGCAGACAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                          <  2:460486/151‑1 (MQ=255)
                                                                                                                          cgcattgccccccAGCCCCACACCAACCAAAAAAAACACCTCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                         <  2:931455/142‑1 (MQ=255)
                                                                                                                          cagcaTTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                         >  1:118496/1‑151 (MQ=255)
                                                                                                                          caccttcccccccAGAGCCAAACAAAACAAAAAAAACACCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                         <  1:1086252/142‑1 (MQ=255)
                                                                                                                           agcaTTACCCCAAGACGAAAACCAAACGAAAATACCAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACg                        <  2:485014/151‑1 (MQ=255)
                                                                                                                              aTTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCt                     >  1:1447383/1‑151 (MQ=255)
                                                                                                                                 acccccCGACGCAGACCAAGGAAAAAAAACCTCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtca                  <  1:1076757/149‑1 (MQ=255)
                                                                                                                                 aGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtc                   >  2:1146176/1‑150 (MQ=255)
                                                                                                                                    cccAGCCGAACACAAAACGAAAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtcatca               <  1:655628/151‑1 (MQ=255)
                                                                                                                                     ccAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAc              >  1:1014127/1‑151 (MQ=255)
                                                                                                                                      cAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAcc             >  1:276467/1‑151 (MQ=255)
                                                                                                                                           gCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGcc         >  1:499908/1‑150 (MQ=255)
                                                                                                                                              gACCAAACAAAAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCa     <  1:1133293/151‑1 (MQ=255)
                                                                                                                                                  aaaCGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAc  >  2:1118629/1‑150 (MQ=255)
                                                                                                                                                  aaaCGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAc  >  1:1314089/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAC  >  NC_000913/3849859‑3850151

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: