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

GTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAC  >  NC_000913/3849859‑3850151
                                                                                                                                                       |                                                                                                                                                 
gTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑a                                                                                                                                                 >  1:136200/1‑150 (MQ=255)
  cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑aaa                                                                                                                                               <  1:1452848/150‑1 (MQ=255)
  cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAAt                                                                                                                                              >  2:1453262/1‑151 (MQ=255)
  cAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAAt                                                                                                                                              <  1:1329535/151‑1 (MQ=255)
      aaaCTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAAc                                                                                                                                           <  2:874488/150‑1 (MQ=255)
      aaaCTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAAc                                                                                                                                           >  1:248159/1‑150 (MQ=255)
          tACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCAc                                                                                                                                      >  2:1086253/1‑151 (MQ=255)
           aCCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCAcc                                                                                                                                     >  1:369427/1‑151 (MQ=255)
             cGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGc                                                                                                                                   >  1:460495/1‑151 (MQ=255)
               atcatcAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGcc                                                                                                                                  >  1:34888/1‑150 (MQ=255)
                tcatcaATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCat                                                                                                                                >  2:598036/1‑151 (MQ=255)
                 catcaATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCat                                                                                                                                >  1:703125/1‑150 (MQ=255)
                     aaTAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCAACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATc                                                                                                                             >  2:198732/1‑149 (MQ=255)
                       tAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAggg                                                                                                                         >  1:510181/1‑151 (MQ=255)
                              aaaaGGGTTAAACCCGCCCGCC‑TGCGTTAGCAAGAAAACCAGCAGCCCCCCCGCGCCCAGAAATACGGGCCAGTACCCGCGCGGGCA‑TGGCAGCATTAGCCCCAGACGCAGACCAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCa                                                                                                                   <  2:787267/151‑1 (MQ=255)
                              aaaaGGGTTAAACCCGCCCCCC‑TGGGTTACCCAGAAAACCAGCAGCCCCCCCGCCCCCAGAAATACGGGCAAGTACCCGCGCGGGCA‑TTGCAGCATTACCCCCAGCCGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCa                                                                                                                   <  2:181978/151‑1 (MQ=255)
                                aTGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAAcc                                                                                                                >  2:1038934/1‑151 (MQ=255)
                                   ggTTAAA‑CCCCCCCCCTTGCTTTAGCCAGAAAATCAGCAGCCACCCCCCCCCCAGCAATACGGGCAGGAACCCGCG‑GGGCATTTGCAGCATTAGCCCCAGACGCACCCAAAAAGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAg                                                                                                              <  1:645869/151‑1 (MQ=255)
                                    gTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAAcc                                                                                                                >  1:762841/1‑147 (MQ=255)
                                      tAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                          >  1:1188571/1‑151 (MQ=255)
                                       aaa‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                          >  1:485022/1‑150 (MQ=255)
                                               ccgCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGAt                                                                                                  >  1:234563/1‑151 (MQ=255)
                                                 gCC‑TGCGTTAGCAAAAAAACCAGCACCCCCCCCGCGCCCAGAAATCCGGGCCAGTACCCGCGCGGGCA‑TTGCAGCATCACCCCCAGACGCAGACAAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATAc                                                                                                <  1:862844/151‑1 (MQ=255)
                                                        gTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCcac                                                                                           >  2:1213939/1‑150 (MQ=255)
                                                               aGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCaaa                                                                                    >  2:1419643/1‑150 (MQ=255)
                                                                      cagcacCCCCCCCGCGCCCAGAAATGCGGGCAAGTACCCGCGCGGGCA‑TTGCACCATTAGCCCCAGACGCAGCCAAAACGA‑AAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCg                                                                            <  1:597981/151‑1 (MQ=255)
                                                                       agccgcccccccccccccAGCAATCCGGCCCAGTACCCGCGGGGGCA‑TGGCAGCATTACCCCCAGCCGAAAACAAAAGGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCgg                                                                           <  2:879523/136‑1 (MQ=255)
                                                                         cagcCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                          >  1:1055477/1‑150 (MQ=255)
                                                                           gcCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                          >  2:488667/1‑148 (MQ=255)
                                                                             ccccccGCGCCAAGAAATCCGGGCCAGTACCCGCGGGGGCA‑TTGCAGCATTACCCCCAGACGCAGACAAAACGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                     <  2:732032/147‑1 (MQ=255)
                                                                             ccccccGCCCCCAGCAATGCGGGCCAGTACCCGCGCGGGCA‑TTGCAGCATTAGCCCCGGACGCAGACAAAAGGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                     <  2:1308149/147‑1 (MQ=255)
                                                                               ccccGCCCCCAGCAATACGGGCCAGTACCCGCGCGGGCA‑TTGCAGCATTACCCCCAGACGCAAACAAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                                                                   <  2:665853/149‑1 (MQ=255)
                                                                                 ccGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAgaga                                                                     >  1:1311977/1‑147 (MQ=255)
                                                                                     gcCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaa                                                              >  1:715744/1‑150 (MQ=255)
                                                                                     ccccAGAAAACCGGGCCAGTACCGGCGCGGGCA‑TTGCAGCATTACCCCCAGACGCAGACAAAACGA‑AAAAAACGGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaa                                                             <  1:1453266/150‑1 (MQ=255)
                                                                                      cccagaaaccgggcccAGTACCCGCGGGGGAA‑TTGCAGCATTACCCCCAGCCGAAGACAAAACGA‑AAAAAACCGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaa                                                            <  1:1204382/137‑1 (MQ=255)
                                                                                      cccAGAAATCGGGCCCAGTACCCGGGGGGGCA‑TTGCACCATGAGCCCAAGCCGCAGCCAAAAGGA‑AAATAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaa                                                            <  2:136193/151‑1 (MQ=255)
                                                                                       ccgccaaggcGGGCCAGTACCCGCGGGGGCA‑TTGCAGCATAGGCCCGAGACGAAAAACAAGCGA‑AAAAACCGGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaaa                                                           <  1:938557/142‑1 (MQ=255)
                                                                                         aGCAATACGGCCCATTACCCGCGGGGGCA‑TTGCAGCATTACCCCCAGACGCAGCCAAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGc                                                         <  2:507434/151‑1 (MQ=255)
                                                                                           cAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGc                                                         >  1:520503/1‑149 (MQ=255)
                                                                                            aaTACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAgg                                                      >  1:1432974/1‑151 (MQ=255)
                                                                                                 gggCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAAt                                                  >  1:454017/1‑150 (MQ=255)
                                                                                                 gggCCAGTACCCGCGGGGGCA‑TTGCAGCTTTAGCCCAAGACGCAAACCAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATa                                                 <  2:40927/151‑1 (MQ=255)
                                                                                                    ccattccccgggggggCA‑TTGAATCTTTACCCCCAGACGCAGACCAAACGA‑AAAAACCAGCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACg                                              <  1:435485/138‑1 (MQ=255)
                                                                                                      agtccccGCGCGGGCA‑TTGCAGAATTACCCCCAGACGCAGACAAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGt                                            <  1:1291505/146‑1 (MQ=255)
                                                                                                       gTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGt                                            >  1:729431/1‑150 (MQ=255)
                                                                                                         ccccGCGCGGGCT‑TGGCAG‑AATAGCCCCCAGCCCAAGCCAAACGGAAAAAAAAAACACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAAt                                         <  1:635167/149‑1 (MQ=255)
                                                                                                         acccGCGCGGGCA‑TTGCAGCATTACCCCCAGACGCAGACAAAACGA‑AAAAAACAGCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTaa                                          <  1:1177260/148‑1 (MQ=255)
                                                                                                           ccGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACAAAACGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTa                                       <  2:8464/151‑1 (MQ=255)
                                                                                                            cgcgggGGCT‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTaa                                      <  2:369420/151‑1 (MQ=255)
                                                                                                               gcgGGCA‑TTGAAGCATGAGCCCCAGACGAAAACAAAACGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCg                                   <  1:486973/151‑1 (MQ=255)
                                                                                                                 gggCA‑TTGCAGCATTACCCCGAGACGCAGACCAAACGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                 <  2:1156143/151‑1 (MQ=255)
                                                                                                                   gCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                 >  2:734480/1‑149 (MQ=255)
                                                                                                                   gCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCggg                                 >  1:992763/1‑149 (MQ=255)
                                                                                                                    cA‑TGGCAGCATCAGCCCCAGACGCAGACAAAACGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgaga                              <  1:1452181/151‑1 (MQ=255)
                                                                                                                      atagcagatttccccccAGCCGCAGACAAAAGGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaa                             <  2:1220276/138‑1 (MQ=255)
                                                                                                                       gtGCGGCATTACCCCCAGACGCAAACAAAAAGA‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                            <  2:223089/150‑1 (MQ=255)
                                                                                                                        ggcagcaTTACCCCCAGCCGCAGACCAAACGG‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                            <  2:510175/149‑1 (MQ=255)
                                                                                                                        ggcagcaTTACCCCCAGCCGAAGACCAAACGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaac                            <  2:1071053/149‑1 (MQ=255)
                                                                                                                         gcaccaTTACCCCCAGACGCAGACAAAAGGA‑AAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                          <  2:460486/151‑1 (MQ=255)
                                                                                                                          cgcattgccccccAGCCCCACACCAACCAA‑AAAAAACACCTCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                         <  2:931455/142‑1 (MQ=255)
                                                                                                                          cagcaTTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                         >  1:118496/1‑151 (MQ=255)
                                                                                                                          caccttcccccccAGAGCCAAACAAAACAA‑AAAAAACACCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                         <  1:1086252/142‑1 (MQ=255)
                                                                                                                           agcaTTACCCCAAGACGAAAACCAAACGA‑AAATACCAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACg                        <  2:485014/151‑1 (MQ=255)
                                                                                                                              aTTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCt                     >  1:1447383/1‑151 (MQ=255)
                                                                                                                                 acccccCGACGCAGACCAAGGAA‑AAAAAACCTCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtca                  <  1:1076757/149‑1 (MQ=255)
                                                                                                                                 aGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtc                   >  2:1146176/1‑150 (MQ=255)
                                                                                                                                    cccAGCCGAACACAAAACGA‑AAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtcatca               <  1:655628/151‑1 (MQ=255)
                                                                                                                                     ccAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAc              >  1:1014127/1‑151 (MQ=255)
                                                                                                                                      cAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCAcc             >  1:276467/1‑151 (MQ=255)
                                                                                                                                           gCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGcc         >  1:499908/1‑150 (MQ=255)
                                                                                                                                              gACCAAACAA‑AAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCa     <  1:1133293/151‑1 (MQ=255)
                                                                                                                                                  aaaCGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAc  >  2:1118629/1‑150 (MQ=255)
                                                                                                                                                  aaaCGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAc  >  1:1314089/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAA‑CCGACCGCC‑TGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCA‑TTGCAGCATTAGCCCCAGACGCAGACCAAACGG‑AAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAC  >  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: ‑