Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG43.5% 143.8 / 18.0 85T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (9/39);  new base G (37/0);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 7.20e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  NC_000913/4346527‑4346822
                                                                                                                                                    |                                                                                                                                                     
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑c                                                                                                                                                    <  2:1084683/151‑1 (MQ=255)
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑c                                                                                                                                                    <  1:457477/151‑1 (MQ=255)
   gCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGTCAGCACCGAGGGGA‑c                                                                                                                                                    >  1:397202/1‑148 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGGTGCAGGACCGCGGGGA‑AGCActtt                                                                                                                                             >  1:736433/1‑151 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCAttt                                                                                                                                              <  1:571233/150‑1 (MQ=255)
          cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCCGGACCCAGGGGA‑CGCAtttt                                                                                                                                             >  1:1080081/1‑148 (MQ=255)
             gttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGA‑CGCCTTTTTTcgac                                                                                                                                       >  1:418852/1‑151 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGA‑CGCCTTTTCTcgac                                                                                                                                       >  1:1102495/1‑150 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTcgac                                                                                                                                       <  2:989688/150‑1 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAg                                                                                                                                <  2:467836/151‑1 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGA‑CGCATTTTTTTGGCGAgg                                                                                                                                   >  1:47183/1‑148 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTACCGAGGGGA‑CGCATTTTTTCGACGAgg                                                                                                                                   >  2:580742/1‑147 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTACCGAGGGGA‑AGCACTTTTTCGACGAGGGAg                                                                                                                                >  1:361820/1‑150 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGt                                                                                                                               <  1:870572/151‑1 (MQ=255)
                      cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCCGTACCGAGGGGA‑CGCATTTTTTCGACGAGGGAGtt                                                                                                                              >  1:457072/1‑151 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGGGGAGTACCGGGGGGA‑CGCATTTTTTCGACGAGGGAGgtttt                                                                                                                           >  1:242829/1‑150 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGA‑CGCATTTTTTCGACGAGGAAGgttttt                                                                                                                          >  2:635780/1‑150 (MQ=255)
                           gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGA‑CGCATTTTTTTGACGAGGGAGgtttt                                                                                                                           >  1:1073719/1‑149 (MQ=255)
                            tCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCa                                                                                                                         <  2:540121/150‑1 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCAGTACCCAGGGGA‑CGCATTTTTTCGACGAGGGAGTTTTTTTGttt                                                                                                                     >  1:347066/1‑151 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCg                                                                                                                 <  2:403706/150‑1 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGCGGGGA‑CGCATTTTCTCCACGAGGGGGGTTTTTAGTTTTCCgg                                                                                                                >  1:907357/1‑151 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGGAGAACCCAGGGGA‑CGCCTTTTTTCGAGGAGGGAGtttttttagtttccggg                                                                                                               >  1:1113807/1‑138 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGA‑CGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCggggg                                                                                                             >  2:91141/1‑147 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGA‑CGCAATTTCTCGACGAGGGAGGTTTTTTGTTTGCCGGGAg                                                                                                             >  2:577178/1‑149 (MQ=255)
                                             tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGAACGAGGGGA‑CGCATTTTTTCGACGAGGAAGGTTTTTAGGTTTCCGGGGGTcagc                                                                                                        >  2:89189/1‑150 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca‑g                                                                                                     <  2:408599/151‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca                                                                                                       <  2:75396/150‑1 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑gg                                                                                                    <  1:421047/151‑1 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTACCGCGGGGG‑CGCATTTTCTCGGCGAGGGAGTTTTTTTGTTTGCCGGGGGGCAGCA‑ggg                                                                                                   >  1:734660/1‑150 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGCGGGGA‑CGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGGCAGCA‑ggg                                                                                                   >  1:931533/1‑148 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTg                                                                                              <  1:635717/151‑1 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTg                                                                                              <  1:956480/151‑1 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGt                                                                                            <  1:218581/150‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGTA‑AGCATTTTTTTGACGAGGAAGGTTTTTAGTTTGCCGGGAGGCAGCA‑GGggtttggttt                                                                                          >  2:31068/1‑151 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCCAGGGGA‑CGCATTTTCTCGACGAGGAAGGTTTTTAGTTTTTCGGGGGTCAGCA‑GGGCTTTGGgttt                                                                                         >  1:1083465/1‑148 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTc                                                                                         <  1:1014980/149‑1 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACCAGGGGA‑CGCATTTTTTCGACGAGGAAGTTTTTTTGTTTTCCGGGAGGCAGAAGGGGTTTTGGTTTc                                                                                         >  2:270962/1‑150 (MQ=255)
                                                             acgtGACGGTGTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGTGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCt                                                                                       <  1:18123/147‑1 (MQ=255)
                                                               ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGAACGAGGGTA‑CGCATTTTTTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGGCAGGA‑GGGGTTTGGGTTTCTgg                                                                                     >  2:510881/1‑150 (MQ=255)
                                                                tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTACCGAGGGTG‑CGCATTTTCTCGACGAGGGAGGTTTTCAGGTTGCCGGGGGTCAGCC‑GGGGTTTGGTTTTCt                                                                                       >  2:386051/1‑147 (MQ=255)
                                                                  aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTaga                                                                                 <  1:91120/151‑1 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGA‑CGCATTTTTTCGACGAGGGAGGTTTTCCGTTTGCCGGGGGGCAGCA‑GGGGTTTTGGTTTCTgg                                                                                     >  2:436697/1‑145 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGGGGGGA‑CGCCTTTTTTTGACGAGGGAGTTTTTTCGTTTTCCGGGGGTCAGCA‑GGGCTTTTGTTTTCTGGTagag                                                                                >  2:451050/1‑149 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTGGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGAt                                                                              <  1:612360/151‑1 (MQ=255)
                                                                      ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTTCCGGGGGGA‑CGCCTTTTCTCCGCGAGGAAGTTTTTCATATTGCCCGGGGTCAGCA‑GGGCttttggtt                                                                                          >  2:1103884/1‑133 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGGAGTGCCGAGGGGG‑AGCATTTTCTCCGCGAGGGGGTTTTTTAGTTTTCCGGGGGGGAGCA‑GGGGTTTGGTTTTCCGGTAGGGATACGttt                                                                        >  1:128778/1‑151 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTAACGAGGGGA‑CGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGTCAGCA‑GGGCTTTTGGTTTCTTTTAGAGATACGtttt                                                                       >  2:320551/1‑150 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCGGCAGTTCCCAGGGGAGGGCATTTT‑TTGGCGGGGGAGTTTTTTAGTTTGTCGGGAGTCAGGA‑GGGCTTTGGGTTTCTGGTGGAGATTCGtttt                                                                       >  1:460504/1‑150 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTttg                                                                      <  1:495321/151‑1 (MQ=255)
                                                                              ttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGA‑CGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGGGGCAGCA‑GGGGTTTGGTTTTCTGGTAGAGATACGtttt                                                                       >  2:328686/1‑149 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                    <  2:69508/150‑1 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTgg                                                                  <  1:83208/151‑1 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGA‑CGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCA‑GGGGTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                    >  2:155011/1‑149 (MQ=255)
                                                                                    tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCa                                                                <  1:362075/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                              <  2:1094730/151‑1 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCagag                                                             <  1:342445/151‑1 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGA‑CGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCGGGGGTCAGCA‑GGGCTTTGGTTTTCTGGTAGAGATACGTTTTGTTGGCag                                                               >  2:510615/1‑149 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGTACCGGGGGGG‑CGCAATTTTTCGGCGAGGGAGGTTTTTAGTTTGTCGGGAGGCAGCAGGGGTTTTGGTTTCTTGG‑AGAGGGACGGTTTGTTGGGag                                                               >  1:712167/1‑147 (MQ=255)
                                                                                          gACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGA‑CGCATTTTCTCCACGAGGAAGGTTTTTAGTTTGTCGGGGGGCAGCA‑GGGGGTTGGTTTTCTGGGAGAGATACGTTTTGTTGGGagag                                                             >  1:38286/1‑147 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGCGGGGA‑CGCATTTTTTTGACGAGGAAGTTTTTTAGTTTGCCGGGGGGCAGCA‑GGGGTTTGGGTTTCTGGTAGGGATACGGTTTTTTTGCCGAggcgcc                                                        >  1:607121/1‑150 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCt                                                    <  2:640645/151‑1 (MQ=255)
                                                                                                ccGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCtt                                                   <  2:991128/151‑1 (MQ=255)
                                                                                                 cGCCAATCACAAACGCGATATGGTACGGCGGGCCAGGTGCAGTAACGAGGGGA‑CGCATTTTTTCTACGGGGAAGTTTTTAAGTTTGTCTGGAGTCAGCA‑GGGGTTTGGTTTTCTGGGAGGGATACGTTTTTTTTGGAGAgccgccgcc                                                     >  2:400329/1‑148 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGAACCGAGGGGA‑CGCATTTTCTCCACGAGGAAGTTTTTTAATTTGCCGGGGGGCAGCA‑GGGGTTTGGTTTCCTGGTAGAGATAAGGTTTTTTGGCAGAGGCGCCCCCCTTtgc                                               >  2:675386/1‑151 (MQ=255)
                                                                                                       tCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTTTCGACGAGGAAGTTTTTTAGTTTTCCGGGGGGCAGCA‑GGGCTTTGGGTTCCTGGTAGAGATACGTTTTTTTTGCAGAGCCGCCGCCTTTcgcaa                                             >  2:870623/1‑150 (MQ=255)
                                                                                                          caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca                                         <  2:713111/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                        <  2:847276/151‑1 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGTA‑CGCCTTTTCTTGACGAGGGAGGTTTTTAAGTTGGCGGGGGTCAGCA‑GGGGTTTGGTTTCTTGGTAGAGATAAGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaaccc                                          >  2:71137/1‑146 (MQ=255)
                                                                                                                cgATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGA‑CGCATTTTCTCGACGAGGGAGGTTTTTTGTTTGCCGGGGGGCAGCA‑GGGCCTTTGTTTTCTGGTGGAGATACGTTTTGTTGGCGGAGCCGCCGCCTTTCGCCACGCAAAgg                                     >  1:934077/1‑149 (MQ=255)
                                                                                                                               ggCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGcc                    <  2:654961/151‑1 (MQ=255)
                                                                                                                                 cAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                  >  1:403685/1‑151 (MQ=255)
                                                                                                                                   ggCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCt                <  2:412969/151‑1 (MQ=255)
                                                                                                                                   ggCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCt                <  1:343144/151‑1 (MQ=255)
                                                                                                                                            aCCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc             <  1:679139/145‑1 (MQ=255)
                                                                                                                                            aCCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc             >  2:679202/1‑145 (MQ=255)
                                                                                                                                            aCCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa        <  1:877576/150‑1 (MQ=255)
                                                                                                                                              cGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAg     <  1:451396/151‑1 (MQ=255)
                                                                                                                                              cGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAg     <  1:675323/151‑1 (MQ=255)
                                                                                                                                               gAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg    <  2:610571/151‑1 (MQ=255)
                                                                                                                                               gAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg    <  2:19797/151‑1 (MQ=255)
                                                                                                                                                aGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg    >  1:424809/1‑150 (MQ=255)
                                                                                                                                                 gggTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  <  2:742705/151‑1 (MQ=255)
                                                                                                                                                 gggTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  <  2:1083564/151‑1 (MQ=255)
                                                                                                                                                 gggTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  >  1:876166/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTA‑CGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCA‑GGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  NC_000913/4346527‑4346822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑