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

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  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)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCAATTTTTCGAAGAGGGAGTTTTTTAGttt                                                                                                                           >  2:1125042/1‑150 (MQ=255)
                                 aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGttt                                                                                                                           >  2:177853/1‑149 (MQ=255)
                                   ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCCAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGcc                                                                                                                        >  1:266148/1‑150 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCg                                                                                                                       <  2:699168/150‑1 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCg                                                                                                                       <  1:659024/150‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGa                                                                                                                    <  1:693065/151‑1 (MQ=255)
                                        cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGt                                                                                                                  <  1:459393/151‑1 (MQ=255)
                                        cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTACGCAGACGTCCTGCAAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCACGACCCCGGGGGCGCATCTTCTCGACTCGGGAGttttt                                                                                                                                  >  2:320654/1‑136 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGCACCGAGGGGACGCATTTTTCCTCCGAGGGAGTTTTTCAGTTTGTCggggg                                                                                                                   >  2:1259158/1‑147 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAATACCGAGGGGAGGCAATTTTTCGACGGGGAAGGTTTTCAGTTTGCCGGGAg                                                                                                                   >  1:51770/1‑147 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGGAGGACCCAGGGGAAGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGGCgggg                                                                                                                    >  1:1053322/1‑145 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTca                                                                                                                <  1:1135950/150‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagcag                                                                                                           <  1:1256212/151‑1 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAgg                                                                                                          <  2:1156690/151‑1 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGc                                                                                                       <  2:577070/151‑1 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTTAGTTTGGCGGGGGTCAGCAggg                                                                                                         >  1:1189663/1‑148 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGGGCCGGGGGGACGCATTTTTTCGGCGAGGAAGGTTTTTTGGTTTCCGGGGGGCAGCCGGGGtt                                                                                                     >  1:355589/1‑151 (MQ=255)
                                                    cgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCttt                                                                                                    <  1:418522/151‑1 (MQ=255)
                                                      cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGAGGCAGCAGGGGttt                                                                                                    >  2:996932/1‑149 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGttt                                                                                               <  1:567728/151‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGtt                                                                                                <  1:744478/150‑1 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGGAGGAACCAGGGGAAGCATTTTTTCGGCGGGGGAGGTTTTTTGTTTGGCGGGGGTCAGCAGGGGCTTGGGTTTc                                                                                              >  1:257093/1‑150 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                       >  2:153188/1‑148 (MQ=255)
                                                                       tttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATa                                                                                  <  1:334261/150‑1 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACg                                                                                <  1:758163/150‑1 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGGTTTGGTTTCCTGGGAGAGATACg                                                                                >  2:1337035/1‑150 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTTCTGGGAGAGATACGtt                                                                              >  2:458855/1‑151 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGt                                                                               <  1:276156/149‑1 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTTCGTTTGCCGGGGGGCAGCAGGGGTTTTGGTTCCTGGGAGAGATACGttt                                                                             >  1:717682/1‑151 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTTAGTTTTCCGGGGGGCAGCAGGGCTTTGGTTTTCTGGGAGAGATACGttt                                                                             >  2:114864/1‑151 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTCGGGTACCCCGGGGAAGCACTTTTTCCCCGAGGAACTTTTTTTGTTTGCCCGGGGTGAGCAGGGGTTTGGTTTTCTGGTCGAGGTACGtttt                                                                            >  2:592449/1‑150 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgttg                                                                        <  2:1335391/151‑1 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTgg                                                                       <  2:333543/151‑1 (MQ=255)
                                                                                    tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTTGTTTTCTGGGAGAGATACGTTttgtt                                                                         >  2:297204/1‑146 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGTCAGCAGGGCGTTGGTTTCCTGGTAGAGATACGTTTTGTTgg                                                                       >  2:1063704/1‑147 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTtttttt                                                                         >  2:1201192/1‑144 (MQ=255)
                                                                                       gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTTCCGGGGGGCAGCAGGGCTTTGGTTTCCTGGGAGAGATACGTTTTGTTGGCagag                                                                 >  2:618262/1‑150 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCCGGACCGAGGGGACGCAATTTCTCCACGAGGAAGGTTTTTAGTTTTCCGGGGGGCAGCAGGGGTTTGGTTTTCTGGTAGAGGTAAGGTTTTGTGGCag                                                                    >  1:921276/1‑147 (MQ=255)
                                                                                           aCGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgccgcc                                                            <  1:101332/151‑1 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTCGTTTGCCGGGGGTCAGCAGGGCTTTTGTTTTCTGGTAGAGATACGTTTTGTTGGCAGAgccgcc                                                            >  2:759303/1‑150 (MQ=255)
                                                                                             gTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTACCGCGGGGAAGCATTTTCTCCACGAGGAAGGTTTTTAGGTTGCCGGGGGGCAGCAGGGGTTTGGGTTTCTGGGAGAGATACGGTTTTTTGGGCGGgccgcc                                                            >  1:1325488/1‑149 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGGCCGAGGGGAAGCATTTTCTCGACGAGGAAGTTTTTCCGTTTGGCGGGAGCTAGCAGGGCTTTTGTTTCCTGGGAGAGGTACGTTTTGTTGTCAGAgccgccgcc                                                         >  2:685268/1‑151 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgccgccgcc                                                         <  2:72699/150‑1 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCt                                                        <  2:937177/151‑1 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTGGTTTCCTGGGAGAGATACGTTTTTTTTGCAGAgccgccgcc                                                         >  2:874511/1‑150 (MQ=255)
                                                                                                ccGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCtt                                                       <  2:903977/151‑1 (MQ=255)
                                                                                                 cGCCAATCACAAACGCGATATGGTACGGCGGACCGGGTGGAAAACCGAGGGCACGTATTTTTTCCGCGTGGGAGTTTTTTACTTTTGCCGGGGGTCGCAAGGGCTTTGTTTTCTCGTTGCGATACGTTTTTTTTGCCGAAGCGGCGGCCtt                                                       >  2:643801/1‑151 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTc                                                     <  1:216030/151‑1 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGAGGGGACGCATTTTCTTGACGAGGGAGGTTTTTAGTTTGTCGGGGGTCAGCCGGGGTTTTGTTTTCTGGGAGAGAGACGTTTTTTTGGCCGAGCCGCCCCCttt                                                      >  1:3346/1‑149 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAATACCCAGGGGAAGCATTTTTTTGGCGAGGGAGGTTTTTTGGTTTCCGGGGGGCAGCCGGGGTTTGGTTTTCTGGGAGAGATACGTTTTTTTGGCAGAGCCGCCGCCttt                                                      >  1:1178694/1‑149 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGAAGCATTTTTTCGACGAGGAAGATTTTCAGTTTTCCGGGGGTCAGCAGGGGTTTTGGTTTCTGGTAGAGATAAGGTTTTTTGGGAGAGCCGCCGCCTTTccc                                                   >  2:1297857/1‑149 (MQ=255)
                                                                                                     aaTCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTTTCGACGAGGAAGTTTTTTTGTTTTCCGGGGGTCAGCAGGGCTTTGGTTTCCTGGGAGGGATACGGTTTGTTTGCAGAGGCGCCGCCCTTcgc                                                   >  2:154137/1‑150 (MQ=255)
                                                                                                     aaTCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgca                                                  <  2:942026/151‑1 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGAATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaac                                                <  2:430524/150‑1 (MQ=255)
                                                                                                          caAACGCGATATGGTACGGCGGGCAAGCCGCAGTACCGGGGGGACGCACTTTTTCGAAGAGGGAGGTTTTTAGTTTTCCGGGGGGCAGCAGGGCCTTTGGTTTCTGGTAGGGGTAAGTTTTGTTGGCAGGAGCGCCGCCTTTtgccacgc                                             >  1:839451/1‑150 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                           <  1:1006864/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                           <  2:747132/151‑1 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCaaa                                          <  1:1216909/151‑1 (MQ=255)
                                                                                                             aCGCGATATGGTACGGCGGGCGGGGGGCAGTACCGAGGGGATGCATTTTTTTGATGAGGGAGATTTTTTGCTTGACCGGAGGAAGATGGGCTTTGGGTTTCTGGTTGAAATAGGTTTTTTTTGTAGAAGCGCCGCCTTTctgaaggcaa                                           >  1:322858/1‑149 (MQ=255)
                                                                                                                    atGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAttt                                   <  1:189082/150‑1 (MQ=255)
                                                                                                                       gTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTAc                               <  2:1141686/151‑1 (MQ=255)
                                                                                                                        tACGGCGGGCAGGGCGGAGGACCGCGGGGACGCATTTTTTTGTCGAGGAGGTTTTTTAGTTTGCCCGGGGGCCGCGGGGCGTTGGTTTCCTGGTAGAGATACGTTTTGTTGGGAGAGCCGCCGCCTTTCGCCACGCCAAAGAAAttttaat                               >  1:114139/1‑147 (MQ=255)
                                                                                                                        tACGGCGGGCAGGCTGCAGTACCGAGGGGGCGCATTTTTTCGTACAGGGAGGTTTTTAGTTTTCCGGGGGTCAGCAGGGGTTTTGGTTTCTGGGAGAGATACGTTTTTTTGGGGGAGCCGCCGCCCTTTCCCACGCAAAGGGATTTGttat                              >  1:78450/1‑148 (MQ=255)
                                                                                                                           ggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCAt                           >  2:1131277/1‑151 (MQ=255)
                                                                                                                                       gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAc                 <  2:1007956/149‑1 (MQ=255)
                                                                                                                                       gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAc                 <  1:708236/149‑1 (MQ=255)
                                                                                                                                         aGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgcg             <  2:721637/151‑1 (MQ=255)
                                                                                                                                             ccGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc          >  2:932390/1‑150 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       >  1:58925/1‑151 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       <  2:1052567/151‑1 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       >  1:627875/1‑151 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc     <  2:9180/151‑1 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc     <  2:633284/151‑1 (MQ=255)
                                                                                                                                                    tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAt  >  2:265975/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                           
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  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: