Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 646,966 A→T 6.5% V116E (GTA→GAA)  citC ← [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273646,9660AT 6.5% 563.2 / 23.8 185V116E (GTA→GAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base A (95/78);  new base T (6/6);  total (101/84)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 186 total aligned reads displayed.

ATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGCGGACG  >  CP009273/646819‑647112
                                                                                                                                                   |                                                                                                                                                  
aTCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCTAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTa                                                                                                                                                    >  1:1893525/1‑148 (MQ=255)
 tCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTAc                                                                                                                                                   >  1:2414753/1‑148 (MQ=255)
  cAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACg                                                                                                                                                  >  2:2840633/1‑148 (MQ=255)
  cAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACg                                                                                                                                                  >  2:2417017/1‑148 (MQ=255)
    gATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCt                                                                                                                                                >  1:1167619/1‑148 (MQ=255)
      tAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTgg                                                                                                                                              >  2:2184029/1‑148 (MQ=255)
             gACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGt                                                                                                                                       <  2:1066253/148‑1 (MQ=255)
              aCCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTg                                                                                                                                      >  1:1214964/1‑148 (MQ=255)
                   tCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGaaaa                                                                                                                                 >  1:1917594/1‑148 (MQ=255)
                      tAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGACAAAcc                                                                                                                              >  2:1446250/1‑148 (MQ=255)
                       aaaGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCg                                                                                                                             >  2:2935057/1‑148 (MQ=255)
                       aaaGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCg                                                                                                                             >  2:1458094/1‑148 (MQ=255)
                        aaGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGc                                                                                                                            <  1:2075308/148‑1 (MQ=255)
                          gggATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCAc                                                                                                                          <  2:937215/148‑1 (MQ=255)
                           ggATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACt                                                                                                                         <  2:3061877/148‑1 (MQ=255)
                            gATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTg                                                                                                                        <  1:1248657/148‑1 (MQ=255)
                              ttGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGc                                                                                                                      >  1:1287981/1‑148 (MQ=255)
                               tGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCg                                                                                                                     >  2:859878/1‑148 (MQ=255)
                                ggCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCgg                                                                                                                    >  1:2931095/1‑148 (MQ=255)
                                 gCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGa                                                                                                                   >  1:658089/1‑148 (MQ=255)
                                    ttCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACa                                                                                                                >  1:288967/1‑148 (MQ=255)
                                            aaTGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCg                                                                                                        >  2:2679843/1‑148 (MQ=255)
                                             aTGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGt                                                                                                       <  1:2532778/148‑1 (MQ=255)
                                               gCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTAt                                                                                                     <  2:2305488/148‑1 (MQ=255)
                                                   ccAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCgg                                                                                                 >  2:1197111/1‑148 (MQ=255)
                                                    cAAACTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACAGGCGGAACAGCGCCTCGTATTCGGt                                                                                                >  2:160484/1‑148 (MQ=255)
                                                      aTCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGttt                                                                                              >  2:2003700/1‑148 (MQ=255)
                                                      aTCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGttt                                                                                              >  2:665618/1‑148 (MQ=255)
                                                       tCTTGGTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGtttt                                                                                             <  2:2556477/148‑1 (MQ=255)
                                                        cTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTg                                                                                            <  1:1255871/148‑1 (MQ=255)
                                                        cTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTg                                                                                            >  2:457924/1‑148 (MQ=255)
                                                          tGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGt                                                                                          <  2:1021632/148‑1 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtata                                                                                       >  2:583441/1‑148 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtata                                                                                       >  1:1674932/1‑148 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtata                                                                                       >  1:2372962/1‑148 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtata                                                                                       >  2:3158131/1‑148 (MQ=255)
                                                                 tGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaata                                                                                   <  1:620444/148‑1 (MQ=255)
                                                                 tGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaata                                                                                   <  1:620165/148‑1 (MQ=255)
                                                                 tGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaata                                                                                   <  1:2824623/148‑1 (MQ=255)
                                                                 tGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaata                                                                                   <  1:2824650/148‑1 (MQ=255)
                                                                  ggATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaataa                                                                                  <  1:3208998/148‑1 (MQ=255)
                                                                   gATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTAtaaataaa                                                                                 <  2:1898893/148‑1 (MQ=255)
                                                                   gATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTAtaaataaa                                                                                 >  1:850826/1‑148 (MQ=255)
                                                                   gATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTAtaaataaa                                                                                 >  1:850128/1‑148 (MQ=255)
                                                                   gATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTAtaaataaa                                                                                 <  1:521502/148‑1 (MQ=255)
                                                                      gACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAg                                                                              >  2:3199103/1‑148 (MQ=255)
                                                                      gACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAg                                                                              >  1:1674374/1‑148 (MQ=255)
                                                                      gACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAg                                                                              <  2:2451278/148‑1 (MQ=255)
                                                                        cGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACCGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGAt                                                                            >  2:1557151/1‑148 (MQ=255)
                                                                         gAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATg                                                                           >  1:695206/1‑148 (MQ=255)
                                                                         gAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATg                                                                           <  2:1213844/148‑1 (MQ=255)
                                                                          aaaTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGc                                                                          <  1:1480997/148‑1 (MQ=255)
                                                                           aaTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCg                                                                         <  2:1315226/148‑1 (MQ=255)
                                                                            aTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGt                                                                        >  1:3121369/1‑148 (MQ=255)
                                                                             ttttttCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTg                                                                       <  2:2412790/148‑1 (MQ=255)
                                                                              tttttCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGc                                                                      <  2:1186041/148‑1 (MQ=255)
                                                                              tttttCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGc                                                                      <  2:657015/148‑1 (MQ=255)
                                                                                tttCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCtg                                                                    <  2:701751/148‑1 (MQ=255)
                                                                                      cGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGc                                                              >  1:1371546/1‑148 (MQ=255)
                                                                                      cGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGc                                                              >  2:1789231/1‑148 (MQ=255)
                                                                                       gATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCt                                                             <  2:1907039/148‑1 (MQ=255)
                                                                                       gATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCt                                                             <  1:1531916/148‑1 (MQ=255)
                                                                                        aTTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTc                                                            <  2:2665656/148‑1 (MQ=255)
                                                                                          tCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCAt                                                          >  1:262727/1‑148 (MQ=255)
                                                                                          tCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCAt                                                          <  1:2651658/148‑1 (MQ=255)
                                                                                              cATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGc                                                      <  1:2184029/148‑1 (MQ=255)
                                                                                                tAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGa                                                    >  2:2292455/1‑148 (MQ=255)
                                                                                                 aGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAg                                                   >  1:234972/1‑148 (MQ=255)
                                                                                                  gCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAgg                                                  <  2:1640038/148‑1 (MQ=255)
                                                                                                     tttCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGttt                                               <  2:2535728/148‑1 (MQ=255)
                                                                                                      ttCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTTGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTa                                              <  1:1333568/148‑1 (MQ=255)
                                                                                                      ttCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTa                                              <  2:1888731/148‑1 (MQ=255)
                                                                                                      ttCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTa                                              >  1:1797128/1‑148 (MQ=255)
                                                                                                      ttCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTa                                              >  1:2723674/1‑148 (MQ=255)
                                                                                                      ttCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTa                                              >  1:382959/1‑148 (MQ=255)
                                                                                                       tCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTAt                                             <  2:604189/148‑1 (MQ=255)
                                                                                                        cAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATc                                            >  1:1394641/1‑148 (MQ=255)
                                                                                                        cAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATc                                            >  1:1499519/1‑148 (MQ=255)
                                                                                                               gTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAg                                     <  1:3218614/148‑1 (MQ=255)
                                                                                                               gTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAg                                     <  1:1603850/148‑1 (MQ=255)
                                                                                                                      tGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTaa                              <  1:449295/148‑1 (MQ=255)
                                                                                                                      tGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTaa                              >  1:3039227/1‑148 (MQ=255)
                                                                                                                       gTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAAt                             >  1:660188/1‑148 (MQ=255)
                                                                                                                        ttttCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATg                            <  1:3199103/148‑1 (MQ=255)
                                                                                                                         tttCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGt                           <  2:1167619/148‑1 (MQ=255)
                                                                                                                          ttCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTc                          >  2:620998/1‑148 (MQ=255)
                                                                                                                           tCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCa                         <  2:1893525/148‑1 (MQ=255)
                                                                                                                           tCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCa                         <  2:2585865/148‑1 (MQ=255)
                                                                                                                            ccATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAg                        <  2:1845792/148‑1 (MQ=255)
                                                                                                                            ccATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAg                        <  2:50962/148‑1 (MQ=255)
                                                                                                                                cAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGcc                    >  2:2033255/1‑148 (MQ=255)
                                                                                                                                      cATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCt              >  1:1646664/1‑148 (MQ=255)
                                                                                                                                      cATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCt              >  1:2943373/1‑148 (MQ=255)
                                                                                                                                      cATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCt              >  2:1292658/1‑148 (MQ=255)
                                                                                                                                        tCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTc            >  1:1075054/1‑148 (MQ=255)
                                                                                                                                        tCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTc            <  2:1582760/148‑1 (MQ=255)
                                                                                                                                           cGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCAcc         >  2:3264051/1‑148 (MQ=255)
                                                                                                                                                ggTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcgga    >  2:356969/1‑148 (MQ=255)
                                                                                                                                                 gTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCTTCTTCGGTTTTGGTAGACATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcggac   >  1:3249055/1‑148 (MQ=255)
                                                                                                                                                  tACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcggacg  >  2:2967074/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGCGGACG  >  CP009273/646819‑647112

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: ‑