Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG40.2% 143.1 / 20.2 87T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (9/43);  new base G (35/0);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 3.38e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCG  >  NC_000913/4346525‑4346823
                                                                                                                                                       |                                                                                                                                                    
ggAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGt                                                                                                                                                      <  2:573106/151‑1 (MQ=255)
   aCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACg                                                                                                                                                   <  2:1224831/151‑1 (MQ=255)
    cGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑AAGGGTGGAGAACCGCGGGGAc                                                                                                                                                    >  1:430679/1‑149 (MQ=255)
        tGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCAt                                                                                                                                                <  1:586193/149‑1 (MQ=255)
         gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCAtttt                                                                                                                                             <  2:1213463/151‑1 (MQ=255)
         gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCAtt                                                                                                                                               <  1:681295/149‑1 (MQ=255)
          aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTtc                                                                                                                                            <  1:823128/151‑1 (MQ=255)
          aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTGGAAGGTTTGCAGGTTGGTCTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTGtc                                                                                                                                            <  1:460875/151‑1 (MQ=255)
            cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                          <  2:1222069/151‑1 (MQ=255)
                  cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑GCGGCTGCCGTACCGGGGGGACGCCTTTTCTccacggg                                                                                                                                    >  1:1223163/1‑149 (MQ=255)
                  cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAg                                                                                                                                    <  2:1190822/151‑1 (MQ=255)
                   gttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGCACCGAGGGGACGCATTTTCTCGACGCgg                                                                                                                                   >  1:292213/1‑151 (MQ=255)
                   gttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGAAGCATTTTTTCGACGAgg                                                                                                                                   >  2:1129455/1‑151 (MQ=255)
                    ttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACCAGGGGACGCATTTTCTCGACGAgg                                                                                                                                   >  2:850538/1‑150 (MQ=255)
                     tcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGa                                                                                                                                  <  1:68127/150‑1 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGCTGCAGAAACGGGGGGACGCATTTTTTCGACGAGGAAg                                                                                                                                >  2:744552/1‑150 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGCGGCAGGACCGGGGGGACGCAATTTTTCGCAGAggg                                                                                                                                  >  1:1155995/1‑147 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAGGt                                                                                                                               <  2:1111296/151‑1 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGt                                                                                                                               <  2:1111241/151‑1 (MQ=255)
                           ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACGAGGGTACGCCTTTTCTCGACGAGGAAGgtttt                                                                                                                           >  2:70875/1‑151 (MQ=255)
                            cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCGGCAGGACACCGGGGGCGCATTTTCTCGCCGGGGGAGgttttt                                                                                                                          >  1:227178/1‑150 (MQ=255)
                             gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGCTGCAGGACCGAGGGGACGCCTTTTCTCCACGAGGAAGttttt                                                                                                                           >  2:951881/1‑149 (MQ=255)
                                  cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTg                                                                                                                    <  2:781212/151‑1 (MQ=255)
                                   aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCAGGC‑GCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGgttttttggtt                                                                                                                     >  2:910147/1‑143 (MQ=255)
                                   aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGGTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTTGttt                                                                                                                     >  1:112056/1‑149 (MQ=255)
                                   aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGGAGGACCGAGGGGACGCATTTTTTTGACGAGGGAGgtttt                                                                                                                           >  1:1218373/1‑143 (MQ=255)
                                    gTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc                                                                                                                  <  2:1079033/151‑1 (MQ=255)
                                     ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTCGAGGCCCGAGGGGACGCATTTTTTTGCCGCGGAAGTGTTTTCGttt                                                                                                                     >  1:488057/1‑147 (MQ=255)
                                      tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGCCgg                                                                                                                >  2:110145/1‑151 (MQ=255)
                                       cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTCGGCGAGGAAGGTTTTTAGTTTTGCggg                                                                                                               >  2:259460/1‑151 (MQ=255)
                                       cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCCGCACAGCCCAGGGTAACCACTTTCTCGGCGACGCAGTTTTTCCGCTTTCCggg                                                                                                               >  1:822036/1‑151 (MQ=255)
                                           gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTc                                                                                                           <  1:1032590/151‑1 (MQ=255)
                                            tAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCCGCAGGACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGGCggg                                                                                                               >  2:941176/1‑146 (MQ=255)
                                             aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcag                                                                                                         <  1:1050524/151‑1 (MQ=255)
                                               tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGCTGCAGTACCCGGGGGACGCATTTTTTCCGCCAGGGAGGTTTTTACTTTTCCggggggca                                                                                                          >  1:1141136/1‑143 (MQ=255)
                                               tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTTGTTTGCCGGGAGTcagc                                                                                                        >  2:1069486/1‑150 (MQ=255)
                                                aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagcag                                                                                                      <  2:926950/151‑1 (MQ=255)
                                                  tGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGGCAGCAggg                                                                                                    >  2:511393/1‑151 (MQ=255)
                                                   gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTCGTTTGCCGGGAGTCAGCAggg                                                                                                    >  2:259798/1‑150 (MQ=255)
                                                     gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCtt                                                                                                 <  2:328643/151‑1 (MQ=255)
                                                     gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCtt                                                                                                 <  2:165055/151‑1 (MQ=255)
                                                       gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCt                                                                                                  <  1:864674/148‑1 (MQ=255)
                                                         ttGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATACGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCGGGGGTCAGCAGGGCttt                                                                                                >  2:7030/1‑148 (MQ=255)
                                                          tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTgg                                                                                              <  2:225809/149‑1 (MQ=255)
                                                           gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGtt                                                                                            <  2:359614/150‑1 (MQ=255)
                                                              aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGCTGCAGTACCGAGGGGACGCATTTTTTCCACGAGGAAGTTTTTTCGTTTGCCGGGGGGCAGCAGGGCTTTGGttt                                                                                           >  2:1120419/1‑148 (MQ=255)
                                                                 ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCCGTAACGCGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCGGGGGGGAGCAGGGGTTTGGTTTTCTgg                                                                                      >  2:6468/1‑150 (MQ=255)
                                                                  tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCGGGGGCAGCACCGCGGGGACGCCTTTTTTCGACGAGGAAGGTTTTTAGTTTGTCGGGGGTCAGCCGGGCTCTGGgtt                                                                                           >  1:969787/1‑144 (MQ=255)
                                                                   gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                   <  2:966898/151‑1 (MQ=255)
                                                                     cGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTTTGTTTTCCGGGGGGCAGCAGGGCTTTTGGTTCCTGGGag                                                                                   >  2:768817/1‑149 (MQ=255)
                                                                       gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                                                               <  1:291592/151‑1 (MQ=255)
                                                                        ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTCAGTTTTCCGGGGGGCAGCAGGGCCTTGGTTTTCt                                                                                        >  1:1100649/1‑141 (MQ=255)
                                                                              ggttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGTACCGAGGGGACGCATTTTCTCGGCGAGGGAGGTTTTCAGTTTTCCGGGGGTCAGCCGGGGCTTGGTTTTCTGGTAGGGATAAGtttt                                                                        >  1:634292/1‑151 (MQ=255)
                                                                                 tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgt                                                                      <  1:282891/150‑1 (MQ=255)
                                                                                  ggttTCCGCAGACGTGCCGCCAATCACAAGCGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgttg                                                                    <  2:769775/151‑1 (MQ=255)
                                                                                  ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                     <  1:974247/150‑1 (MQ=255)
                                                                                  ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                     >  2:245331/1‑150 (MQ=255)
                                                                                     ttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc                                                                  <  2:1161690/150‑1 (MQ=255)
                                                                                       ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                               <  1:234739/151‑1 (MQ=255)
                                                                                        cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCCGTACCGAGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCGGGGGTGAGCAGGGGTTTTGTTTCCTGGTAGAGATACGTTTTTTTgg                                                                   >  2:544408/1‑146 (MQ=255)
                                                                                        cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCCGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGCCGGGGGGCAGCAGGGCTTTGGGTTTCTGGTAGGGATACGTTTTGTTGGCag                                                                >  2:644584/1‑149 (MQ=255)
                                                                                             aCGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGGGCCGTACCCAGGGGACGCAATTTTTTGCCGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGGTTGCTGGTAGAGATACGTTTTGTTTGCAGAgccg                                                           >  1:651235/1‑149 (MQ=255)
                                                                                                tACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTTGACGAGGGAGGTTTTTAGTTTGCCGGGAGGGAGGAGGGCGTTGGTTTCCTTGGAGAGGTACGGTTTTTTTGGAGGgccgccgcc                                                      >  1:1006389/1‑151 (MQ=255)
                                                                                                   cGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                   <  1:135607/151‑1 (MQ=255)
                                                                                                     ccAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACCCGGGGACGCATTTTTTTCACGAGGGAGGTTTTTAGTTTTCCGGGGGTCAGCCGGGGTTTGGTTTCCTGGTAGAGATACGTGTTGTTGGGAGAGCCGCCGCCttt                                                   >  1:1159370/1‑149 (MQ=255)
                                                                                                       aaccacaAACCCGATATGGTACGGCGGG‑CAGGTTGCATTCCACAGGTTAGGCATTTTCTCGACTAGGAAGTTTTTTAGTTTGCCGGGAGTCAGCAGGTCATTCGTTTCCTGCTAGAGATACGTTTTGGTGGCAGCggcgccgcccctctca                                               <  1:43320/148‑8 (MQ=255)
                                                                                                         tCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGTACCGAGGGTACGCCTTTTCTCGGCGAGGGAGGTTTTTAGGTTGCCGGGGGGCAGCAGGGGTTTGGGTTTCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaa                                              >  1:531414/1‑150 (MQ=255)
                                                                                                          cacaAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacg                                            <  1:1175371/151‑1 (MQ=255)
                                                                                                           acaAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTGTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGGTGGCAGAGCCGCCGCCTTTcgcaacgc                                           <  1:1029313/151‑1 (MQ=255)
                                                                                                           acaAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACGGGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTTCCGGGGGGCAGCAGGGGTTTGGTTTTCTGGTAGAGATACGTTTTGTTGGCAGAACCGCCGCCTTTtgcaa                                              >  2:43015/1‑148 (MQ=255)
                                                                                                                   gATATGGTACGGCGGG‑GAGGGGGCATTACCGAGGGGGCGCATTTTTTCGACGGGGAAGATTTTCAATTTGGCGGGGGTGAGCAGGGGTTTGGGTTTCTGGTAGAGATACGTTTTGTTGGCAGCGCCGCCGCCTTTCTCAACGCACAGGaa                                    >  1:17994/1‑150 (MQ=255)
                                                                                                                   gATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTGTTTCAGTGTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATGCGTTTTGTGGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAt                                   <  1:388847/151‑1 (MQ=255)
                                                                                                                      atGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAttt                                 <  1:890964/150‑1 (MQ=255)
                                                                                                                       tGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCGCGACGAGAAAGTTGTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGATATACGTTTTGTTGGCAGAGCCGACGCCTTTCGCAACGCAAAGGAATTTg                                <  1:710884/150‑1 (MQ=255)
                                                                                                                         gTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTAc                             <  1:1046997/151‑1 (MQ=255)
                                                                                                                             ggcggG‑CAGGTTGCAGTGCCGAGGGTACGCATTTTCTCGACGAGGAAAATTTTCAGTTTGCCGGGAGGCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTGTCGCCACGCAAAGGAATTTGTACGCa                          <  1:603215/150‑1 (MQ=255)
                                                                                                                                    cAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                   >  2:511512/1‑151 (MQ=255)
                                                                                                                                     aGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                   >  1:93133/1‑150 (MQ=255)
                                                                                                                                     aGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                  >  2:470699/1‑151 (MQ=255)
                                                                                                                                       gCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTa                <  1:850302/151‑1 (MQ=255)
                                                                                                                                             gTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGt          <  2:430806/151‑1 (MQ=255)
                                                                                                                                              tACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGt          >  2:501267/1‑150 (MQ=255)
                                                                                                                                               aCCGAGGGTACGCATTTTCTCGGCGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc        <  1:1098773/151‑1 (MQ=255)
                                                                                                                                                ccGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc        <  1:704442/150‑1 (MQ=255)
                                                                                                                                                   aGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGt    >  2:1161622/1‑151 (MQ=255)
                                                                                                                                                    gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc   >  1:569082/1‑151 (MQ=255)
                                                                                                                                                     ggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCg  <  2:674636/151‑1 (MQ=255)
                                                                                                                                                     ggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCg  <  2:79648/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCG  >  NC_000913/4346525‑4346823

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