Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG40.7% 98.5 / 16.3 59T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (5/30);  new base G (24/0);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 5.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGT  >  NC_000913/4346527‑4346821
                                                                                                                                                    |                                                                                                                                                     
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                     <  2:1028153/151‑1 (MQ=255)
    ccTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAGCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCAt                                                                                                                                                 <  1:1010287/151‑1 (MQ=255)
        aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACCAGGGGACGCAtttt                                                                                                                                              >  2:905112/1‑150 (MQ=255)
         ccATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtct                                                                                                                                            <  1:1086458/151‑1 (MQ=255)
               tcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCGATTCCGGGGGGGGGCAtttt                                                                                                                                              >  2:505624/1‑143 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGTAAGCATTTTCTCGACG‑A‑gg                                                                                                                                  >  1:207226/1‑148 (MQ=255)
                      cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCTGCCGGACCCGCGGGACGCACTTTTTCGGCG‑A‑GGAAgg                                                                                                                              >  1:970383/1‑149 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGtt                                                                                                                             <  1:1033093/149‑1 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTc                                                                                                                         <  2:688736/151‑1 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGtt                                                                                                                     <  2:1059624/151‑1 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGtt                                                                                                                     <  2:570756/151‑1 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCGGCAGGACCGAGGGGAAGCATTTTCTCGACG‑A‑GGGAGGTTTTCCGttt                                                                                                                    >  1:943317/1‑151 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGttt                                                                                                                    <  1:655773/151‑1 (MQ=255)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGttt                                                                                                                    <  1:791488/150‑1 (MQ=255)
                                 aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGc                                                                                                                  <  1:317322/151‑1 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCgg                                                                                                               <  2:931058/150‑1 (MQ=255)
                                          tAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTc                                                                                                          <  1:882803/150‑1 (MQ=255)
                                            aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGAAGCAACTCCGCGGGTCAGCATTTTCTCTGCG‑A‑GGAAACTTTTCGCTTTACTGGGCGTcacc                                                                                                       >  1:821259/1‑149 (MQ=255)
                                             tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTAACGAGGGGACGCATTTTTTCGACG‑A‑GGAAGTTTTTTAGTTTGCCGGGGGGcagc                                                                                                       >  2:1020808/1‑150 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTcagcag                                                                                                     <  2:1075962/151‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTcagca                                                                                                      <  2:525577/150‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGTACGCCTTTTCTCGCAG‑A‑GGAAGGTTTTTTGTTTGCCggg                                                                                                              >  1:872171/1‑142 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGGGCAGGACCCAGGGGACGCCTTTTCTCGCAG‑A‑GGAAGGTTTTTCGGTTGGCGGGGGGCAGCAgg                                                                                                    >  1:333833/1‑151 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCCCAGCACCCAGGGGACGAATTTTTTCGACG‑G‑GGGAATTTTTCCGTTTTTCGGGGGTgag                                                                                                        >  1:83099/1‑147 (MQ=255)
                                                tGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGGAGTACCGAGGGGAGACCTTTTCTTGACG‑A‑GGGGGGTTTTTTGTGTGTCGGGAGGCAGCAggg                                                                                                   >  1:97693/1‑151 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAggg                                                                                                   <  1:317535/150‑1 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGAAGCATTTTTTCGACG‑G‑GGAAGTTTTTTAGTTTTCCGGGGGGCAGCAggg                                                                                                   >  2:823000/1‑150 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTTGACG‑A‑GGAAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGG‑tt                                                                                               >  2:565370/1‑151 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCt                                                                                      <  2:178901/151‑1 (MQ=255)
                                                                tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTgg                                                                                    <  1:841711/149‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTag                                                                                 <  2:718980/151‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTa                                                                                  <  1:109953/150‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGACGGGTGCAGGACCGAGGGGACGCATTTTCTCGACG‑A‑GGGAGTTTTTCTGTTTTTCGGGAGTCAGCAGGGG‑CTTGGTTTCCTg                                                                                     >  2:1110771/1‑147 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAGCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAgaga                                                                              <  2:364691/151‑1 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACG‑A‑GGGAGTTTTTTAGTTTGCCGGGAGTCAGCAGGGC‑CTTTGTTTTCTGGTag                                                                                 >  2:918478/1‑148 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGACCGAGGGGACGCATTTTTTCGACG‑A‑GGGAGGTTTTTAGTTTGCCGGGAGGCAGCAGGGG‑CTTGGTTTTCTGGTag                                                                                 >  1:546567/1‑147 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGAAGCCTTTTTTCCACG‑AGGGAGGTTTTTTGGTTTGCGGGGGGT‑AGGAGGGC‑TTTTGttttc                                                                                       >  1:476429/1‑135 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACG‑A‑GGAAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGtt                                                                        >  2:376117/1‑151 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGCACCGAGGGGAAGCCTTTTCTCGACG‑A‑GGGAGGTTTTCCGTTTGCCCGGAGGCAGCAGGGC‑TTTTGGTTCCTGGTAGAGATACGTTttgtt                                                                   >  1:121376/1‑151 (MQ=255)
                                                                                   ttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc                                                                <  2:421509/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                             <  1:752350/151‑1 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgcc                                                          <  2:1090703/151‑1 (MQ=255)
                                                                                          gACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCTACG‑G‑GGGAGTTTTTTTGTTTTGCGGGGGGCAGCAGGGG‑TTTGGTTTTCTGGGTGAGATTAGTTTTTTTGGCAGAgccg                                                         >  1:515661/1‑150 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGCGGGGACGCCTTTTTTCGACG‑A‑GGAAGGTTTTTAGTTTGCCGGGAGGCAGCAGGGC‑TTTGGGTTCCTGGGAGAGATACGTTTTTTTTGGAGAgccgcc                                                       >  1:483120/1‑150 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCGTCAAGAACGAGGGGACGCAGTTTTTCCACGAA‑GGAATTTTTTATGTTT‑ACGGGGGGCAGCAGGGCGTTTGTTTTCATGG‑AGGGAAACGTGTTGTTTGTTAAgccgacgcc                                                    >  2:314090/1‑150 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGACGCATTTTCTCGACG‑A‑GGGAGGTTTTTAGTTTGGCGGGAGGCAGCAGGGG‑TTTTGGTTTCTGGGAGAGATACGGTTTTTTTGCAGGggcgccgcc                                                    >  2:338482/1‑150 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCCAGGGGACGCATTTTCTCGACG‑A‑GGAAGTTTTTCCGTTTTCCGGGGGTCAGCAGGGC‑TTTGGGTTTCTGGTAGGGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                 >  2:786291/1‑150 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcg                                               <  1:303714/151‑1 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaac                                           <  2:2271/150‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                       <  2:932825/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCCAGGGGACGCAATTTTTCGACG‑A‑GGAAGGTTTTCAGTTTGCCCGGGGTCAGCAGGGC‑TTTTGTTTTCTGGTAGAGATACGTTTTTTTGGCAGAGCCGCCGGCCTTtgcacagc                                         >  1:1048223/1‑145 (MQ=255)
                                                                                                             aCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACG‑A‑GGAAGGTTTTCAGTTTTCCGGGAGGCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCaaa                                      >  2:1000610/1‑150 (MQ=255)
                                                                                                                      ggTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTa                            <  1:35307/151‑1 (MQ=255)
                                                                                                                        tACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGTCG‑A‑GGGAGGTTTTTTGTTTTCCGGGGGGCAGCAGGGG‑TTTTGTTTTCTGGTAGAGATACGGTTTTTTGGGAGAGCCGCCGCCTTTCGCCACGCAAAGGAAatt                               >  1:351939/1‑146 (MQ=255)
                                                                                                                                 cAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                 <  2:1081432/151‑1 (MQ=255)
                                                                                                                                    gCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCGGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                 <  1:268530/148‑1 (MQ=255)
                                                                                                                                        cAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgc          >  2:167056/1‑151 (MQ=255)
                                                                                                                                            aCCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa       <  2:758144/150‑1 (MQ=255)
                                                                                                                                                aGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGt  >  2:612618/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACG‑A‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGC‑TTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGT  >  NC_000913/4346527‑4346821

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