Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 646,966 A→T 7.3% 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 7.3% 689.1 / 37.2 232V116E (GTA→GAA) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base A (126/89);  new base T (12/5);  total (138/94)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 236 total aligned reads displayed.

ATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGCGGACGG  >  CP009273/646819‑647113
                                                                                                                                                   |                                                                                                                                                   
    gATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCt                                                                                                                                                 >  2:3350446/1‑148 (MQ=255)
    gATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCt                                                                                                                                                 >  2:1588970/1‑148 (MQ=255)
      tAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTgg                                                                                                                                               >  2:3030475/1‑148 (MQ=255)
        aCGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTc                                                                                                                                             >  1:2902362/1‑148 (MQ=255)
        aCGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTc                                                                                                                                             >  1:1854526/1‑148 (MQ=255)
         cGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACAATCACGCCGGGTACGCTGGTCa                                                                                                                                            >  2:2242184/1‑148 (MQ=255)
            tGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCg                                                                                                                                         >  1:2794006/1‑148 (MQ=255)
            tGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCg                                                                                                                                         >  1:2792666/1‑148 (MQ=255)
              aCCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTg                                                                                                                                       >  2:3186644/1‑148 (MQ=255)
                 aTTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGaa                                                                                                                                    >  2:481051/1‑148 (MQ=255)
                     gTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAAc                                                                                                                                <  1:702206/148‑1 (MQ=255)
                      tAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAAcc                                                                                                                               >  1:1562760/1‑148 (MQ=255)
                      tAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAAcc                                                                                                                               >  1:2656700/1‑148 (MQ=255)
                        aaGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGc                                                                                                                             <  2:283907/148‑1 (MQ=255)
                           ggATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACt                                                                                                                          >  1:2934262/1‑148 (MQ=255)
                            gATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTg                                                                                                                         <  1:2838422/148‑1 (MQ=255)
                                 gCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGa                                                                                                                    >  2:3261425/1‑148 (MQ=255)
                                  cGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGaa                                                                                                                   <  2:2484573/148‑1 (MQ=255)
                                    ttCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACa                                                                                                                 >  2:595620/1‑148 (MQ=255)
                                     tcatcaCAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAg                                                                                                                >  2:1632594/1‑148 (MQ=255)
                                     tcatcaCAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAg                                                                                                                <  2:2562805/148‑1 (MQ=255)
                                     tcatcaCAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAg                                                                                                                <  2:2959117/148‑1 (MQ=255)
                                     tcatcaCAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAg                                                                                                                >  2:782664/1‑148 (MQ=255)
                                     tcatcaCAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAg                                                                                                                >  2:1765326/1‑148 (MQ=255)
                                         cacaCTGCAGCAAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGcc                                                                                                            <  2:621042/148‑1 (MQ=255)
                                         cacaATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGcc                                                                                                            <  2:1076468/148‑1 (MQ=255)
                                         cacaATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGcc                                                                                                            >  2:3212235/1‑148 (MQ=255)
                                              tGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTa                                                                                                       >  1:667963/1‑148 (MQ=255)
                                                cAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTAtt                                                                                                     >  2:557745/1‑148 (MQ=255)
                                                cAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTAtt                                                                                                     <  1:564546/148‑1 (MQ=255)
                                                      aTCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGttt                                                                                               >  2:1309740/1‑148 (MQ=255)
                                                       tCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGtttt                                                                                              >  2:23429/1‑148 (MQ=255)
                                                       tCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGtttt                                                                                              >  1:3538820/1‑148 (MQ=255)
                                                        cTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTg                                                                                             >  1:1055456/1‑148 (MQ=255)
                                                        cTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTg                                                                                             >  1:1904784/1‑148 (MQ=255)
                                                         ttGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTgg                                                                                            <  2:4508/148‑1 (MQ=255)
                                                          tGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGt                                                                                           >  2:1489771/1‑148 (MQ=255)
                                                          tGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGt                                                                                           >  1:2341929/1‑148 (MQ=255)
                                                          tGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGt                                                                                           >  1:3313110/1‑148 (MQ=255)
                                                            ttCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtat                                                                                         >  2:757811/1‑148 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTCTTCGGTTTTGGtata                                                                                        >  2:2404060/1‑148 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtata                                                                                        >  1:2627445/1‑148 (MQ=255)
                                                             tCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGtata                                                                                        >  2:832767/1‑148 (MQ=255)
                                                              cccTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataa                                                                                       <  2:731000/148‑1 (MQ=255)
                                                              cccTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataa                                                                                       >  1:837273/1‑148 (MQ=255)
                                                                cTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaat                                                                                     >  2:2814223/1‑148 (MQ=255)
                                                                cTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaat                                                                                     <  2:2791161/148‑1 (MQ=255)
                                                                cTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTataaat                                                                                     <  2:1726923/148‑1 (MQ=255)
                                                                   gATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCACGTATTCGGTTTTGGTAtaaataaa                                                                                  >  2:1221322/1‑148 (MQ=255)
                                                                      gACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAg                                                                               <  1:11411/148‑1 (MQ=255)
                                                                       aCGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGa                                                                              <  2:1040163/148‑1 (MQ=255)
                                                                       aCGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGa                                                                              <  2:1071081/148‑1 (MQ=255)
                                                                        cGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGAt                                                                             <  1:1926256/148‑1 (MQ=255)
                                                                         gAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATg                                                                            >  1:3437997/1‑148 (MQ=255)
                                                                         gAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATg                                                                            <  1:2436022/148‑1 (MQ=255)
                                                                            aTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGt                                                                         >  2:218233/1‑148 (MQ=255)
                                                                              tttttCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGc                                                                       >  2:452872/1‑148 (MQ=255)
                                                                              tttttCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGc                                                                       >  2:835372/1‑148 (MQ=255)
                                                                                        aTTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTc                                                             >  1:1191695/1‑148 (MQ=255)
                                                                                         ttCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCa                                                            >  2:2541708/1‑148 (MQ=255)
                                                                                          tCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGGCGCTCAt                                                           <  2:3236930/148‑1 (MQ=255)
                                                                                           cGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATa                                                          <  2:1176912/148‑1 (MQ=255)
                                                                                           cGGCATAGCGTTTCAGTCGCGTGGCGCTGTTGTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATa                                                          <  2:1176858/148‑1 (MQ=255)
                                                                                              cATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGc                                                       >  1:2171213/1‑148 (MQ=255)
                                                                                                tAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGa                                                     >  1:3278036/1‑148 (MQ=255)
                                                                                                  gCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAgg                                                   >  2:1222360/1‑148 (MQ=255)
                                                                                                  gCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAgg                                                   >  2:1222491/1‑148 (MQ=255)
                                                                                                   cGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGt                                                  <  1:3197562/148‑1 (MQ=255)
                                                                                                    gTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGtt                                                 <  1:1418768/148‑1 (MQ=255)
                                                                                                         aGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCa                                            >  2:2068398/1‑148 (MQ=255)
                                                                                                         aGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCa                                            >  2:273135/1‑148 (MQ=255)
                                                                                                          gTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCaa                                           <  2:2199377/148‑1 (MQ=255)
                                                                                                          gTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCaa                                           >  1:2422880/1‑148 (MQ=255)
                                                                                                           tCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAAt                                          >  2:1381167/1‑148 (MQ=255)
                                                                                                           tCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAAt                                          <  2:3083783/148‑1 (MQ=255)
                                                                                                            cgcgTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAAtt                                         >  2:1562724/1‑148 (MQ=255)
                                                                                                             gcgTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTc                                        >  1:2731772/1‑148 (MQ=255)
                                                                                                                 ggCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTg                                    <  2:3330971/148‑1 (MQ=255)
                                                                                                                 ggCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTg                                    <  2:1675484/148‑1 (MQ=255)
                                                                                                                  gcgcTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTgg                                   <  2:2866619/148‑1 (MQ=255)
                                                                                                                     cTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTa                                >  2:3515782/1‑148 (MQ=255)
                                                                                                                     cTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTa                                <  1:3186644/148‑1 (MQ=255)
                                                                                                                     cTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTa                                >  2:311974/1‑148 (MQ=255)
                                                                                                                     cTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTa                                >  2:1665219/1‑148 (MQ=255)
                                                                                                                         tttCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGt                            >  1:315870/1‑148 (MQ=255)
                                                                                                                           tCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCa                          <  1:2996384/148‑1 (MQ=255)
                                                                                                                               tCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAgcgc                      >  1:1615363/1‑148 (MQ=255)
                                                                                                                                cAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGcc                     >  1:1016478/1‑148 (MQ=255)
                                                                                                                                     ccATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTCATGTCAGCGCCAGTcc                >  2:1372062/1‑148 (MQ=255)
                                                                                                                                      cATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCt               >  2:2011851/1‑148 (MQ=255)
                                                                                                                                       aTCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCtt              <  1:1998738/148‑1 (MQ=255)
                                                                                                                                       aTCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCtt              >  2:1758403/1‑148 (MQ=255)
                                                                                                                                        tCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTc             >  2:1972265/1‑148 (MQ=255)
                                                                                                                                         cACGCCGGGTTCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCa            <  2:2656700/148‑1 (MQ=255)
                                                                                                                                           cGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCAcc          >  2:1938695/1‑148 (MQ=255)
                                                                                                                                              cGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACcgcg       <  1:1289456/148‑1 (MQ=255)
                                                                                                                                               gggTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcgg      <  1:2542334/148‑1 (MQ=255)
                                                                                                                                                ggTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcgga     >  2:2613089/1‑148 (MQ=255)
                                                                                                                                                  tACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcggacg   <  1:2115551/148‑1 (MQ=255)
                                                                                                                                                   aCGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGcggacgg  >  1:625864/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACAGCGCCTCGTATTCGGTTTTGGTATAAATAAACAGATGCGTGCTGTGCCGCTCATAGGCGAGGTTTATCAATTCAGTGGCTAATGTCAGCGCCAGTCCTTCACCGCGGACGG  >  CP009273/646819‑647113

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