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

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  NC_000913/4346527‑4346822
                                                                                                                                                    |                                                                                                                                                     
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                    <  2:1084683/151‑1 (MQ=255)
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                    <  1:457477/151‑1 (MQ=255)
   gCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGTCAGCACCGAGGGGAc                                                                                                                                                    >  1:397202/1‑148 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGGTGCAGGACCGCGGGGAAGCActtt                                                                                                                                             >  1:736433/1‑151 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCAttt                                                                                                                                              <  1:571233/150‑1 (MQ=255)
          cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCCGGACCCAGGGGACGCAtttt                                                                                                                                             >  1:1080081/1‑148 (MQ=255)
             gttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTcgac                                                                                                                                       >  1:418852/1‑151 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGACGCCTTTTCTcgac                                                                                                                                       >  1:1102495/1‑150 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgac                                                                                                                                       <  2:989688/150‑1 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAg                                                                                                                                <  2:467836/151‑1 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGACGCATTTTTTTGGCGAgg                                                                                                                                   >  1:47183/1‑148 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAgg                                                                                                                                   >  2:580742/1‑147 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTACCGAGGGGAAGCACTTTTTCGACGAGGGAg                                                                                                                                >  1:361820/1‑150 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGt                                                                                                                               <  1:870572/151‑1 (MQ=255)
                      cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCCGTACCGAGGGGACGCATTTTTTCGACGAGGGAGtt                                                                                                                              >  1:457072/1‑151 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGGGGAGTACCGGGGGGACGCATTTTTTCGACGAGGGAGgtttt                                                                                                                           >  1:242829/1‑150 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGgttttt                                                                                                                          >  2:635780/1‑150 (MQ=255)
                           gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTTGACGAGGGAGgtttt                                                                                                                           >  1:1073719/1‑149 (MQ=255)
                            tCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCa                                                                                                                         <  2:540121/150‑1 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCAGTACCCAGGGGACGCATTTTTTCGACGAGGGAGTTTTTTTGttt                                                                                                                     >  1:347066/1‑151 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCg                                                                                                                 <  2:403706/150‑1 (MQ=255)
                                    tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGCGGGGACGCATTTTCTCCACGAGGGGGGTTTTTAGTTTTCCgg                                                                                                                >  1:907357/1‑151 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGGAGAACCCAGGGGACGCCTTTTTTCGAGGAGGGAGtttttttagtttccggg                                                                                                               >  1:1113807/1‑138 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCggggg                                                                                                             >  2:91141/1‑147 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGACGCAATTTCTCGACGAGGGAGGTTTTTTGTTTGCCGGGAg                                                                                                             >  2:577178/1‑149 (MQ=255)
                                             tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGAACGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGGTTTCCGGGGGTcagc                                                                                                        >  2:89189/1‑150 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagcag                                                                                                     <  2:408599/151‑1 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca                                                                                                       <  2:75396/150‑1 (MQ=255)
                                               gTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAgg                                                                                                    <  1:421047/151‑1 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTACCGCGGGGGCGCATTTTCTCGGCGAGGGAGTTTTTTTGTTTGCCGGGGGGCAGCAggg                                                                                                   >  1:734660/1‑150 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGCGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGGCAGCAggg                                                                                                   >  1:931533/1‑148 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTg                                                                                              <  1:635717/151‑1 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTg                                                                                              <  1:956480/151‑1 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt                                                                                            <  1:218581/150‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGTAAGCATTTTTTTGACGAGGAAGGTTTTTAGTTTGCCGGGAGGCAGCAGGggtttggttt                                                                                          >  2:31068/1‑151 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCCAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTTTCGGGGGTCAGCAGGGCTTTGGgttt                                                                                         >  1:1083465/1‑148 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTc                                                                                         <  1:1014980/149‑1 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACCAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTTGTTTTCCGGGAGGCAGAAGGGGTTTTGGTTTc                                                                                         >  2:270962/1‑150 (MQ=255)
                                                             acgtGACGGTGTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGTGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                       <  1:18123/147‑1 (MQ=255)
                                                               ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGAACGAGGGTACGCATTTTTTCGACGAGGAAGGTTTTCAGTTTGCCGGGGGGCAGGAGGGGTTTGGGTTTCTgg                                                                                     >  2:510881/1‑150 (MQ=255)
                                                                tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTACCGAGGGTGCGCATTTTCTCGACGAGGGAGGTTTTCAGGTTGCCGGGGGTCAGCCGGGGTTTGGTTTTCt                                                                                       >  2:386051/1‑147 (MQ=255)
                                                                  aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTaga                                                                                 <  1:91120/151‑1 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTCCGTTTGCCGGGGGGCAGCAGGGGTTTTGGTTTCTgg                                                                                     >  2:436697/1‑145 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGGGGGGACGCCTTTTTTTGACGAGGGAGTTTTTTCGTTTTCCGGGGGTCAGCAGGGCTTTTGTTTTCTGGTagag                                                                                >  2:451050/1‑149 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTGGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                                                              <  1:612360/151‑1 (MQ=255)
                                                                      ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTTCCGGGGGGACGCCTTTTCTCCGCGAGGAAGTTTTTCATATTGCCCGGGGTCAGCAGGGCttttggtt                                                                                          >  2:1103884/1‑133 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGGAGTGCCGAGGGGGAGCATTTTCTCCGCGAGGGGGTTTTTTAGTTTTCCGGGGGGGAGCAGGGGTTTGGTTTTCCGGTAGGGATACGttt                                                                        >  1:128778/1‑151 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTTGGTTTCTTTTAGAGATACGtttt                                                                       >  2:320551/1‑150 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCGGCAGTTCCCAGGGGAGGGCATTTTTTGGCGGGGGAGTTTTTTAGTTTGTCGGGAGTCAGGAGGGCTTTGGGTTTCTGGTGGAGATTCGtttt                                                                       >  1:460504/1‑150 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttg                                                                      <  1:495321/151‑1 (MQ=255)
                                                                              ttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTGGTTTTCTGGTAGAGATACGtttt                                                                       >  2:328686/1‑149 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                    <  2:69508/150‑1 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTgg                                                                  <  1:83208/151‑1 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                    >  2:155011/1‑149 (MQ=255)
                                                                                    tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCa                                                                <  1:362075/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga                                                              <  2:1094730/151‑1 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCagag                                                             <  1:342445/151‑1 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCGGGGGTCAGCAGGGCTTTGGTTTTCTGGTAGAGATACGTTTTGTTGGCag                                                               >  2:510615/1‑149 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGTACCGGGGGGGCGCAATTTTTCGGCGAGGGAGGTTTTTAGTTTGTCGGGAGGCAGCAGGGGTTTTGGTTTCTTGGAGAGGGACGGTTTGTTGGGag                                                               >  1:712167/1‑147 (MQ=255)
                                                                                          gACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGAGGGGACGCATTTTCTCCACGAGGAAGGTTTTTAGTTTGTCGGGGGGCAGCAGGGGGTTGGTTTTCTGGGAGAGATACGTTTTGTTGGGagag                                                             >  1:38286/1‑147 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGCGGGGACGCATTTTTTTGACGAGGAAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTGGGTTTCTGGTAGGGATACGGTTTTTTTGCCGAggcgcc                                                        >  1:607121/1‑150 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCt                                                    <  2:640645/151‑1 (MQ=255)
                                                                                                ccGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCtt                                                   <  2:991128/151‑1 (MQ=255)
                                                                                                 cGCCAATCACAAACGCGATATGGTACGGCGGGCCAGGTGCAGTAACGAGGGGACGCATTTTTTCTACGGGGAAGTTTTTAAGTTTGTCTGGAGTCAGCAGGGGTTTGGTTTTCTGGGAGGGATACGTTTTTTTTGGAGAgccgccgcc                                                     >  2:400329/1‑148 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGAACCGAGGGGACGCATTTTCTCCACGAGGAAGTTTTTTAATTTGCCGGGGGGCAGCAGGGGTTTGGTTTCCTGGTAGAGATAAGGTTTTTTGGCAGAGGCGCCCCCCTTtgc                                               >  2:675386/1‑151 (MQ=255)
                                                                                                       tCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTTCCGGGGGGCAGCAGGGCTTTGGGTTCCTGGTAGAGATACGTTTTTTTTGCAGAGCCGCCGCCTTTcgcaa                                             >  2:870623/1‑150 (MQ=255)
                                                                                                          caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca                                         <  2:713111/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                        <  2:847276/151‑1 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGTACGCCTTTTCTTGACGAGGGAGGTTTTTAAGTTGGCGGGGGTCAGCAGGGGTTTGGTTTCTTGGTAGAGATAAGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaaccc                                          >  2:71137/1‑146 (MQ=255)
                                                                                                                cgATATGGTACGGCGGGCAGGCCGCAGGACCGAGGGGACGCATTTTCTCGACGAGGGAGGTTTTTTGTTTGCCGGGGGGCAGCAGGGCCTTTGTTTTCTGGTGGAGATACGTTTTGTTGGCGGAGCCGCCGCCTTTCGCCACGCAAAgg                                     >  1:934077/1‑149 (MQ=255)
                                                                                                                               ggCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGcc                    <  2:654961/151‑1 (MQ=255)
                                                                                                                                 cAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                  >  1:403685/1‑151 (MQ=255)
                                                                                                                                   ggCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCt                <  2:412969/151‑1 (MQ=255)
                                                                                                                                   ggCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCt                <  1:343144/151‑1 (MQ=255)
                                                                                                                                            aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc             <  1:679139/145‑1 (MQ=255)
                                                                                                                                            aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc             >  2:679202/1‑145 (MQ=255)
                                                                                                                                            aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTa        <  1:877576/150‑1 (MQ=255)
                                                                                                                                              cGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAg     <  1:451396/151‑1 (MQ=255)
                                                                                                                                              cGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAg     <  1:675323/151‑1 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg    <  2:610571/151‑1 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg    <  2:19797/151‑1 (MQ=255)
                                                                                                                                                aGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg    >  1:424809/1‑150 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  <  2:742705/151‑1 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  <  2:1083564/151‑1 (MQ=255)
                                                                                                                                                 gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  >  1:876166/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  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: