Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,932,701 C→A 53.1% Q266K (CAG→AAG)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,932,7010CA53.1% 29.3 / 273.3 162Q266K (CAG→AAG) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base A (38/48);  total (73/89)
Fisher's exact test for biased strand distribution p-value = 8.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 166 total aligned reads displayed.

GATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTG  >  CP009273/1932554‑1932848
                                                                                                                                                   |                                                                                                                                                   
gATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTc                                                                                                                                                     <  1:1744336/148‑1 (MQ=255)
     gAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGaaa                                                                                                                                                <  2:2685613/148‑1 (MQ=255)
      aaaaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAg                                                                                                                                               >  2:1975414/1‑148 (MQ=255)
      aaaaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAg                                                                                                                                               <  2:2301107/148‑1 (MQ=255)
      aaaaTTGTTGCCAAGGGTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAg                                                                                                                                               <  2:3454589/148‑1 (MQ=255)
        aaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCg                                                                                                                                             >  1:2053872/1‑148 (MQ=255)
        aaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCg                                                                                                                                             <  2:303884/148‑1 (MQ=255)
        aaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCg                                                                                                                                             >  1:2053862/1‑148 (MQ=255)
        aaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGAGACCTCTGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCg                                                                                                                                             <  1:3399111/148‑1 (MQ=255)
         aTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGt                                                                                                                                            >  1:3233874/1‑148 (MQ=255)
          ttgttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGtt                                                                                                                                           >  1:2446255/1‑148 (MQ=255)
          ttgttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGtt                                                                                                                                           >  1:2938670/1‑148 (MQ=255)
           tgttgCCAAGGTTGAACGTGCGGACGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTg                                                                                                                                          >  2:1249545/1‑148 (MQ=255)
            gttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGa                                                                                                                                         <  1:3133087/148‑1 (MQ=255)
            gttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGa                                                                                                                                         <  2:1126012/148‑1 (MQ=255)
             ttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGAt                                                                                                                                        >  1:1124332/1‑148 (MQ=255)
             ttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGAt                                                                                                                                        <  1:1534070/148‑1 (MQ=255)
             ttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGAt                                                                                                                                        <  2:2388753/148‑1 (MQ=255)
                ccAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCcg                                                                                                                                     >  2:665647/1‑148 (MQ=255)
                 cAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgt                                                                                                                                    <  1:2771798/148‑1 (MQ=255)
                 cAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCcgt                                                                                                                                    >  2:1756099/1‑148 (MQ=255)
                    ggTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcgt                                                                                                                                 >  2:756025/1‑148 (MQ=255)
                     gTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTACTGGTTGCACGTGGCGACCTCGGTGTGGACATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTg                                                                                                                                >  2:3274482/1‑148 (MQ=255)
                            gTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCa                                                                                                                         <  1:477777/148‑1 (MQ=255)
                             tGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAg                                                                                                                        <  2:1809918/148‑1 (MQ=255)
                              gCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGc                                                                                                                       >  2:869401/1‑148 (MQ=255)
                                  aaGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTaaa                                                                                                                   >  1:3022332/1‑148 (MQ=255)
                                  aaGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTaaa                                                                                                                   <  1:1756099/148‑1 (MQ=255)
                                     ccGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCg                                                                                                                >  1:733158/1‑148 (MQ=255)
                                      cGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGa                                                                                                               >  1:2878432/1‑148 (MQ=255)
                                      cGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGGGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGa                                                                                                               >  2:1735256/1‑148 (MQ=255)
                                         ttGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCg                                                                                                            <  2:1968056/148‑1 (MQ=255)
                                           gCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGt                                                                                                          <  1:3535921/148‑1 (MQ=255)
                                             aGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTaa                                                                                                        <  1:3518243/148‑1 (MQ=255)
                                               ccAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATc                                                                                                      >  2:451090/1‑148 (MQ=255)
                                               ccAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATc                                                                                                      <  1:209707/148‑1 (MQ=255)
                                                cAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCa                                                                                                     >  2:2634818/1‑148 (MQ=255)
                                                  ggATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACg                                                                                                   <  2:1748186/148‑1 (MQ=255)
                                                  ggATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACg                                                                                                   <  2:3336926/148‑1 (MQ=255)
                                                  ggATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACg                                                                                                   <  2:588036/148‑1 (MQ=255)
                                                     tGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCg                                                                                                >  2:1531075/1‑148 (MQ=255)
                                                         aTGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGAccc                                                                                            >  2:2946989/1‑148 (MQ=255)
                                                         aTGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATACGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGAccc                                                                                            <  2:1988423/148‑1 (MQ=255)
                                                          tGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCa                                                                                           <  1:2634818/148‑1 (MQ=255)
                                                           ggATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAg                                                                                          <  1:577433/148‑1 (MQ=255)
                                                           ggATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAg                                                                                          <  1:2860746/148‑1 (MQ=255)
                                                                aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAgatgat                                                                                     >  1:2710172/1‑148 (MQ=255)
                                                                 catcatCCTCGCCTCTGAGGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGATCTGTTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGttgatg                                                                                    <  2:3364636/148‑1 (MQ=255)
                                                                 catcatCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGTTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGatgatg                                                                                    <  2:3364490/148‑1 (MQ=255)
                                                                   tcatcCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGa                                                                                  >  1:926641/1‑148 (MQ=255)
                                                                     atcCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGt                                                                                <  1:1975922/148‑1 (MQ=255)
                                                                       ccTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCa                                                                              <  2:3464673/148‑1 (MQ=255)
                                                                       ccTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCa                                                                              <  2:2256669/148‑1 (MQ=255)
                                                                       ccTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCa                                                                              <  1:3185918/148‑1 (MQ=255)
                                                                       ccTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCa                                                                              >  1:1185367/1‑148 (MQ=255)
                                                                        cTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCaa                                                                             <  2:80986/148‑1 (MQ=255)
                                                                         tCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAAt                                                                            <  2:2495663/148‑1 (MQ=255)
                                                                         tCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAAt                                                                            >  1:2375702/1‑148 (MQ=255)
                                                                          cGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATg                                                                           >  2:323045/1‑148 (MQ=255)
                                                                           gCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGa                                                                          <  2:2907715/148‑1 (MQ=255)
                                                                             ctctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGAtt                                                                        <  2:3312969/148‑1 (MQ=255)
                                                                             ctctGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGAtt                                                                        >  1:2937265/1‑148 (MQ=255)
                                                                                 gACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTa                                                                    >  2:253384/1‑148 (MQ=255)
                                                                                 gACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTa                                                                    >  1:3054907/1‑148 (MQ=255)
                                                                                  aCGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTaa                                                                   >  2:2894191/1‑148 (MQ=255)
                                                                                   cGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAAc                                                                  >  1:2905602/1‑148 (MQ=255)
                                                                                    gTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAAcc                                                                 <  1:3316064/148‑1 (MQ=255)
                                                                                      ggTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCg                                                               >  1:1735808/1‑148 (MQ=255)
                                                                                      ggTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCg                                                               >  1:1735728/1‑148 (MQ=255)
                                                                                           tGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGcc                                                          >  1:1205747/1‑148 (MQ=255)
                                                                                                 cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcgcg                                                    >  2:2130380/1‑148 (MQ=255)
                                                                                                 cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcgcg                                                    >  1:168399/1‑148 (MQ=255)
                                                                                                 cACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcgcg                                                    >  1:168367/1‑148 (MQ=255)
                                                                                                   cGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTg                                                  >  2:2988988/1‑148 (MQ=255)
                                                                                                     tGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCa                                                >  1:3174924/1‑148 (MQ=255)
                                                                                                          aCCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGt                                           <  1:503697/148‑1 (MQ=255)
                                                                                                          aCCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGt                                           <  2:2740262/148‑1 (MQ=255)
                                                                                                          aCCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGt                                           <  1:329323/148‑1 (MQ=255)
                                                                                                           ccTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTc                                          <  1:1437158/148‑1 (MQ=255)
                                                                                                             tCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCAt                                        <  2:2053872/148‑1 (MQ=255)
                                                                                                             tCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCAt                                        <  2:2053862/148‑1 (MQ=255)
                                                                                                               ggTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATgg                                      <  1:3116517/148‑1 (MQ=255)
                                                                                                               ggTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATgg                                      <  1:16329/148‑1 (MQ=255)
                                                                                                                  gtgGATATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACg                                   <  2:2431339/148‑1 (MQ=255)
                                                                                                                   tgGAAATTGGCGACCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGt                                  <  1:1989514/148‑1 (MQ=255)
                                                                                                                       aaTTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCa                              <  1:1975414/148‑1 (MQ=255)
                                                                                                                               aCCCGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGt                      >  1:1326381/1‑148 (MQ=255)
                                                                                                                                cccGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGtt                     >  2:1576596/1‑148 (MQ=255)
                                                                                                                                cccGGAACTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGtt                     >  1:1404032/1‑148 (MQ=255)
                                                                                                                                     aaCTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGa                <  1:884512/148‑1 (MQ=255)
                                                                                                                                       cTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATg              <  2:1287294/148‑1 (MQ=255)
                                                                                                                                       cTGGTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATg              >  1:744278/1‑148 (MQ=255)
                                                                                                                                         ggTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGATTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGt            >  1:3500810/1‑148 (MQ=255)
                                                                                                                                         ggTCGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGt            <  2:542928/148‑1 (MQ=255)
                                                                                                                                          gTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTa           >  1:1419783/1‑148 (MQ=255)
                                                                                                                                            cGGCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACt         <  1:1523915/148‑1 (MQ=255)
                                                                                                                                              gCATTAAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGa       >  1:2120667/1‑148 (MQ=255)
                                                                                                                                                aTTAAGCAAGCGATGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTTATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACg     <  2:3445275/148‑1 (MQ=255)
                                                                                                                                                 ttCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGc    >  2:1142536/1‑148 (MQ=255)
                                                                                                                                                   cAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCtg  >  1:1041246/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTG  >  CP009273/1932554‑1932848

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