Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG43.3% 157.7 / 27.8 99T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (11/44);  new base G (42/0);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 1.89e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  NC_000913/4346527‑4346825
                                                                                                                                                    |                                                                                                                                                          
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc                                                                                                                                                          <  1:96454/151‑1 (MQ=255)
   gCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCa                                                                                                                                                      <  1:481894/151‑1 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCAttt                                                                                                                                                   <  1:1160170/150‑1 (MQ=255)
        aCCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTtc                                                                                                                                                 <  1:660925/151‑1 (MQ=255)
            tGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGGACCGCGGGGACGCCtttttt                                                                                                                                                >  1:1247934/1‑146 (MQ=255)
            tGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCAGTACCGAGGGGACGCCtttttt                                                                                                                                                >  2:124206/1‑146 (MQ=255)
             gttcgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTcga                                                                                                                                             <  2:1177602/150‑1 (MQ=255)
                cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTACCGCGGGGACGCATTTTTTcgacg                                                                                                                                           >  1:935429/1‑149 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAg                                                                                                                                    <  2:58035/151‑1 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGTGCCGGACCGAGGGGAAGCATTTtcttgcccaggga                                                                                                                                      >  1:151770/1‑139 (MQ=255)
                     ccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGt                                                                                                                                   <  2:559983/151‑1 (MQ=255)
                       ttCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCCGGACCGGGGGGACGCACTTTTTTGGCGAggga                                                                                                                                      >  1:587384/1‑145 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGGAGtttt                                                                                                                                >  2:914050/1‑151 (MQ=255)
                        tCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGGAGgttt                                                                                                                                >  1:72434/1‑151 (MQ=255)
                         ccGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGtttt                                                                                                                                <  2:399501/150‑1 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTc                                                                                                                              <  1:264528/151‑1 (MQ=255)
                           gTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCAGTAACGAGGGGACGCATTTTTTCGACGAGGGAGtttttt                                                                                                                              >  2:1083728/1‑149 (MQ=255)
                             cGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGGGCAGGACCGAGGGGACGGAATTTTTCGATGAGGAAGgttttttgt                                                                                                                           >  1:705907/1‑147 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtt                                                                                                                          <  2:163581/151‑1 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCCGCCCGACCGGGGGGACGCCTTTTTTCCACGGGGGAGgtttt                                                                                                                               >  1:988373/1‑146 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTGTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTg                                                                                                                        <  2:1008019/151‑1 (MQ=255)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGTGCAGTACCGAGGGGCCGCATTTTTTCGGCGAGGAAGGTTTTTAGttt                                                                                                                         >  1:882068/1‑150 (MQ=255)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTg                                                                                                                        <  2:255669/151‑1 (MQ=255)
                                cAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTg                                                                                                                        <  1:161655/151‑1 (MQ=255)
                                  gTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc                                                                                                                      <  2:755967/151‑1 (MQ=255)
                                   ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCAGGACCGCGGGGACGCATTTTTTCGACGAGGGAGATTTTTAGttttgc                                                                                                                      >  1:988644/1‑147 (MQ=255)
                                   ttCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGAAGGTTTTCAGTTTGcc                                                                                                                      >  2:911539/1‑150 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCgg                                                                                                                    <  2:1116664/150‑1 (MQ=255)
                                          tAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGTAACGAGGGGAGGACATTTTTCCACGCGGAAGGTTTTTCGGTTTCCCg                                                                                                                     >  1:707091/1‑145 (MQ=255)
                                            aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCCAGGGGACGCCTTTTTTTTGCGAGGGAGGTTTTTAGTTTGCCggggg                                                                                                                 >  2:951645/1‑144 (MQ=255)
                                              aGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGGGGGGACGCATTTTTTCGGCGAGGGAGGTTTTTCGTTTTCCggggggca                                                                                                              >  1:823021/1‑142 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCtt                                                                                                   <  1:576988/151‑1 (MQ=255)
                                                   gcgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCtt                                                                                                   <  2:430795/151‑1 (MQ=255)
                                                    cgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTTCCGGGGGGCAGCAGGGCttt                                                                                                  >  2:1175740/1‑151 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGTAACGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTTAGTTTTGCGGGGGGCAGCAggggctt                                                                                                   >  2:417850/1‑146 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTg                                                                                                 <  2:371369/151‑1 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTg                                                                                                 <  2:462491/151‑1 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGGCGGGGGGCAGCAGGGGTTTgggt                                                                                              >  2:96335/1‑149 (MQ=255)
                                                          cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGCCGGTGCCGGGGGGACGCCTTTTTTCCACGAGGAAGGTTTTTCATTTGCCGGGGGGGAGGAGGGCTTTTGGttt                                                                                             >  1:505046/1‑151 (MQ=255)
                                                            aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt                                                                                          <  2:437465/151‑1 (MQ=255)
                                                             aCTTGACGGTTTTCAGTTTGGTTTCCGCAGACGTACCGCCAATCTCAAACGCTATATGGTACGGCGGGCAGGCTGGAGTACCGAGGCTACGCATTTTCTCGACCAGGAAGTTTTTCAGTTTGCCGGGAGTGAGTAGGGCTTTGGtt                                                                                              <  1:368404/146‑1 (MQ=255)
                                                             aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTg                                                                                         <  1:451850/151‑1 (MQ=255)
                                                               ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTGGTTTCCTgg                                                                                        >  2:301032/1‑150 (MQ=255)
                                                               ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCCAGGGGACGGATTTTTTTCCCGAGGAAGTTTTTTTGTTTGCCGGGGGGCAGCAGGGCTTTTGttttc                                                                                           >  1:782747/1‑145 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCACTCCCAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTa                                                                                      <  2:553441/150‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                     <  1:553698/151‑1 (MQ=255)
                                                                 gACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag                                                                                     <  2:1277379/151‑1 (MQ=255)
                                                                    gggtttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGAGCGAGGGCACGCATTTTCTTTAAGAGGAAGGTTTTCATTTTGCCGGGGGACAGCAGGGGCTTGGTTTCCTg                                                                                         >  2:1206737/4‑144 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTAACGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGAGGCAGCAGGGCTTTGGTTTTCTGGGagag                                                                                   >  2:16161/1‑150 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                                                                 <  1:436526/151‑1 (MQ=255)
                                                                      ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCAGCCGCAGGAACGGGGGGACGCCATTTTTCGACGGGGGAGGGTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTTGTTTTCTGGTAgaga                                                                                  >  2:802939/1‑149 (MQ=255)
                                                                        ttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGCCAGGGTCCACTACCGAGGGGACGCAATTTTTCGAAGAGGAAGTTTTTCAGTTTTGCGGGGGGCAGCCAGGGTTTGGTTTCCTCGTAGAGATACg                                                                              >  2:1253515/1‑151 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACg                                                                              <  1:559981/150‑1 (MQ=255)
                                                                          cAGGTTGGTTTCCGCGGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGGAGTTTTTTCGTTTGCCGGGGGGCAGCCGGGCGTTGGGTTCCTGGGAGAGATACGtt                                                                            >  1:354002/1‑151 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGGGGGGACGCATTTTCTCGACGAGGGAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGCGTTGGGTTTCTGGTAGAGATACGtt                                                                            >  2:1031176/1‑151 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGt                                                                             <  2:121167/150‑1 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCCCGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGGGAGCAGGGCTTTTGTTTTCTGGTGGAGATAAGtt                                                                            >  1:756742/1‑151 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGACGCATTTTTTCGGCGAGGGAGTTTTTTAGTTTGCCGGGAGTCAGCAGGGGTTTGGTTTTCTGGTAGAGATACGgtt                                                                           >  1:573023/1‑151 (MQ=255)
                                                                           aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCCGCAGTACCGGGGGGACGCATTTTCTCGACGAGGGAGTTTTTTAGTTTGCCGGGAGGCAGCAGGGGCTTGGTTTCCTGGGAGAGATACGttt                                                                           >  1:546447/1‑151 (MQ=255)
                                                                            ggttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGGAGGAACGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGGCGGGGGGCAGCAGGGCTTTGGTTTCCTGGGAGAGATACGtttt                                                                          >  2:1242773/1‑151 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTCAGTTTGCCGGGGGGCAGCAGGGCTTTTGTTTCCTGGGTGAGATACGtttt                                                                          >  2:639191/1‑150 (MQ=255)
                                                                             gttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAAGCCGCTGCAACGGGGGGACGGAATTTCCATTCGGGGGAGTTTTTTTGATTGCCGGGAGGCAGTAGGGGCTTGGGTTTCTTggagtggtaagtt                                                                             >  2:1090695/1‑137 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCCGGACCGGGGGGACGCCTTTTTTCGGCGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTGGGTTTCTGGTAGAGAAACGGtttttt                                                                       >  2:525484/1‑151 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgt                                                                        <  1:426965/150‑1 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgtt                                                                       <  1:726054/150‑1 (MQ=255)
                                                                                   ttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTTCCGGGGGGCAGCAGGGGTTTGGGTTCCTGGTAGAGATACGTTttgtttg                                                                     >  2:1094208/1‑147 (MQ=255)
                                                                                      cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGAGGCAGCAGGGCTTTGGTTTTCTTGGAGAGATACGTTTTTTTGGCag                                                                  >  2:75764/1‑149 (MQ=255)
                                                                                          gACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGGACCGAGGGGACGCCTTTTCTCGACGAGGAAGGTTTTCAGTTTGCCGGGAGGCAGCAGGGCTTTGGTTTTCTGGTAGAGATACGTTTTTTTGGGAGAgccgc                                                            >  2:791098/1‑151 (MQ=255)
                                                                                           aCGTACCGCCAATCACAAACGCGATATGGTACGGCGGGAAGGGTGCAGGACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAAGTTTTGTTGGGAGAgccgcc                                                           >  2:1018865/1‑151 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgccgcc                                                           <  2:746412/150‑1 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGGGAAGCTGCAGGACCCAGGGGACGCCTTTTCTCTACGAGGAAGGTTTTTAGTTTTGCGGGGGGGAGCAGGGCTTTGGTTTTCTGGGAGAGATACGGTTTGTTGGCCGAgccgccgcc                                                        >  1:770330/1‑151 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGCTGCAGTACCGAGGGGGCGCCTTTTCTCGACGAGGAAGTTTTTTAGTTTGCCGGGAGGCAGCAGGGCTTTGGGTTTCTGGTAGAGGTAAGTTTTGTTTGCAGAgccgccgcc                                                        >  1:1107344/1‑151 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGTTTTTCAGTTTGCCGGGGGTCAGCAGGGCTTTTGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgccgccgcc                                                        >  2:24398/1‑151 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGAAGCATTTTCTCGACGAGGAAGGTTTTTAGTTTTCCGGGAGTCAGCAGGGGTTTGGTTTCCTGGTAGAGATACGTTTTTTTGGCAGAgccgccgcc                                                        >  2:508418/1‑151 (MQ=255)
                                                                                                ccGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCtt                                                      <  2:648704/151‑1 (MQ=255)
                                                                                                  gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTc                                                    <  2:587483/151‑1 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcg                                                   <  1:278706/151‑1 (MQ=255)
                                                                                                    cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgc                                                  <  1:104076/151‑1 (MQ=255)
                                                                                                      aTCACAAACGCGATATGGTACGGCGGGCAGGGGGGAAGACCGAGGGAACGCATTTTTTCGAAGAGGAAGGTTTTCCGGTTGCCGGGGGTCAGGAGGGGCTTTGTTTTCTTGTAGAGGTACGTTTTGTTGGCAGGGCCGGCGCCCTTccta                                                 >  2:129712/1‑147 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaac                                               <  2:375807/150‑1 (MQ=255)
                                                                                                                cgATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTTAGTTTTCCGGGGGGCAGCAGGGCTTTGGGTTTCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCTCAACGCCAAgg                                        >  1:1032855/1‑149 (MQ=255)
                                                                                                                cgATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTTTGTTTGCCGGGGGTCAGCAGGGCTTTGGGTTCCTGGGAGAGATACGGTTTGTTGGGCGAGCCGCCGCCCTTTGCAACGCAAAgg                                        >  1:199482/1‑149 (MQ=255)
                                                                                                                      ggTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGt                                 <  1:442998/150‑1 (MQ=255)
                                                                                                                       gTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTa                                >  1:390724/1‑150 (MQ=255)
                                                                                                                         aCGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTc                             <  2:28346/151‑1 (MQ=255)
                                                                                                                          cggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTc                             >  2:526194/1‑150 (MQ=255)
                                                                                                                           ggcggGCAGGCTGCAGTACCGCGGGTACGCACTTTCTCGACGAGGAAGTTTTTAAGTTTGCCGGGAGCCAGCCGGGCTTTGGTTTCCTGGTAGAGATTCGTTTTGTCGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGCATTTGTACTc                             >  1:938216/1‑149 (MQ=255)
                                                                                                                           ggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCAt                           >  1:21145/1‑151 (MQ=255)
                                                                                                                              gggCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTTAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGc                        >  1:1103756/1‑151 (MQ=255)
                                                                                                                              gggCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCg                         <  1:1087203/150‑1 (MQ=255)
                                                                                                                                gCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCa                      <  2:78012/151‑1 (MQ=255)
                                                                                                                                  aGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                    >  2:1053840/1‑151 (MQ=255)
                                                                                                                                       gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAc                 <  2:1070268/149‑1 (MQ=255)
                                                                                                                                       gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcg               <  2:905416/151‑1 (MQ=255)
                                                                                                                                        cAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgc              >  2:513233/1‑151 (MQ=255)
                                                                                                                                           tACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGt            >  2:936035/1‑150 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       <  1:1214085/151‑1 (MQ=255)
                                                                                                                                               gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg       >  1:336735/1‑151 (MQ=255)
                                                                                                                                                   gTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGa   <  2:768380/151‑1 (MQ=255)
                                                                                                                                                    tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAt  >  2:1041574/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                          
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  NC_000913/4346527‑4346825

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: