Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG40.9% 111.2 / 10.8 66T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (8/31);  new base G (27/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  NC_000913/4346527‑4346822
                                                                                                                                                     |                                                                                                                                                    
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTAc                                                                                                                                                    <  2:76178/151‑1 (MQ=255)
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTAc                                                                                                                                                    <  2:1139702/151‑1 (MQ=255)
    ccTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTA‑CCGGGGGGACg                                                                                                                                                   >  2:587094/1‑148 (MQ=255)
            tGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGCCGGAGGA‑CCCAGGGGACGCGtttttt                                                                                                                                           >  2:444328/1‑146 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGTA‑CCGAGGGGACGCCTTTTCTcgacg                                                                                                                                      >  1:777944/1‑151 (MQ=255)
               tcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTcgacg                                                                                                                                      <  2:944512/150‑1 (MQ=255)
                 gttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGCA‑CCGAGGGGACGCCATTTTTCGGCG‑Agg                                                                                                                                  >  1:868690/1‑151 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCCCCGTA‑CCGAGGGGACGCATTTTTTTCACG‑AGGaa                                                                                                                                >  1:314816/1‑149 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGt                                                                                                                              <  1:304008/151‑1 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCGTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGtt                                                                                                                             <  2:993638/149‑1 (MQ=255)
                         ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGGAAGA‑CCGAGGGGACGCCTTTTTTCGACG‑AGGAAAttttt                                                                                                                          >  1:1135315/1‑151 (MQ=255)
                                 aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGA‑CCGAGGGGACGCATTTTTTCGCCG‑AGGAAGGTTTTTAGttt                                                                                                                    >  2:472723/1‑149 (MQ=255)
                                   ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGCCAGGGTGCAGCA‑ACCCGGGTAAGCCTTTTTTCCAAC‑AGGAAGattctta                                                                                                                        >  2:1061674/1‑139 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGGCGTA‑CCAGGGGGACGCCTTTTCTCGGCG‑GGGGAGatttttgtgttttcggg                                                                                                              >  1:628399/1‑139 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCgg                                                                                                               <  2:214139/150‑1 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCgg                                                                                                               <  2:679876/150‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCGGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTG‑CCGAGGGTACGCATTTTCTCGTCG‑AGGAAGTTTTTCAGTTTGCCGGGa                                                                                                             <  1:101958/151‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGa                                                                                                             <  2:343465/151‑1 (MQ=255)
                                       tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAATA‑CCGAGGGTAAGCATTTTTTCGACG‑AGGAAGGTTTTTAGTTTGCCggg                                                                                                              >  1:188243/1‑149 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGGACGCCTTTTTTCGACG‑AGGGAGGTTTTCAGTTTGCCggggg                                                                                                            >  1:394975/1‑147 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGGACGCATTTTTTCGACG‑AGGGAGTTTTTTAGTTTGCCggggggc                                                                                                          >  1:1154623/1‑145 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCACGA‑CCGAGGGTACGCAATTTCTCTGCGCGGGAAGTTTTTCGTTTTACCgg                                                                                                               >  2:670232/1‑142 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCCGCAGGA‑ACGAGGGGACGCCTTTTCTCGACG‑AAGGAGGTTTTTAGTTTTGCGGGGGGGAGCAggg                                                                                                   >  2:1059125/1‑150 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCCGCCGGAGGA‑CCGAGGGGACGCATTTTTTTCCCG‑AGGGAGATTTTTCGTTTTCCGGGGGTGAGCAggg                                                                                                   >  1:1017650/1‑150 (MQ=255)
                                                      cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGCCCAGAC‑CCGAGGGTAAGCCTTTTCTCTGCC‑AGGGAGTTTTTTAGTTTGCCGGGGGTGATCTCGGCttt                                                                                               >  2:63187/1‑149 (MQ=255)
                                                       ttGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGCGGGTACGCATTTTTTCGACG‑AGGAAGATTTTTAGTTTGCCGGGGGTCAGCAGGGCCTTgg                                                                                             >  1:983684/1‑150 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCAGCTGAGGTT‑CCCGGGGGGCGCATTTTCTCGCAG‑AGGAAGTGTTTTTGCTTGCCGGGGGACAGCAGGGACTTGGt                                                                                            >  1:770225/1‑150 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt                                                                                            <  1:907401/150‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTTCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt                                                                                            <  2:828104/149‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGtt                                                                                           <  2:106570/150‑1 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                       <  2:480607/151‑1 (MQ=255)
                                                             aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCAGTA‑CCGAGGGGACGCATTTTTTCGGAG‑GGGGAGGTTTTCAGCTTGGCGGGGGGCAGGAGGGGTTTTGTTTcc                                                                                        >  1:159095/1‑149 (MQ=255)
                                                                tGACGGATTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGGA‑CCGAGGGGAGGCATTTTTTTGACG‑GGGGAGGTTTTTAGTTTTCCGGGGGGTAGGAGGGCTTTTGGTTCCt                                                                                       >  2:247238/1‑147 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                  <  1:566214/151‑1 (MQ=255)
                                                                  aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGGACGCATTTTCTCGACG‑AGGGAGGTTTTCAGTTTTCCGGGAGGCAGCAGGGCTTTGGTTTCCTgg                                                                                     >  2:447157/1‑147 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCGGCAAGA‑CCGAGGGGAAGCATTTTCTCGAGG‑AGGAAGAGTTTTCGGTTGCCGCGGGTTGGCCGGGGTTTggggtttt                                                                                       >  1:759379/1‑137 (MQ=255)
                                                                      ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGGGCCATATCCGGGGGGTGGC‑TTTTTTCGGAG‑AGGAAGTTTTTTAGTGTTTCGGGGGTTAGCCGGGCTTTGGGTTCATgg                                                                                     >  2:884614/1‑143 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACg                                                                           <  2:80869/149‑1 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGttt                                                                        <  1:922353/151‑1 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgt                                                                     <  2:261700/150‑1 (MQ=255)
                                                                                  tttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTA‑ACCCGGGGACGCATTTTCTCCACG‑AGGGAGGTTTTTAGTTTGCCCGGGGTCAGCCGGGCTTTGGTTTTCTGGGTGAGATACGTTttgttt                                                                   >  2:1107446/1‑148 (MQ=255)
                                                                                    tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGAT‑CCGAGGGGACGCATTTTCTTTTAG‑AGGGAGTTTTTCAGCTTGCCGGGAGAGAGCAGGGGTTTGGGTTTCTGGGAGAGAGATCTTttgtttg                                                                  >  1:503144/1‑146 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag                                                               <  2:82827/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCa                                                                <  2:705704/149‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGGACGCCTTTTTTCGACG‑AGGAAGTTTTTTAGTTTTCCGGGGGGCCGCAGGGCTTTGGTTTTCTGGTAGAGATACGTtttttt                                                                    >  1:768932/1‑145 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCGGCAGGA‑CCGAGGGGAAGCAATTTTTGCGCG‑AAGAAGGGTTTTTTTTTGCCGGGGGTTTAAAGGGGTTTTGTTTTCTGGGAGAGGTATGTTTTTTTggga                                                                >  2:1176481/1‑146 (MQ=255)
                                                                                       gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCAGGA‑CCGAGGGGACGCATTTTTTCGACG‑AGGGAGATTTTTAGTTTTCCGGGGGTCAGCCGGGGTTTTGTTTTCTGGGAGAGATAAGTTTTTTTGGCaga                                                              >  2:930014/1‑149 (MQ=255)
                                                                                             gTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGGACGCATTTTTTCGACG‑AGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTTGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAggcgccgc                                                      >  1:1136629/1‑151 (MQ=255)
                                                                                                 cGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                  <  1:1099858/151‑1 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTc                                                 <  1:167415/151‑1 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcg                                                <  2:373794/151‑1 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc                                               <  1:877790/151‑1 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGGACGCATTTTTTCGACG‑AGGAAGGTTTTTAGTTTGGCGGGGGTCAGCAGGGCTTTGGTTTTCTGGTAGAGATACGTTTTTTTGGGAGAGCCGCCGCCctt                                                  >  1:505813/1‑148 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCCGCAGTA‑ACGAGGGGACGCATTTTTTCGACG‑AGGGAGGTTTTTAGTTTGGCGGGGGGCAGCAGGGCTTTTGTTTTCTGGTAGAGGTACGTTTTTTTGGGAGAGCCGCCGCCCTTtgcaac                                            >  2:703561/1‑150 (MQ=255)
                                                                                                         acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgc                                          <  2:59912/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                        <  2:977687/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                        <  2:843918/151‑1 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGGTGCAGTA‑ACGAGGGGACGCATTTTTTCGACG‑AGGGAGGTTTTTAGTTTGCCGGGGGGCAGGAGGGGTTTGGGTTTCTGGGAGAGATACGTTTTTTTTGCAGAGGCGCCGCCTTTtgcaacgcaa                                        >  2:1042914/1‑150 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCaaa                                       <  2:1115319/151‑1 (MQ=255)
                                                                                                                      ggTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCGTTCGCAACGCAAAGGAATTTGTa                             <  2:999497/151‑1 (MQ=255)
                                                                                                                           ggcggGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCAt                        <  1:57174/151‑1 (MQ=255)
                                                                                                                             cggGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCg                      >  2:791763/1‑151 (MQ=255)
                                                                                                                             cggGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCg                      <  2:1156391/151‑1 (MQ=255)
                                                                                                                                  aGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                 >  2:104721/1‑151 (MQ=255)
                                                                                                                                      tGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc             >  2:81412/1‑151 (MQ=255)
                                                                                                                                                  gggTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  <  1:1176554/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTA‑CCGAGGGTACGCATTTTCTCGACG‑AGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑