Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG46.0% 143.9 / 24.8 88T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (6/41);  new base G (40/0);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 9.42e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  NC_000913/4346527‑4346825
                                                                                                                                                    |                                                                                                                                                           
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                           <  2:955992/151‑1 (MQ=255)
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGAc                                                                                                                                                           >  2:358541/1‑151 (MQ=255)
     cTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCGAGGGGACGCAtt                                                                                                                                                      >  1:391041/1‑151 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAAGACCGAGGGGACGCAtttt                                                                                                                                                    >  2:827646/1‑151 (MQ=255)
        aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCCGGACCGAGGGGACGCAtttt                                                                                                                                                    >  2:792477/1‑150 (MQ=255)
         ccATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtct                                                                                                                                                  <  2:529222/151‑1 (MQ=255)
          cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                                 <  1:758866/151‑1 (MQ=255)
          cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtct                                                                                                                                                  >  2:1105316/1‑150 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGCGCAGTAACCAGGGGACGCATTTTTTcgacg                                                                                                                                             >  1:405336/1‑151 (MQ=255)
               tcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgacg                                                                                                                                             <  2:894627/150‑1 (MQ=255)
                cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCCGGACCGAGGGGACGCAtttttttcac                                                                                                                                              >  2:189085/1‑142 (MQ=255)
                cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgacga                                                                                                                                            <  2:917981/150‑1 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAg                                                                                                                                       <  2:1177874/151‑1 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGtt                                                                                                                                     <  2:650857/149‑1 (MQ=255)
                           gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTAACCCGGGGACGCATTTTTTTGGCCAGGGAGgtttt                                                                                                                                  >  1:407787/1‑149 (MQ=255)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCAATTTTTCGAAGAGGGAGTTTTTTA‑Gttt                                                                                                                           >  2:1125042/1‑150 (MQ=255)
                                 aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCA‑Gttt                                                                                                                           >  2:177853/1‑149 (MQ=255)
                                   ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCCAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCA‑GTTTGcc                                                                                                                        >  1:266148/1‑150 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCg                                                                                                                       <  2:699168/150‑1 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCg                                                                                                                       <  1:659024/150‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGa                                                                                                                    <  1:693065/151‑1 (MQ=255)
                                        cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGt                                                                                                                  <  1:459393/151‑1 (MQ=255)
                                        cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTACGCAGACGTCCTGCAAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCACGACCCCGGGGGCGCATCTTCTCGACTCGGGAGttttt                                                                                                                                  >  2:320654/1‑136 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGCACCGAGGGGACGCATTTTTCCTCCGAGGGAGTTTTTCA‑GTTTGTCggggg                                                                                                                   >  2:1259158/1‑147 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAATACCGAGGGGAGGCAATTTTTCGACGGGGAAGGTTTTCA‑GTTTGCCGGGAg                                                                                                                   >  1:51770/1‑147 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGGAGGACCCAGGGGAAGCATTTTTTCGACGAGGGAGGTTTTTA‑GTTTGGCgggg                                                                                                                    >  1:1053322/1‑145 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTca                                                                                                                <  1:1135950/150‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTcagca‑g                                                                                                           <  1:1256212/151‑1 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑gg                                                                                                          <  2:1156690/151‑1 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑c                                                                                                       <  2:577070/151‑1 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTTA‑GTTTGGCGGGGGTCAGCA‑ggg                                                                                                         >  1:1189663/1‑148 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGGGCCGGGGGGACGCATTTTTTCGGCGAGGAAGGTTTTTT‑GGTTTCCGGGGGGCAGCC‑GGG‑Gtt                                                                                                     >  1:355589/1‑151 (MQ=255)
                                                    cgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑Cttt                                                                                                    <  1:418522/151‑1 (MQ=255)
                                                      cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTA‑GTTTGCCGGGAGGCAGCA‑GGG‑Gttt                                                                                                    >  2:996932/1‑149 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGttt                                                                                               <  1:567728/151‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGtt                                                                                                <  1:744478/150‑1 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGGAGGAACCAGGGGAAGCATTTTTTCGGCGGGGGAGGTTTTTT‑GTTTGGCGGGGGTCAGCAGGGG‑CTTGGGTTTc                                                                                              >  1:257093/1‑150 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTag                                                                                       >  2:153188/1‑148 (MQ=255)
                                                                       tttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATa                                                                                  <  1:334261/150‑1 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACg                                                                                <  1:758163/150‑1 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTA‑GTTTGCCGGGGGTCAGCA‑GGG‑GTTTGGTTTCCTGGGAGAGATACg                                                                                >  2:1337035/1‑150 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTA‑GTTTGCCGGGGGGCAGCA‑GGG‑CTTTGGTTTTCTGGGAGAGATACGtt                                                                              >  2:458855/1‑151 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGt                                                                               <  1:276156/149‑1 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTTC‑GTTTGCCGGGGGGCAGCA‑GGG‑GTTTTGGTTCCTGGGAGAGATACGttt                                                                             >  1:717682/1‑151 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTTA‑GTTTTCCGGGGGGCAGCA‑GGG‑CTTTGGTTTTCTGGGAGAGATACGttt                                                                             >  2:114864/1‑151 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTCGGGTACCCCGGGGAAGCACTTTTTCCCCGAGGAACTTTTTTT‑GTTTGCCCGGGGTGAGCA‑GGG‑GTTTGGTTTTCTGGTCGAGGTACGtttt                                                                            >  2:592449/1‑150 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTttgttg                                                                        <  2:1335391/151‑1 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTgg                                                                       <  2:333543/151‑1 (MQ=255)
                                                                                    tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTA‑GTTTGCCGGGGGTCAGCA‑GGG‑CTTTTGTTTTCTGGGAGAGATACGTTttgtt                                                                         >  2:297204/1‑146 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCA‑GTTTGCCGGGGGTCAGCA‑GGG‑CGTTGGTTTCCTGGTAGAGATACGTTTTGTTgg                                                                       >  2:1063704/1‑147 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTtttttt                                                                         >  2:1201192/1‑144 (MQ=255)
                                                                                       gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTA‑GTTTTCCGGGGGGCAGCA‑GGG‑CTTTGGTTTCCTGGGAGAGATACGTTTTGTTGGCag‑ag                                                                 >  2:618262/1‑150 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCCGGACCGAGGGGACGCAATTTCTCCACGAGGAAGGTTTTTA‑GTTTTCCGGGGGGCAGCA‑GGG‑GTTTGGTTTTCTGGTAGAGGTAAGGTTTTGTGGCag                                                                    >  1:921276/1‑147 (MQ=255)
                                                                                           aCGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑Agccgcc                                                            <  1:101332/151‑1 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTC‑GTTTGCCGGGGGTCAGCA‑GGG‑CTTTTGTTTTCTGGTAGAGATACGTTTTGTTGGCAG‑Agccgcc                                                            >  2:759303/1‑150 (MQ=255)
                                                                                             gTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTACCGCGGGGAAGCATTTTCTCCACGAGGAAGGTTTTTA‑GGTTGCCGGGGGGCAGCA‑GGG‑GTTTGGGTTTCTGGGAGAGATACGGTTTTTTGGGCG‑Ggccgcc                                                            >  1:1325488/1‑149 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGGCCGAGGGGAAGCATTTTCTCGACGAGGAAGTTTTTCC‑GTTTGGCGGGAGCTAGCA‑GGG‑CTTTTGTTTCCTGGGAGAGGTACGTTTTGTTGTCAG‑Agccgccgcc                                                         >  2:685268/1‑151 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑Agccgccgcc                                                         <  2:72699/150‑1 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCt                                                        <  2:937177/151‑1 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTA‑GTTTGCCGGGGGTCAGCA‑GGG‑CTTTGGTTTCCTGGGAGAGATACGTTTTTTTTGCAG‑Agccgccgcc                                                         >  2:874511/1‑150 (MQ=255)
                                                                                                ccGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCtt                                                       <  2:903977/151‑1 (MQ=255)
                                                                                                 cGCCAATCACAAACGCGATATGGTACGGCGGACCGGGTGGAAAACCGAGGGCACGTATTTTTTCCGCGTGGGAGTTTTTTACTTTTGCCGGGGGTCGCAA‑GGG‑CTTT‑GTTTTCTCGTTGCGATACGTTTTTTTTGCCGAAGCGGCGGCCtt                                                       >  2:643801/1‑151 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTc                                                     <  1:216030/151‑1 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGAGGGGACGCATTTTCTTGACGAGGGAGGTTTTTA‑GTTTGTCGGGGGTCAGCC‑GGG‑GTTTTGTTTTCTGGGAGAGAGACGTTTTTTTGGCCG‑AGCCGCCCCCttt                                                      >  1:3346/1‑149 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAATACCCAGGGGAAGCATTTTTTTGGCGAGGGAGGTTTTTT‑GGTTTCCGGGGGGCAGCC‑GGG‑GTTTGGTTTTCTGGGAGAGATACGTTTTTTTGGCAG‑AGCCGCCGCCttt                                                      >  1:1178694/1‑149 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGAAGCATTTTTTCGACGAGGAAGATTTTCA‑GTTTTCCGGGGGTCAGCA‑GGG‑GTTTTGGTTTCTGGTAGAGATAAGGTTTTTTGGGAG‑AGCCGCCGCCTTTccc                                                   >  2:1297857/1‑149 (MQ=255)
                                                                                                     aaTCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTTTCGACGAGGAAGTTTTTTT‑GTTTTCCGGGGGTCAGCA‑GGG‑CTTTGGTTTCCTGGGAGGGATACGGTTTGTTTGCAG‑AGGCGCCGCCCTTcgc                                                   >  2:154137/1‑150 (MQ=255)
                                                                                                     aaTCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTcgca                                                  <  2:942026/151‑1 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGAATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTcgcaac                                                <  2:430524/150‑1 (MQ=255)
                                                                                                          caAACGCGATATGGTACGGCGGGCAAGCCGCAGTACCGGGGGGACGCACTTTTTCGAAGAGGGAGGTTTTTA‑GTTTTCCGGGGGGCAGCA‑GGGCCTTTGGTTT‑CTGGTAGGGGTAAGTTTTGTTGGCAG‑GAGCGCCGCCTTTtgccacg‑c                                             >  1:839451/1‑150 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTcgcaacg‑caa                                           <  1:1006864/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTcgcaacg‑caa                                           <  2:747132/151‑1 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑Caaa                                          <  1:1216909/151‑1 (MQ=255)
                                                                                                             aCGCGATATGGTACGGCGGGCGGGGGGCAGTACCGAGGGGATGCATTTTTTTGATGAGGGAGATTTTTT‑GCTTGACCGGAGGAAGAT‑GGG‑CTTTGGGTTTCTGGTTGAAATAGGTTTTTTTTGTAG‑AAGCGCCGCCTTTctgaagg‑caa                                           >  1:322858/1‑149 (MQ=255)
                                                                                                                    atGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAAttt                                   <  1:189082/150‑1 (MQ=255)
                                                                                                                       gTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTAc                               <  2:1141686/151‑1 (MQ=255)
                                                                                                                        tACGGCGGGCAGGGCGGAGGACCGCGGGGACGCATTTTTTTGTCGAGGAGGTTTTTTA‑GTTTGCCCGGGGGCCGCG‑GGG‑CGTTGGTTTCCTGGTAGAGATACGTTTTGTTGGGAG‑AGCCGCCGCCTTTCGCCACGCCAAAGAAAttttaat                               >  1:114139/1‑147 (MQ=255)
                                                                                                                        tACGGCGGGCAGGCTGCAGTACCGAGGGGGCGCATTTTTTCGTACAGGGAGGTTTTTA‑GTTTTCCGGGGGTCAGCA‑GGG‑GTTTTGGTTTCTGGGAGAGATACGTTTTTTTGGGGG‑AGCCGCCGCCCTTTCCCACG‑CAAAGGGATTTGttat                              >  1:78450/1‑148 (MQ=255)
                                                                                                                           ggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCAt                           >  2:1131277/1‑151 (MQ=255)
                                                                                                                                       gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTAc                 <  2:1007956/149‑1 (MQ=255)
                                                                                                                                       gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTAc                 <  1:708236/149‑1 (MQ=255)
                                                                                                                                         aGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACcgcg             <  2:721637/151‑1 (MQ=255)
                                                                                                                                             ccGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc          >  2:932390/1‑150 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       >  1:58925/1‑151 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       <  2:1052567/151‑1 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       >  1:627875/1‑151 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc     <  2:9180/151‑1 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc     <  2:633284/151‑1 (MQ=255)
                                                                                                                                                    tACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAt  >  2:265975/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                           
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCA‑GTTTGCCGGGAGTCAGCA‑GGG‑CTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAG‑AGCCGCCGCCTTTCGCAACG‑CAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  NC_000913/4346527‑4346825

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: ‑