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

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGT  >  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)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGTAAGCATTTTCTCGACGAgg                                                                                                                                  >  1:207226/1‑148 (MQ=255)
                      cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCTGCCGGACCCGCGGGACGCACTTTTTCGGCGAGGAAgg                                                                                                                              >  1:970383/1‑149 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGtt                                                                                                                             <  1:1033093/149‑1 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTc                                                                                                                         <  2:688736/151‑1 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtt                                                                                                                     <  2:1059624/151‑1 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtt                                                                                                                     <  2:570756/151‑1 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCGGCAGGACCGAGGGGAAGCATTTTCTCGACGAGGGAGGTTTTCCGttt                                                                                                                    >  1:943317/1‑151 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGttt                                                                                                                    <  1:655773/151‑1 (MQ=255)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGttt                                                                                                                    <  1:791488/150‑1 (MQ=255)
                                 aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGc                                                                                                                  <  1:317322/151‑1 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCgg                                                                                                               <  2:931058/150‑1 (MQ=255)
                                          tAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTc                                                                                                          <  1:882803/150‑1 (MQ=255)
                                            aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGAAGCAACTCCGCGGGTCAGCATTTTCTCTGCGAGGAAACTTTTCGCTTTACTGGGCGTcacc                                                                                                       >  1:821259/1‑149 (MQ=255)
                                             tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGGGGcagc                                                                                                       >  2:1020808/1‑150 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagcag                                                                                                     <  2:1075962/151‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca                                                                                                      <  2:525577/150‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGTACGCCTTTTCTCGCAGAGGAAGGTTTTTTGTTTGCCggg                                                                                                              >  1:872171/1‑142 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGGGCAGGACCCAGGGGACGCCTTTTCTCGCAGAGGAAGGTTTTTCGGTTGGCGGGGGGCAGCAgg                                                                                                    >  1:333833/1‑151 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCCCAGCACCCAGGGGACGAATTTTTTCGACGGGGGAATTTTTCCGTTTTTCGGGGGTgag                                                                                                        >  1:83099/1‑147 (MQ=255)
                                                tGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGGAGTACCGAGGGGAGACCTTTTCTTGACGAGGGGGGTTTTTTGTGTGTCGGGAGGCAGCAggg                                                                                                   >  1:97693/1‑151 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAggg                                                                                                   <  1:317535/150‑1 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGAAGCATTTTTTCGACGGGGAAGTTTTTTAGTTTTCCGGGGGGCAGCAggg                                                                                                   >  2:823000/1‑150 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTTGACGAGGAAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGGtt                                                                                               >  2:565370/1‑151 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                      <  2:178901/151‑1 (MQ=255)
                                                                tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg                                                                                    <  1:841711/149‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                 <  2:718980/151‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTa                                                                                  <  1:109953/150‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGACGGGTGCAGGACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTCTGTTTTTCGGGAGTCAGCAGGGGCTTGGTTTCCTg                                                                                     >  2:1110771/1‑147 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAGCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAgaga                                                                              <  2:364691/151‑1 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCGGGAGTCAGCAGGGCCTTTGTTTTCTGGTag                                                                                 >  2:918478/1‑148 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGAGGCAGCAGGGGCTTGGTTTTCTGGTag                                                                                 >  1:546567/1‑147 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGAAGCCTTTTTTCCACGAGGGAGGTTTTTTGGTTTGCGGGGGGTAGGAGGGCTTTTGttttc                                                                                       >  1:476429/1‑135 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGtt                                                                        >  2:376117/1‑151 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGCACCGAGGGGAAGCCTTTTCTCGACGAGGGAGGTTTTCCGTTTGCCCGGAGGCAGCAGGGCTTTTGGTTCCTGGTAGAGATACGTTttgtt                                                                   >  1:121376/1‑151 (MQ=255)
                                                                                   ttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc                                                                <  2:421509/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                             <  1:752350/151‑1 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgcc                                                          <  2:1090703/151‑1 (MQ=255)
                                                                                          gACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCTACGGGGGAGTTTTTTTGTTTTGCGGGGGGCAGCAGGGGTTTGGTTTTCTGGGTGAGATTAGTTTTTTTGGCAGAgccg                                                         >  1:515661/1‑150 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGCGGGGACGCCTTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGAGGCAGCAGGGCTTTGGGTTCCTGGGAGAGATACGTTTTTTTTGGAGAgccgcc                                                       >  1:483120/1‑150 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCGTCAAGAACGAGGGGACGCAGTTTTTCCACGAAGGAATTTTTTATGTTTACGGGGGGCAGCAGGGCGTTTGTTTTCATGGAGGGAAACGTGTTGTTTGTTAAgccgacgcc                                                    >  2:314090/1‑150 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGACGCATTTTCTCGACGAGGGAGGTTTTTAGTTTGGCGGGAGGCAGCAGGGGTTTTGGTTTCTGGGAGAGATACGGTTTTTTTGCAGGggcgccgcc                                                    >  2:338482/1‑150 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCCAGGGGACGCATTTTCTCGACGAGGAAGTTTTTCCGTTTTCCGGGGGTCAGCAGGGCTTTGGGTTTCTGGTAGGGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                 >  2:786291/1‑150 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcg                                               <  1:303714/151‑1 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaac                                           <  2:2271/150‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                       <  2:932825/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCCAGGGGACGCAATTTTTCGACGAGGAAGGTTTTCAGTTTGCCCGGGGTCAGCAGGGCTTTTGTTTTCTGGTAGAGATACGTTTTTTTGGCAGAGCCGCCGGCCTTtgcacagc                                         >  1:1048223/1‑145 (MQ=255)
                                                                                                             aCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTCAGTTTTCCGGGAGGCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCaaa                                      >  2:1000610/1‑150 (MQ=255)
                                                                                                                      ggTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTa                            <  1:35307/151‑1 (MQ=255)
                                                                                                                        tACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGTCGAGGGAGGTTTTTTGTTTTCCGGGGGGCAGCAGGGGTTTTGTTTTCTGGTAGAGATACGGTTTTTTGGGAGAGCCGCCGCCTTTCGCCACGCAAAGGAAatt                               >  1:351939/1‑146 (MQ=255)
                                                                                                                                 cAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                 <  2:1081432/151‑1 (MQ=255)
                                                                                                                                    gCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCGGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                 <  1:268530/148‑1 (MQ=255)
                                                                                                                                        cAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgc          >  2:167056/1‑151 (MQ=255)
                                                                                                                                            aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa       <  2:758144/150‑1 (MQ=255)
                                                                                                                                                aGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGt  >  2:612618/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGT  >  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: