Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,073,741 A→T 26.9% L283Q (CTG→CAG)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,073,7410AT26.9% 360.7 / 155.2 200L283Q (CTG→CAG) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (69/77);  new base T (21/33);  total (90/110)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 210 total aligned reads displayed.

TGGTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTC  >  CP009273/4073595‑4073887
                                                                                                                                                  |                                                                                                                                                  
tGGTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCag                                                                                                                                                   <  2:1628588/148‑1 (MQ=255)
 ggTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCaga                                                                                                                                                  >  1:2185680/1‑148 (MQ=255)
   tgtCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCagaga                                                                                                                                                <  1:1488696/148‑1 (MQ=255)
   tgtCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCagaga                                                                                                                                                <  2:366757/148‑1 (MQ=255)
       gtcACACCCTATGAGCTGTACTCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCg                                                                                                                                            <  1:631162/148‑1 (MQ=255)
       gtcACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCg                                                                                                                                            <  2:3028867/148‑1 (MQ=255)
        tcACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGa                                                                                                                                           >  2:3348608/1‑148 (MQ=255)
         cacaCCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAg                                                                                                                                          >  2:2264592/1‑148 (MQ=255)
         cacaCCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAg                                                                                                                                          >  2:471207/1‑148 (MQ=255)
         cacaCCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAg                                                                                                                                          >  2:2062122/1‑148 (MQ=255)
           cacCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGt                                                                                                                                        >  2:3222786/1‑148 (MQ=255)
            acCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGAGTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTc                                                                                                                                       <  2:2886206/148‑1 (MQ=255)
             cccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCa                                                                                                                                      <  2:1650778/148‑1 (MQ=255)
              ccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGAGTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCAt                                                                                                                                     <  2:923295/148‑1 (MQ=255)
                   gAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCt                                                                                                                                >  2:1748919/1‑148 (MQ=255)
                   gAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCt                                                                                                                                >  2:1748947/1‑148 (MQ=255)
                      cTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACg                                                                                                                             >  1:24143/1‑148 (MQ=255)
                       tGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACgg                                                                                                                            <  2:113559/148‑1 (MQ=255)
                       tGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACgg                                                                                                                            <  1:360729/148‑1 (MQ=255)
                       tGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACgg                                                                                                                            <  2:1193697/148‑1 (MQ=255)
                        gTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGc                                                                                                                           >  1:1697997/1‑148 (MQ=255)
                        gTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGc                                                                                                                           >  1:2799091/1‑148 (MQ=255)
                        gTAATCCGGGCGACCTGAGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGc                                                                                                                           <  1:2369511/148‑1 (MQ=255)
                          aaTCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCtt                                                                                                                         <  1:1097700/148‑1 (MQ=255)
                          aaTCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCtt                                                                                                                         <  1:3430652/148‑1 (MQ=255)
                            tCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCg                                                                                                                       <  2:3131876/148‑1 (MQ=255)
                              cGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAt                                                                                                                     <  1:2852423/148‑1 (MQ=255)
                                ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAttt                                                                                                                   >  1:1783572/1‑148 (MQ=255)
                                 gCGACCTGCGCCACACGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGCTTGATGGAACTGAGGGCATAGCCTTCTTGCTCCATACGTGCGTCCAGTACCAGAGAGGCGAGCTAATATGCTACGGCTTCGaagta                                                                                                                  >  2:3281440/1‑144 (MQ=255)
                                  cGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAttttt                                                                                                                 <  2:2977557/148‑1 (MQ=255)
                                       tGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATc                                                                                                            <  2:1995828/148‑1 (MQ=255)
                                        gcgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCt                                                                                                           <  1:2889716/148‑1 (MQ=255)
                                        gcgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCt                                                                                                           <  1:177456/148‑1 (MQ=255)
                                         cgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtt                                                                                                          <  2:1465398/148‑1 (MQ=255)
                                         cgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtt                                                                                                          <  1:2718097/148‑1 (MQ=255)
                                         cgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtt                                                                                                          <  2:1761392/148‑1 (MQ=255)
                                         cgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtt                                                                                                          <  2:272884/148‑1 (MQ=255)
                                          gcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCttt                                                                                                         <  1:1341995/148‑1 (MQ=255)
                                             aCCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTc                                                                                                      <  1:3091682/148‑1 (MQ=255)
                                              cccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     >  1:319031/1‑148 (MQ=255)
                                              cccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     <  2:2442450/148‑1 (MQ=255)
                                              cccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     <  1:28725/148‑1 (MQ=255)
                                              cccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     <  1:3534357/148‑1 (MQ=255)
                                                cGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGTACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTg                                                                                                   <  1:322915/148‑1 (MQ=255)
                                                 ggAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGa                                                                                                  <  1:585801/148‑1 (MQ=255)
                                                 ggAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGa                                                                                                  <  2:1546224/148‑1 (MQ=255)
                                                 ggAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGa                                                                                                  <  2:1381500/148‑1 (MQ=255)
                                                  gAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGAt                                                                                                 <  1:1795594/148‑1 (MQ=255)
                                                  gAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGAt                                                                                                 >  2:590215/1‑148 (MQ=255)
                                                       cGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAgaga                                                                                            <  2:1538579/148‑1 (MQ=255)
                                                       cGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAgaga                                                                                            <  1:1044163/148‑1 (MQ=255)
                                                         tCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAt                                                                                          <  2:702580/148‑1 (MQ=255)
                                                            gagaTTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAtttt                                                                                       <  1:2282985/148‑1 (MQ=255)
                                                            gagaTTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAtttt                                                                                       >  1:1091411/1‑148 (MQ=255)
                                                                ttACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAgg                                                                                   >  1:1246830/1‑148 (MQ=255)
                                                                 tACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGt                                                                                  >  2:1235644/1‑148 (MQ=255)
                                                                    tctcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAg                                                                               >  1:275404/1‑148 (MQ=255)
                                                                    tctcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAg                                                                               >  2:804149/1‑148 (MQ=255)
                                                                       ccTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAgtgt                                                                            <  2:932478/148‑1 (MQ=255)
                                                                          tCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGc                                                                         >  1:1286357/1‑148 (MQ=255)
                                                                            aCCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAg                                                                       <  2:847289/148‑1 (MQ=255)
                                                                            aCCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAg                                                                       >  2:530287/1‑148 (MQ=255)
                                                                             cccGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGt                                                                      >  1:3466301/1‑148 (MQ=255)
                                                                                ggAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCAt                                                                   >  1:729616/1‑148 (MQ=255)
                                                                                 gAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCTGTCATc                                                                  >  2:1832304/1‑148 (MQ=255)
                                                                                 gAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATc                                                                  >  2:3006580/1‑148 (MQ=255)
                                                                                  aaaCAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCg                                                                 <  1:3348608/148‑1 (MQ=255)
                                                                                     cAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAg                                                              >  1:2247498/1‑148 (MQ=255)
                                                                                     cAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAg                                                              >  1:68693/1‑148 (MQ=255)
                                                                                     cAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAg                                                              >  2:2164104/1‑148 (MQ=255)
                                                                                      aGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCNTCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGc                                                             <  1:1808542/148‑1 (MQ=255)
                                                                                       ggAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCt                                                            >  1:1991092/1‑148 (MQ=255)
                                                                                        gAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCtt                                                           >  1:1694824/1‑148 (MQ=255)
                                                                                        gAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCtt                                                           <  2:34117/148‑1 (MQ=255)
                                                                                             ggTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGc                                                      <  1:961761/148‑1 (MQ=255)
                                                                                              gTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCt                                                     <  1:2744719/148‑1 (MQ=255)
                                                                                               ttGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCtt                                                    <  1:2333993/148‑1 (MQ=255)
                                                                                                tGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCttt                                                   >  1:1083571/1‑148 (MQ=255)
                                                                                                 gATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTa                                                  <  1:276477/148‑1 (MQ=255)
                                                                                                     gAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcg                                              >  1:2649162/1‑148 (MQ=255)
                                                                                                        cTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCgg                                           >  1:1939080/1‑148 (MQ=255)
                                                                                                              gCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGtt                                     <  1:1663139/148‑1 (MQ=255)
                                                                                                                    ccTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCAt                               >  1:3533725/1‑148 (MQ=255)
                                                                                                                     cttcttGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATc                              <  2:1246830/148‑1 (MQ=255)
                                                                                                                          tGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCg                         >  2:2032298/1‑148 (MQ=255)
                                                                                                                           gCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGa                        <  2:2471973/148‑1 (MQ=255)
                                                                                                                            cTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAc                       <  2:665863/148‑1 (MQ=255)
                                                                                                                            cTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAc                       <  1:1748947/148‑1 (MQ=255)
                                                                                                                            cTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAc                       <  1:1748919/148‑1 (MQ=255)
                                                                                                                                aTTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGt                   <  2:3466301/148‑1 (MQ=255)
                                                                                                                                 ttCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTTCGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGGTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  <  1:111746/148‑1 (MQ=255)
                                                                                                                                 ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  >  1:3040093/1‑148 (MQ=255)
                                                                                                                                 ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  >  1:227105/1‑148 (MQ=255)
                                                                                                                                 ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  >  1:1259062/1‑148 (MQ=255)
                                                                                                                                   cGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAAt                <  1:1644674/148‑1 (MQ=255)
                                                                                                                                     tGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGc              <  2:275285/148‑1 (MQ=255)
                                                                                                                                       cGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAg            <  1:2495541/148‑1 (MQ=255)
                                                                                                                                       cGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAg            <  1:708308/148‑1 (MQ=255)
                                                                                                                                             gTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGcc      >  1:2660198/1‑148 (MQ=255)
                                                                                                                                                ccAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGct   <  2:1919285/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGTGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTC  >  CP009273/4073595‑4073887

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: