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

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTC  >  NC_000913/4346527‑4346822
                                                                                                                                                     |                                                                                                                                                    
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                    <  2:76178/151‑1 (MQ=255)
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                    <  2:1139702/151‑1 (MQ=255)
    ccTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGGGGGGACg                                                                                                                                                   >  2:587094/1‑148 (MQ=255)
            tGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGCCGGAGGACCCAGGGGACGCGtttttt                                                                                                                                           >  2:444328/1‑146 (MQ=255)
              ttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGTACCGAGGGGACGCCTTTTCTcgacg                                                                                                                                      >  1:777944/1‑151 (MQ=255)
               tcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcgacg                                                                                                                                      <  2:944512/150‑1 (MQ=255)
                 gttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGCACCGAGGGGACGCCATTTTTCGGCGAgg                                                                                                                                  >  1:868690/1‑151 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCCCCGTACCGAGGGGACGCATTTTTTTCACGAGGaa                                                                                                                                >  1:314816/1‑149 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGt                                                                                                                              <  1:304008/151‑1 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCGTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGtt                                                                                                                             <  2:993638/149‑1 (MQ=255)
                         ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGGAAGACCGAGGGGACGCCTTTTTTCGACGAGGAAAttttt                                                                                                                          >  1:1135315/1‑151 (MQ=255)
                                 aGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGACCGAGGGGACGCATTTTTTCGCCGAGGAAGGTTTTTAGttt                                                                                                                    >  2:472723/1‑149 (MQ=255)
                                   ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGCCAGGGTGCAGCAACCCGGGTAAGCCTTTTTTCCAACAGGAAGattctta                                                                                                                        >  2:1061674/1‑139 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGGCGTACCAGGGGGACGCCTTTTCTCGGCGGGGGAGatttttgtgttttcggg                                                                                                              >  1:628399/1‑139 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCgg                                                                                                               <  2:214139/150‑1 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCgg                                                                                                               <  2:679876/150‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCGGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTGCCGAGGGTACGCATTTTCTCGTCGAGGAAGTTTTTCAGTTTGCCGGGa                                                                                                             <  1:101958/151‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGa                                                                                                             <  2:343465/151‑1 (MQ=255)
                                       tCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAATACCGAGGGTAAGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCggg                                                                                                              >  1:188243/1‑149 (MQ=255)
                                         gTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTCAGTTTGCCggggg                                                                                                            >  1:394975/1‑147 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCggggggc                                                                                                          >  1:1154623/1‑145 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCACGACCGAGGGTACGCAATTTCTCTGCGCGGGAAGTTTTTCGTTTTACCgg                                                                                                               >  2:670232/1‑142 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCCGCAGGAACGAGGGGACGCCTTTTCTCGACGAAGGAGGTTTTTAGTTTTGCGGGGGGGAGCAggg                                                                                                   >  2:1059125/1‑150 (MQ=255)
                                                 gAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCCGCCGGAGGACCGAGGGGACGCATTTTTTTCCCGAGGGAGATTTTTCGTTTTCCGGGGGTGAGCAggg                                                                                                   >  1:1017650/1‑150 (MQ=255)
                                                      cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGCCCAGACCCGAGGGTAAGCCTTTTCTCTGCCAGGGAGTTTTTTAGTTTGCCGGGGGTGATCTCGGCttt                                                                                               >  2:63187/1‑149 (MQ=255)
                                                       ttGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGCGGGTACGCATTTTTTCGACGAGGAAGATTTTTAGTTTGCCGGGGGTCAGCAGGGCCTTgg                                                                                             >  1:983684/1‑150 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCAGCTGAGGTTCCCGGGGGGCGCATTTTCTCGCAGAGGAAGTGTTTTTGCTTGCCGGGGGACAGCAGGGACTTGGt                                                                                            >  1:770225/1‑150 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt                                                                                            <  1:907401/150‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTTCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt                                                                                            <  2:828104/149‑1 (MQ=255)
                                                         gCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGtt                                                                                           <  2:106570/150‑1 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                       <  2:480607/151‑1 (MQ=255)
                                                             aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCAGTACCGAGGGGACGCATTTTTTCGGAGGGGGAGGTTTTCAGCTTGGCGGGGGGCAGGAGGGGTTTTGTTTcc                                                                                        >  1:159095/1‑149 (MQ=255)
                                                                tGACGGATTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGGACCGAGGGGAGGCATTTTTTTGACGGGGGAGGTTTTTAGTTTTCCGGGGGGTAGGAGGGCTTTTGGTTCCt                                                                                       >  2:247238/1‑147 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                  <  1:566214/151‑1 (MQ=255)
                                                                  aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGGTTTTCAGTTTTCCGGGAGGCAGCAGGGCTTTGGTTTCCTgg                                                                                     >  2:447157/1‑147 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCGGCAAGACCGAGGGGAAGCATTTTCTCGAGGAGGAAGAGTTTTCGGTTGCCGCGGGTTGGCCGGGGTTTggggtttt                                                                                       >  1:759379/1‑137 (MQ=255)
                                                                      ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGGGCCATATCCGGGGGGTGGCTTTTTTCGGAGAGGAAGTTTTTTAGTGTTTCGGGGGTTAGCCGGGCTTTGGGTTCATgg                                                                                     >  2:884614/1‑143 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACg                                                                           <  2:80869/149‑1 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGttt                                                                        <  1:922353/151‑1 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgt                                                                     <  2:261700/150‑1 (MQ=255)
                                                                                  tttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTAACCCGGGGACGCATTTTCTCCACGAGGGAGGTTTTTAGTTTGCCCGGGGTCAGCCGGGCTTTGGTTTTCTGGGTGAGATACGTTttgttt                                                                   >  2:1107446/1‑148 (MQ=255)
                                                                                    tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGATCCGAGGGGACGCATTTTCTTTTAGAGGGAGTTTTTCAGCTTGCCGGGAGAGAGCAGGGGTTTGGGTTTCTGGGAGAGAGATCTTttgtttg                                                                  >  1:503144/1‑146 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag                                                               <  2:82827/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCa                                                                <  2:705704/149‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGAAGTTTTTTAGTTTTCCGGGGGGCCGCAGGGCTTTGGTTTTCTGGTAGAGATACGTtttttt                                                                    >  1:768932/1‑145 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGCGGCAGGACCGAGGGGAAGCAATTTTTGCGCGAAGAAGGGTTTTTTTTTGCCGGGGGTTTAAAGGGGTTTTGTTTTCTGGGAGAGGTATGTTTTTTTggga                                                                >  2:1176481/1‑146 (MQ=255)
                                                                                       gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGATTTTTAGTTTTCCGGGGGTCAGCCGGGGTTTTGTTTTCTGGGAGAGATAAGTTTTTTTGGCaga                                                              >  2:930014/1‑149 (MQ=255)
                                                                                             gTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTTGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAggcgccgc                                                      >  1:1136629/1‑151 (MQ=255)
                                                                                                 cGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                  <  1:1099858/151‑1 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTc                                                 <  1:167415/151‑1 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcg                                                <  2:373794/151‑1 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc                                               <  1:877790/151‑1 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGGCGGGGGTCAGCAGGGCTTTGGTTTTCTGGTAGAGATACGTTTTTTTGGGAGAGCCGCCGCCctt                                                  >  1:505813/1‑148 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCCGCAGTAACGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGGCGGGGGGCAGCAGGGCTTTTGTTTTCTGGTAGAGGTACGTTTTTTTGGGAGAGCCGCCGCCCTTtgcaac                                            >  2:703561/1‑150 (MQ=255)
                                                                                                         acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgc                                          <  2:59912/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                        <  2:977687/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgcaa                                        <  2:843918/151‑1 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGGTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGGCAGGAGGGGTTTGGGTTTCTGGGAGAGATACGTTTTTTTTGCAGAGGCGCCGCCTTTtgcaacgcaa                                        >  2:1042914/1‑150 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCaaa                                       <  2:1115319/151‑1 (MQ=255)
                                                                                                                      ggTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCGTTCGCAACGCAAAGGAATTTGTa                             <  2:999497/151‑1 (MQ=255)
                                                                                                                           ggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCAt                        <  1:57174/151‑1 (MQ=255)
                                                                                                                             cggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCg                      >  2:791763/1‑151 (MQ=255)
                                                                                                                             cggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCg                      <  2:1156391/151‑1 (MQ=255)
                                                                                                                                  aGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                 >  2:104721/1‑151 (MQ=255)
                                                                                                                                      tGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcc             >  2:81412/1‑151 (MQ=255)
                                                                                                                                                  gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc  <  1:1176554/151‑1 (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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: