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

GGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCG  >  NC_000913/4346525‑4346823
                                                                                                                                                       |                                                                                                                                                    
ggAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGt                                                                                                                                                      <  2:573106/151‑1 (MQ=255)
   aCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACg                                                                                                                                                   <  2:1224831/151‑1 (MQ=255)
    cGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGAAGGGTGGAGAACCGCGGGGAc                                                                                                                                                    >  1:430679/1‑149 (MQ=255)
        tGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCAt                                                                                                                                                <  1:586193/149‑1 (MQ=255)
         gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCAtttt                                                                                                                                             <  2:1213463/151‑1 (MQ=255)
         gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCAtt                                                                                                                                               <  1:681295/149‑1 (MQ=255)
          aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtc                                                                                                                                            <  1:823128/151‑1 (MQ=255)
          aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTGGAAGGTTTGCAGGTTGGTCTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTGtc                                                                                                                                            <  1:460875/151‑1 (MQ=255)
            cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                          <  2:1222069/151‑1 (MQ=255)
                  cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCTGCCGTACCGGGGGGACGCCTTTTCTccacggg                                                                                                                                    >  1:1223163/1‑149 (MQ=255)
                  cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAg                                                                                                                                    <  2:1190822/151‑1 (MQ=255)
                   gttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGCACCGAGGGGACGCATTTTCTCGACGCgg                                                                                                                                   >  1:292213/1‑151 (MQ=255)
                   gttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGAAGCATTTTTTCGACGAgg                                                                                                                                   >  2:1129455/1‑151 (MQ=255)
                    ttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACCAGGGGACGCATTTTCTCGACGAgg                                                                                                                                   >  2:850538/1‑150 (MQ=255)
                     tcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGa                                                                                                                                  <  1:68127/150‑1 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGAAACGGGGGGACGCATTTTTTCGACGAGGAAg                                                                                                                                >  2:744552/1‑150 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCGGCAGGACCGGGGGGACGCAATTTTTCGCAGAggg                                                                                                                                  >  1:1155995/1‑147 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAGGt                                                                                                                               <  2:1111296/151‑1 (MQ=255)
                       ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGt                                                                                                                               <  2:1111241/151‑1 (MQ=255)
                           ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGTACGCCTTTTCTCGACGAGGAAGgtttt                                                                                                                           >  2:70875/1‑151 (MQ=255)
                            cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGGACACCGGGGGCGCATTTTCTCGCCGGGGGAGgttttt                                                                                                                          >  1:227178/1‑150 (MQ=255)
                             gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGACCGAGGGGACGCCTTTTCTCCACGAGGAAGttttt                                                                                                                           >  2:951881/1‑149 (MQ=255)
                                  cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTg                                                                                                                    <  2:781212/151‑1 (MQ=255)
                                   aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCAGGCGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGgttttttggtt                                                                                                                     >  2:910147/1‑143 (MQ=255)
                                   aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTTGttt                                                                                                                     >  1:112056/1‑149 (MQ=255)
                                   aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGGAGGACCGAGGGGACGCATTTTTTTGACGAGGGAGgtttt                                                                                                                           >  1:1218373/1‑143 (MQ=255)
                                    gTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc                                                                                                                  <  2:1079033/151‑1 (MQ=255)
                                     ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTCGAGGCCCGAGGGGACGCATTTTTTTGCCGCGGAAGTGTTTTCGttt                                                                                                                     >  1:488057/1‑147 (MQ=255)
                                      tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGCCgg                                                                                                                >  2:110145/1‑151 (MQ=255)
                                       cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGGCGAGGAAGGTTTTTAGTTTTGCggg                                                                                                               >  2:259460/1‑151 (MQ=255)
                                       cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCACAGCCCAGGGTAACCACTTTCTCGGCGACGCAGTTTTTCCGCTTTCCggg                                                                                                               >  1:822036/1‑151 (MQ=255)
                                           gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTc                                                                                                           <  1:1032590/151‑1 (MQ=255)
                                            tAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGGCggg                                                                                                               >  2:941176/1‑146 (MQ=255)
                                             aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcag                                                                                                         <  1:1050524/151‑1 (MQ=255)
                                               tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCCGGGGGACGCATTTTTTCCGCCAGGGAGGTTTTTACTTTTCCggggggca                                                                                                          >  1:1141136/1‑143 (MQ=255)
                                               tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTTGTTTGCCGGGAGTcagc                                                                                                        >  2:1069486/1‑150 (MQ=255)
                                                aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagcag                                                                                                      <  2:926950/151‑1 (MQ=255)
                                                  tGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGGCAGCAggg                                                                                                    >  2:511393/1‑151 (MQ=255)
                                                   gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTCGTTTGCCGGGAGTCAGCAggg                                                                                                    >  2:259798/1‑150 (MQ=255)
                                                     gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCtt                                                                                                 <  2:328643/151‑1 (MQ=255)
                                                     gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCtt                                                                                                 <  2:165055/151‑1 (MQ=255)
                                                       gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCt                                                                                                  <  1:864674/148‑1 (MQ=255)
                                                         ttGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATACGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCGGGGGTCAGCAGGGCttt                                                                                                >  2:7030/1‑148 (MQ=255)
                                                          tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTgg                                                                                              <  2:225809/149‑1 (MQ=255)
                                                           gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGtt                                                                                            <  2:359614/150‑1 (MQ=255)
                                                              aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGACGCATTTTTTCCACGAGGAAGTTTTTTCGTTTGCCGGGGGGCAGCAGGGCTTTGGttt                                                                                           >  2:1120419/1‑148 (MQ=255)
                                                                 ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTAACGCGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCGGGGGGGAGCAGGGGTTTGGTTTTCTgg                                                                                      >  2:6468/1‑150 (MQ=255)
                                                                  tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGGGCAGCACCGCGGGGACGCCTTTTTTCGACGAGGAAGGTTTTTAGTTTGTCGGGGGTCAGCCGGGCTCTGGgtt                                                                                           >  1:969787/1‑144 (MQ=255)
                                                                   gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                   <  2:966898/151‑1 (MQ=255)
                                                                     cGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTTTGTTTTCCGGGGGGCAGCAGGGCTTTTGGTTCCTGGGag                                                                                   >  2:768817/1‑149 (MQ=255)
                                                                       gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                                                               <  1:291592/151‑1 (MQ=255)
                                                                        ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTCAGTTTTCCGGGGGGCAGCAGGGCCTTGGTTTTCt                                                                                        >  1:1100649/1‑141 (MQ=255)
                                                                              ggttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTCTCGGCGAGGGAGGTTTTCAGTTTTCCGGGGGTCAGCCGGGGCTTGGTTTTCTGGTAGGGATAAGtttt                                                                        >  1:634292/1‑151 (MQ=255)
                                                                                 tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgt                                                                      <  1:282891/150‑1 (MQ=255)
                                                                                  ggttTCCGCAGACGTGCCGCCAATCACAAGCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgttg                                                                    <  2:769775/151‑1 (MQ=255)
                                                                                  ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                     <  1:974247/150‑1 (MQ=255)
                                                                                  ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                     >  2:245331/1‑150 (MQ=255)
                                                                                     ttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc                                                                  <  2:1161690/150‑1 (MQ=255)
                                                                                       ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                               <  1:234739/151‑1 (MQ=255)
                                                                                        cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTACCGAGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCGGGGGTGAGCAGGGGTTTTGTTTCCTGGTAGAGATACGTTTTTTTgg                                                                   >  2:544408/1‑146 (MQ=255)
                                                                                        cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGCCGGGGGGCAGCAGGGCTTTGGGTTTCTGGTAGGGATACGTTTTGTTGGCag                                                                >  2:644584/1‑149 (MQ=255)
                                                                                             aCGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCCGTACCCAGGGGACGCAATTTTTTGCCGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGGTTGCTGGTAGAGATACGTTTTGTTTGCAGAgccg                                                           >  1:651235/1‑149 (MQ=255)
                                                                                                tACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTTGACGAGGGAGGTTTTTAGTTTGCCGGGAGGGAGGAGGGCGTTGGTTTCCTTGGAGAGGTACGGTTTTTTTGGAGGgccgccgcc                                                      >  1:1006389/1‑151 (MQ=255)
                                                                                                   cGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                   <  1:135607/151‑1 (MQ=255)
                                                                                                     ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACCCGGGGACGCATTTTTTTCACGAGGGAGGTTTTTAGTTTTCCGGGGGTCAGCCGGGGTTTGGTTTCCTGGTAGAGATACGTGTTGTTGGGAGAGCCGCCGCCttt                                                   >  1:1159370/1‑149 (MQ=255)
                                                                                                       aaccacaAACCCGATATGGTACGGCGGGCAGGTTGCATTCCACAGGTTAGGCATTTTCTCGACTAGGAAGTTTTTTAGTTTGCCGGGAGTCAGCAGGTCATTCGTTTCCTGCTAGAGATACGTTTTGGTGGCAGCggcgccgcccctctca                                               <  1:43320/148‑8 (MQ=255)
                                                                                                         tCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGTACGCCTTTTCTCGGCGAGGGAGGTTTTTAGGTTGCCGGGGGGCAGCAGGGGTTTGGGTTTCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaa                                              >  1:531414/1‑150 (MQ=255)
                                                                                                          cacaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacg                                            <  1:1175371/151‑1 (MQ=255)
                                                                                                           acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTGTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGGTGGCAGAGCCGCCGCCTTTcgcaacgc                                           <  1:1029313/151‑1 (MQ=255)
                                                                                                           acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGGGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTTCCGGGGGGCAGCAGGGGTTTGGTTTTCTGGTAGAGATACGTTTTGTTGGCAGAACCGCCGCCTTTtgcaa                                              >  2:43015/1‑148 (MQ=255)
                                                                                                                   gATATGGTACGGCGGGGAGGGGGCATTACCGAGGGGGCGCATTTTTTCGACGGGGAAGATTTTCAATTTGGCGGGGGTGAGCAGGGGTTTGGGTTTCTGGTAGAGATACGTTTTGTTGGCAGCGCCGCCGCCTTTCTCAACGCACAGGaa                                    >  1:17994/1‑150 (MQ=255)
                                                                                                                   gATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTGTTTCAGTGTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATGCGTTTTGTGGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAt                                   <  1:388847/151‑1 (MQ=255)
                                                                                                                      atGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAttt                                 <  1:890964/150‑1 (MQ=255)
                                                                                                                       tGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCGCGACGAGAAAGTTGTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGATATACGTTTTGTTGGCAGAGCCGACGCCTTTCGCAACGCAAAGGAATTTg                                <  1:710884/150‑1 (MQ=255)
                                                                                                                         gTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTAc                             <  1:1046997/151‑1 (MQ=255)
                                                                                                                             ggcggGCAGGTTGCAGTGCCGAGGGTACGCATTTTCTCGACGAGGAAAATTTTCAGTTTGCCGGGAGGCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTGTCGCCACGCAAAGGAATTTGTACGCa                          <  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)
                                                                                                                                                       |                                                                                                                                                    
GGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCG  >  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: