Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,073,741 A→T 43.8% 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,7410AT43.8% 101.7 / 287.5 217L283Q (CTG→CAG) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (59/63);  new base T (38/57);  total (97/120)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

TGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTC  >  CP009273/4073598‑4073887
                                                                                                                                               |                                                                                                                                                  
tgtCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCagaga                                                                                                                                                <  1:719986/148‑1 (MQ=255)
      cacaCCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAg                                                                                                                                          >  1:3049889/1‑148 (MQ=255)
      cacaCCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAg                                                                                                                                          >  1:1274891/1‑148 (MQ=255)
       acacCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAgg                                                                                                                                         >  1:1601218/1‑148 (MQ=255)
        cacCCTTTGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGt                                                                                                                                        >  2:1281905/1‑148 (MQ=255)
         acCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTc                                                                                                                                       <  2:2720171/148‑1 (MQ=255)
          cccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCa                                                                                                                                      >  1:651665/1‑148 (MQ=255)
           ccTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCAt                                                                                                                                     <  2:1008247/148‑1 (MQ=255)
              aTGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTg                                                                                                                                  <  1:292441/148‑1 (MQ=255)
              aTGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTg                                                                                                                                  <  2:353645/148‑1 (MQ=255)
               tGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGc                                                                                                                                 >  1:1660915/1‑148 (MQ=255)
                   cTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACg                                                                                                                             >  1:1169873/1‑148 (MQ=255)
                      tAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCt                                                                                                                          <  1:608355/148‑1 (MQ=255)
                        aTCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTc                                                                                                                        >  1:2979062/1‑148 (MQ=255)
                          ccGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGa                                                                                                                      <  2:2039030/148‑1 (MQ=255)
                          ccGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGa                                                                                                                      <  1:447685/148‑1 (MQ=255)
                           cGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAt                                                                                                                     <  1:1860640/148‑1 (MQ=255)
                             ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGAtgt                                                                                                                   >  2:850607/1‑146 (MQ=255)
                             ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAttt                                                                                                                   <  1:866041/148‑1 (MQ=255)
                             ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAttt                                                                                                                   <  1:3076197/148‑1 (MQ=255)
                             ggCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGGTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGAttt                                                                                                                   <  1:937508/148‑1 (MQ=255)
                              gCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGAtttt                                                                                                                  >  2:1938837/1‑148 (MQ=255)
                              gCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAtttt                                                                                                                  <  2:771762/148‑1 (MQ=255)
                              gCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGAtttt                                                                                                                  <  2:1179418/148‑1 (MQ=255)
                                 aCCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGATCATCTGCTACGGCTTCGATTTTTgg                                                                                                               <  1:247320/148‑1 (MQ=255)
                                  ccTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGa                                                                                                              <  1:510161/148‑1 (MQ=255)
                                   cTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGAt                                                                                                             >  2:1415650/1‑148 (MQ=255)
                                   cTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGAt                                                                                                             >  2:2795856/1‑148 (MQ=255)
                                    tGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATc                                                                                                            <  2:2713185/148‑1 (MQ=255)
                                     gcgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCt                                                                                                           <  2:1873864/148‑1 (MQ=255)
                                      cgcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtt                                                                                                          <  1:180074/148‑1 (MQ=255)
                                       gcCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCttt                                                                                                         <  2:922217/148‑1 (MQ=255)
                                        ccACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtttt                                                                                                        >  1:2039106/1‑148 (MQ=255)
                                        ccACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCtttt                                                                                                        >  1:2746536/1‑148 (MQ=255)
                                           cccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTcc                                                                                                     <  1:3492595/148‑1 (MQ=255)
                                            ccGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCt                                                                                                    >  1:2920339/1‑148 (MQ=255)
                                              ggAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGa                                                                                                  <  2:620960/148‑1 (MQ=255)
                                                   aCGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATagag                                                                                             >  2:445081/1‑148 (MQ=255)
                                                      tCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAt                                                                                          <  2:214618/148‑1 (MQ=255)
                                                       cAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAtt                                                                                         <  2:1048735/148‑1 (MQ=255)
                                                             ttACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAgg                                                                                   >  1:1465599/1‑148 (MQ=255)
                                                                 tctcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAg                                                                               <  1:1666637/148‑1 (MQ=255)
                                                                  ctcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAgt                                                                              <  1:302607/148‑1 (MQ=255)
                                                                   tcCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTTCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAgtg                                                                             <  1:3447188/148‑1 (MQ=255)
                                                                      ttCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCg                                                                          <  1:2164594/148‑1 (MQ=255)
                                                                        cACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCa                                                                        <  1:826367/148‑1 (MQ=255)
                                                                          cccGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGt                                                                      >  1:1361530/1‑148 (MQ=255)
                                                                           ccGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTc                                                                     >  2:2508145/1‑148 (MQ=255)
                                                                            cGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCa                                                                    <  1:2725450/148‑1 (MQ=255)
                                                                            cGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCa                                                                    <  1:1144588/148‑1 (MQ=255)
                                                                                  cAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAg                                                              <  1:699672/148‑1 (MQ=255)
                                                                                   aGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGc                                                             <  2:1601218/148‑1 (MQ=255)
                                                                                   aGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGc                                                             <  2:168276/148‑1 (MQ=255)
                                                                                      aaCGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCttt                                                          <  1:128804/148‑1 (MQ=255)
                                                                                       aCGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTc                                                         >  2:72059/1‑148 (MQ=255)
                                                                                        cGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCg                                                        >  2:1374661/1‑148 (MQ=255)
                                                                                        cGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCg                                                        >  1:353593/1‑148 (MQ=255)
                                                                                          ggTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGc                                                      <  2:2479171/148‑1 (MQ=255)
                                                                                           gTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCt                                                     >  1:2619000/1‑148 (MQ=255)
                                                                                           gTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCt                                                     >  1:2619174/1‑148 (MQ=255)
                                                                                           gTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCt                                                     >  1:223121/1‑148 (MQ=255)
                                                                                               aTGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTaa                                                 >  1:1672221/1‑148 (MQ=255)
                                                                                               aTGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTaa                                                 >  1:2249710/1‑148 (MQ=255)
                                                                                               aTGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTaa                                                 <  1:3272663/148‑1 (MQ=255)
                                                                                               aTGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTaa                                                 >  1:117656/1‑148 (MQ=255)
                                                                                                tGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAAt                                                >  2:2737314/1‑148 (MQ=255)
                                                                                                tGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAAt                                                >  2:2737345/1‑148 (MQ=255)
                                                                                                 ggAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATc                                               >  2:2083970/1‑148 (MQ=255)
                                                                                                 ggAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATc                                               >  2:2321780/1‑148 (MQ=255)
                                                                                                 ggAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATc                                               >  2:289590/1‑148 (MQ=255)
                                                                                                  gAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcg                                              >  1:219203/1‑148 (MQ=255)
                                                                                                    aCTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcgcg                                            >  2:2705588/1‑148 (MQ=255)
                                                                                                     cTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCgg                                           >  2:2269900/1‑148 (MQ=255)
                                                                                                     cTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCgg                                           >  2:2658003/1‑148 (MQ=255)
                                                                                                          ggCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGt                                      <  1:2529415/148‑1 (MQ=255)
                                                                                                                gCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCa                                >  1:3029311/1‑148 (MQ=255)
                                                                                                                 ccTTCTTGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCAt                               >  2:3029085/1‑148 (MQ=255)
                                                                                                                  cttcttGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATc                              <  2:3138019/148‑1 (MQ=255)
                                                                                                                  cttcttGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATc                              <  2:1169873/148‑1 (MQ=255)
                                                                                                                    tcttGCTCCATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCa                            <  2:503308/148‑1 (MQ=255)
                                                                                                                    tcttGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCa                            >  2:1750569/1‑148 (MQ=255)
                                                                                                                          tCCATTCGTGCGTCCAGTACCAGACAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGCTAAGTGTCGCAGTCATCGCAGCTTTAGTCTTTAATCGCGGCCTGTTCGTCATCCAGCGAcc                      >  1:2016872/1‑148 (MQ=255)
                                                                                                                           ccATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCa                     <  2:1274891/148‑1 (MQ=255)
                                                                                                                           ccATTCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCa                     <  2:3049889/148‑1 (MQ=255)
                                                                                                                             aTTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGt                   >  1:469454/1‑148 (MQ=255)
                                                                                                                              ttCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  >  1:950252/1‑148 (MQ=255)
                                                                                                                              ttCGTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTa                  >  1:2726210/1‑148 (MQ=255)
                                                                                                                               tCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTaa                 <  1:72059/148‑1 (MQ=255)
                                                                                                                               tCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTaa                 <  2:1792990/148‑1 (MQ=255)
                                                                                                                               tCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTaa                 >  1:2025394/1‑148 (MQ=255)
                                                                                                                                 gTGCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATg               <  2:2979062/148‑1 (MQ=255)
                                                                                                                                  tGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGc              >  2:3112886/1‑148 (MQ=255)
                                                                                                                                  tGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGc              <  1:2767541/148‑1 (MQ=255)
                                                                                                                                   gCGTCCAGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCa             <  2:29443/148‑1 (MQ=255)
                                                                                                                                       ccAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGttt         >  1:2210931/1‑148 (MQ=255)
                                                                                                                                        cAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGCGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTg        <  2:2039106/148‑1 (MQ=255)
                                                                                                                                        cAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTg        <  2:2746536/148‑1 (MQ=255)
                                                                                                                                         aGTACCTGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGc       >  2:3395986/1‑148 (MQ=255)
                                                                                                                                         aGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGc       <  2:3185492/148‑1 (MQ=255)
                                                                                                                                           tACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCg     >  2:494817/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGTCGTCACACCCTATGAGCTGTAATCCGGGCGACCTGCGCCACCCGGAAAACGTCAGAGATTACTCTCCTTCACCCGGAAACAGGAACGGGTTGATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTC  >  CP009273/4073598‑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: ‑