Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,7910AC44.8% 181.7 / 33.3 105V161G (GTG→GGG) basSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (54/4);  new base C (0/47);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 1.40e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGC  >  NC_000913/4332646‑4332940
                                                                                                                                                    |                                                                                                                                                     
ttGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCAccccc                                                                                                                                                  >  1:103524/1‑151 (MQ=255)
  gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGc                                                                                                                                                <  1:1188998/151‑1 (MQ=255)
    tAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGcc                                                                                                                                               >  2:1516755/1‑150 (MQ=255)
     aaTTCCCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCa                                                                                                                                              <  2:193339/150‑1 (MQ=255)
     aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGcc                                                                                                                                               >  1:1281325/1‑149 (MQ=255)
        tACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTgg                                                                                                                                          >  1:1348639/1‑151 (MQ=255)
         aCCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGc                                                                                                                                         >  1:344277/1‑151 (MQ=255)
          ccGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCg                                                                                                                                        >  1:270613/1‑151 (MQ=255)
          ccGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCg                                                                                                                                        >  2:1039491/1‑151 (MQ=255)
           cGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGt                                                                                                                                       >  1:101546/1‑151 (MQ=255)
               aaaaaaTGCCTGCCGGCCACGCGCCAGTGGAAGCGGCGGGGAGCCGCTCTCCTTCATCTGACCAGGCCGTCCCACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCA‑GCATTTCCAG‑ATGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGc                                                                                                                                   <  2:396610/149‑1 (MQ=255)
                gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                   >  2:227123/1‑150 (MQ=255)
                  aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGt                                                                                                                                >  2:349524/1‑151 (MQ=255)
                  aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGt                                                                                                                                >  2:1763096/1‑151 (MQ=255)
                   aaTGACTGCCCGGCCCGCGCCGGTTGCAGCGGCGGGGAGCCGCTCTCCATCTTCTGATCAAGCCGTCCCACTAACGGTGCTACATCAATGTGATGCTTTTCCGCCA‑GCATTCCCAG‑AGGCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                               <  1:764504/151‑1 (MQ=255)
                          gTCCGGCCCGCGCTAGTTGCACCGGCTGGGAGCCGTCCTCCTTCATCTGATCAGGCCGTCCCACTAAGGGTGCTACATAAATGTGATGCTTTTTCGCCA‑GCAGTTCCAG‑AAAAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGt                                                                                                                            <  1:1087124/147‑1 (MQ=255)
                          gTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                        >  2:1020382/1‑151 (MQ=255)
                          gTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                        >  2:705871/1‑151 (MQ=255)
                            ccGGCGCGCGCAAGTGGCAGCAGCGGGGAGCCGCTCTCCATCATCTGACCAAGCCGTGCCACTAACGGTGCTACACCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGAAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGa                                                                                                                      <  1:1489501/151‑1 (MQ=255)
                            ccGGCCCGCGCCAGTTGCAGCGGCGGGGAGCCGCTCTCCACCATCTGATCAAGCCGTGCCACTAACGGTGCAACACCAATGTGATGCGTTTCCGCCA‑GCAGTCCCAG‑ATGCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                                          <  1:658508/147‑1 (MQ=255)
                            ccGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGcgcg                                                                                                                       >  1:14591/1‑150 (MQ=255)
                              ggcccgccctatttGCAGCAGCGGGGAACCCCTCTCAACCATCTCATCAACCCGTCCCCCAACCGTTGCTACATCAATGTGTTGTTTTTTCGCCA‑GCATTTCAAA‑ATGCAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGt                                                                                                                            <  2:1013531/131‑1 (MQ=255)
                              ggCACGCGCCAGTTGCAGCGGCGGGGAGACGCTCCCCATCATCTGATCAGGCCGTGCCACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCA‑GCAGTTCCAG‑ATGAAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACg                                                                                                                    <  1:76770/151‑1 (MQ=255)
                               gCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGt                                                                                                                   >  1:434947/1‑151 (MQ=255)
                               gCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGt                                                                                                                   >  1:49452/1‑151 (MQ=255)
                                 aCGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                  >  1:959378/1‑150 (MQ=255)
                                   gcgcCATTTGCAGCGGTGGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAAGGGTGCTACATAAATGTGATGCGTTTTCGCCA‑GCAGTCCCAG‑AGGAAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACg                                                                                                                    <  1:1441563/146‑1 (MQ=255)
                                    cgcCAGTGGCGGCGGGGGGGAGCCCCCCCCCATCTTCTAACCAGGCCGTCCCACTAACGGTGCTACATAATTGTGTGGCTTTTCCCCCG‑GCAGTCCCAG‑ATGAAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGAc                                                                                                                     <  2:1567619/144‑1 (MQ=255)
                                      ccAGTTGAAGCAGCTGGGAGCCCCTCTCCATCTTCTGACCAAGCCGTCCCACTACCGGTGCTACATCAATGTGATGGTTTTCCGCCA‑GCATTTCCAG‑ATGAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                            <  1:1070089/151‑1 (MQ=255)
                                       cAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                            >  2:788794/1‑150 (MQ=255)
                                       cAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                            >  1:1005912/1‑150 (MQ=255)
                                        aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                                           >  2:1446872/1‑150 (MQ=255)
                                            gcagcaGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                       >  1:1075348/1‑150 (MQ=255)
                                              agcagcGGGGAGCCCCCCTCCATCATTTGACCAAGCCGTGCCACTAACGGTGCTACATAAATGTGTTGCGTTTTCGCCA‑GCAGTCCCAG‑ATGCAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                        <  1:896545/147‑1 (MQ=255)
                                               gcagcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCt                                                                                                   >  1:1352715/1‑151 (MQ=255)
                                                  gcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                        <  1:865859/143‑1 (MQ=255)
                                                   cTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCg                                                                                               >  1:1453160/1‑151 (MQ=255)
                                                    ggggAGCCCCTCTCCACCTTCTGCCCACGCCGTCCCACTAACGGTTCTACATAAATGTGATGCTTTTCCGCCA‑GCAGTCCCAG‑ATGAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                              <  2:499711/150‑1 (MQ=255)
                                                    ggggAGCCCCTCCCCATCATCTGACCAGGCCGTCCCACTAGCGGTGCTACACCAATGTGATGCTTTTCCGCCA‑GCGGTTCCAG‑ATGAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGt                                                                                              <  1:1080599/150‑1 (MQ=255)
                                                     gggAGCCGCTCCCCTTCATCTGATCAAGCCGTCCCAAAAACGGTGCTACATAAATGTGATGCGTTTCCGCCA‑GCGGTTCCAG‑ATGCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGTGACGTCAGCGGTAAACAACCttt                                                                                                 <  1:638114/147‑1 (MQ=255)
                                                     gggAGCCCCTCCCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATAAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCg                                                                                               <  1:766833/149‑1 (MQ=255)
                                                      ggAGCCGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCTTTTCCGCCA‑GCATTTCCAG‑ATGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTa                                                                                            <  1:1158063/151‑1 (MQ=255)
                                                        agccGCTCCCCATCATCTAATAACGCGGCGCCACTAACGGTGCTAAATAAATGTGTTGCTTTTCCCCCG‑GCATTTCCAG‑ATGAAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATc                                                                                          <  1:636554/148‑1 (MQ=255)
                                                          ccGCTCTCCATCATCTAATCAAGCCGGCCCACTAACGGTGCTACATCAATGGGATGCTTTTCCGCCA‑GCAGTTCCAG‑ACACAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATcc                                                                                         <  2:1077879/149‑1 (MQ=255)
                                                           cGCTCCCCATCATCTGATCAGGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCATTTCCAG‑ATGAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                       <  1:1524275/151‑1 (MQ=255)
                                                            gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                      >  2:1575945/1‑151 (MQ=255)
                                                            gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                      >  1:1114806/1‑151 (MQ=255)
                                                            gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                      >  2:270145/1‑151 (MQ=255)
                                                             ctctCCATCATCTGATCACGCCGTGCCACTAACGGTGCTACATAAATGTGTTGCGTTTCCGCCA‑GCAGTTCCAG‑ATGCAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                      <  2:908795/150‑1 (MQ=255)
                                                              tctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGt                                                                                    >  2:1368649/1‑151 (MQ=255)
                                                               ctcCATCTTCTGATCAAGCCGTGCCACTACCGGTGCTACATAAATGTGATGCTTTTTCCCCA‑GCAGTTCAAG‑ATGAAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATc                                                                                          <  1:1438391/144‑1 (MQ=255)
                                                                tcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGc                                                                                  >  2:1042981/1‑151 (MQ=255)
                                                                   atcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTCCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                                <  1:972430/150‑1 (MQ=255)
                                                                          gATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCTTTTCCGCCA‑GCAGTTCCAG‑AGGAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAg                                                                           <  1:1708902/148‑1 (MQ=255)
                                                                          gATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                        >  2:1514159/1‑151 (MQ=255)
                                                                          aaccagcccGTGCCACTAGGGGTCCAACACAAAGGTGTTGCGTTTCCGCCG‑GCAGTTCCGG‑AGCCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGc                                                                                  <  2:1275429/134‑1 (MQ=255)
                                                                           accaggCCGTGCCACTAACGGTGTTACAACAAGGGGATCGGTTTCCGCCA‑GCTGTCCCGG‑AGCCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                                <  2:1549527/137‑1 (MQ=255)
                                                                            tCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                      >  2:405265/1‑151 (MQ=255)
                                                                            tCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                      >  2:237285/1‑151 (MQ=255)
                                                                             caggcggccccccTAACGGTCCTAATACAATGTGTTTTTTTTGCCCCA‑GCTG‑TCCAGAAGGCAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                         <  2:1472502/135‑1 (MQ=255)
                                                                              aaGCCGGGCCACTAACGGTGCTACATAAATGTGATGCGTTTTCCCCA‑GCAGTCCCAG‑AGGAAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                    <  1:1757060/151‑1 (MQ=255)
                                                                                  cGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                 >  2:141321/1‑150 (MQ=255)
                                                                                     gCCACTAAGGGGGCTACATCAATGGGTTGCGTTTTCGCCA‑GCAGTTCCAG‑ATGAAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGa                                                              <  1:892497/150‑1 (MQ=255)
                                                                                      ccACAACCGGTGCTACATCAAT‑TGATTGCTTTTCCCCAGGCAGTCCCAG‑ATGAAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                                     <  2:1674197/142‑1 (MQ=255)
                                                                                       cACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCA‑GCATTTCCAG‑ATGCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCg                                                           <  1:437537/151‑1 (MQ=255)
                                                                                       cACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCA‑GCAGTTCCAG‑ATGCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                               <  2:755796/147‑1 (MQ=255)
                                                                                           aagggtgtaaCATCAATGTGTTGCTTTTTCGCCA‑GCATTTCCAG‑ATGCAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTa                                                       <  1:1039544/142‑1 (MQ=255)
                                                                                           aaCGGTGATACCACAATGTGATGTTTTTTCCCCA‑GCAGTTCAAA‑AAACAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                            <  1:1399892/146‑1 (MQ=255)
                                                                                               gTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                                   >  2:81575/1‑151 (MQ=255)
                                                                                               gTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                                   >  1:702574/1‑151 (MQ=255)
                                                                                                  cTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                 >  1:198229/1‑150 (MQ=255)
                                                                                                   tACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                                >  1:895533/1‑150 (MQ=255)
                                                                                                     cATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                             >  1:871158/1‑151 (MQ=255)
                                                                                                       tCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                            >  2:1629654/1‑150 (MQ=255)
                                                                                                        cAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                          >  2:581565/1‑151 (MQ=255)
                                                                                                         aaTGTGTTGCTTTTCCGCCA‑GCATTCCAAG‑ATGCAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                         <  2:1621821/151‑1 (MQ=255)
                                                                                                         aaTGTGATGCTTTTTCGCCA‑GCGGTTCAAG‑ATGAAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                          <  1:342908/150‑1 (MQ=255)
                                                                                                              gATGCTTTTCCGCCA‑GCGTTTCAAG‑ATGCAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                                     <  1:1236322/150‑1 (MQ=255)
                                                                                                                tGCGTTTTCGCCA‑GCAGTCCCAG‑ATGCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                                   <  1:2761/150‑1 (MQ=255)
                                                                                                                 gCGTTTCCGCCG‑GCTGTCCGAG‑GTGAAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                                   <  2:1659840/149‑1 (MQ=255)
                                                                                                                   ttttcCGCCA‑GCAGTCCGAG‑ATGCACC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGa                                  <  2:1588/147‑1 (MQ=255)
                                                                                                                   gTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt                               >  1:915212/1‑151 (MQ=255)
                                                                                                                    ttttCGCCA‑GCAGTTCCAG‑ATGCAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                                   <  1:1260997/146‑1 (MQ=255)
                                                                                                                    ttttCCCCA‑GCAGTCCCAG‑ATCAAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                                <  1:1089906/149‑1 (MQ=255)
                                                                                                                     tttCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                             >  1:1068283/1‑151 (MQ=255)
                                                                                                                     ttcccccA‑GCAGTCCCAG‑ATGAAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                              <  2:193396/145‑1 (MQ=255)
                                                                                                                      ttCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                             >  1:590737/1‑150 (MQ=255)
                                                                                                                         gCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                         >  1:1665327/1‑151 (MQ=255)
                                                                                                                           ca‑gcagTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                       >  1:1153051/1‑151 (MQ=255)
                                                                                                                            a‑gcagTTCCAG‑AGGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                        <  1:808963/149‑1 (MQ=255)
                                                                                                                              gcagTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgc                     >  2:10826/1‑151 (MQ=255)
                                                                                                                              gcagTCCAAG‑AGGCAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                        <  2:270600/148‑1 (MQ=255)
                                                                                                                               cagTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcg                    >  1:698992/1‑151 (MQ=255)
                                                                                                                                      caagataaaaccGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt                <  1:1311293/139‑1 (MQ=255)
                                                                                                                                         aTGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg           >  1:225081/1‑151 (MQ=255)
                                                                                                                                          agCAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg           <  2:1313571/150‑1 (MQ=255)
                                                                                                                                             aaa‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCa        >  2:1619374/1‑150 (MQ=255)
                                                                                                                                              aa‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt      >  2:236753/1‑151 (MQ=255)
                                                                                                                                              aa‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt      >  2:861459/1‑151 (MQ=255)
                                                                                                                                                   cACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGc  >  1:40236/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCA‑GCAGTTCCAG‑ATGCAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGC  >  NC_000913/4332646‑4332940

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