Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,013 C→A 5.1% D365Y (GAT→TAT)  gudX ← glucarate dehydratase‑related protein, substrate unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,913,0130CA 5.1% 773.2 / 11.7 201D365Y (GAT→TAT) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (96/94);  new base A (3/8);  total (99/102)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

ATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCA  >  CP009273/2912867‑2913160
                                                                                                                                                  |                                                                                                                                                   
ataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCg                                                                                                                                                    <  1:2330662/148‑1 (MQ=255)
ataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCg                                                                                                                                                    <  1:2015305/148‑1 (MQ=255)
  aGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGAt                                                                                                                                                  <  2:2786807/148‑1 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAg                                                                                                                                                >  2:3156362/1‑148 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAg                                                                                                                                                >  1:400185/1‑148 (MQ=255)
     ccTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAga                                                                                                                                               <  2:1802157/148‑2 (MQ=255)
      cTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAgag                                                                                                                                              <  1:1425591/148‑3 (MQ=255)
        cATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGt                                                                                                                                            <  2:300088/148‑1 (MQ=255)
        cATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGt                                                                                                                                            <  2:1678760/148‑1 (MQ=255)
             gCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGAt                                                                                                                                       <  1:3149507/148‑1 (MQ=255)
             gCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGAt                                                                                                                                       >  1:1712943/1‑148 (MQ=255)
             gCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGAt                                                                                                                                       >  1:1297925/1‑148 (MQ=255)
               cTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTa                                                                                                                                     <  2:372646/148‑1 (MQ=255)
                 tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTAcc                                                                                                                                   >  1:1883047/1‑148 (MQ=255)
                 tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTAcc                                                                                                                                   <  1:376331/148‑1 (MQ=255)
                    gTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCgg                                                                                                                                <  2:1916390/148‑1 (MQ=255)
                     tACCTGTTCCCAGTCCAGGTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGt                                                                                                                               <  1:749108/148‑1 (MQ=255)
                       ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgc                                                                                                                             >  1:1136621/1‑148 (MQ=255)
                        cTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgcc                                                                                                                            >  2:3250100/1‑148 (MQ=255)
                          gTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgccgc                                                                                                                          <  1:790597/148‑1 (MQ=255)
                           ttCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgccgcc                                                                                                                         >  2:1037244/1‑148 (MQ=255)
                             cccAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCcgccgc                                                                                                                       >  2:407772/1‑148 (MQ=255)
                                agtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGcc                                                                                                                    <  2:1503324/148‑1 (MQ=255)
                                  tccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATAATTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGACCTCCAGCCAAATCAAGTGGGAATCGAAAGATGAAGGATTAAACGGTGACGCCGCGCCca                                                                                                                  >  1:1065980/1‑148 (MQ=255)
                                    cagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCcaca                                                                                                                <  2:519605/148‑1 (MQ=255)
                                     agtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACAt                                                                                                               >  2:2607105/1‑148 (MQ=255)
                                     agtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACAt                                                                                                               >  2:2383016/1‑148 (MQ=255)
                                      gtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATg                                                                                                              <  2:83/148‑1 (MQ=255)
                                      gtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATg                                                                                                              <  1:3517665/148‑1 (MQ=255)
                                      gtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATg                                                                                                              <  2:2179236/148‑1 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGTAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             >  1:2424202/1‑148 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             <  2:3597983/148‑1 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             <  2:2315649/148‑1 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             <  2:3597688/148‑1 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             >  2:71867/1‑148 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             >  1:1087487/1‑148 (MQ=255)
                                            cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                        <  1:1696380/148‑1 (MQ=255)
                                            cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                        <  1:602913/148‑1 (MQ=255)
                                              cccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACa                                                                                                      <  1:2576814/148‑1 (MQ=255)
                                               ccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACAt                                                                                                     <  1:1778398/148‑1 (MQ=255)
                                                   aCCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGcc                                                                                                 <  2:2467073/148‑1 (MQ=255)
                                                     cGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCag                                                                                               <  2:3508302/148‑1 (MQ=255)
                                                          gcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAgagaga                                                                                          <  2:708242/148‑1 (MQ=255)
                                                           cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                         >  1:410074/1‑148 (MQ=255)
                                                           cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                         >  1:2332772/1‑148 (MQ=255)
                                                           cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                         <  1:2009911/148‑1 (MQ=255)
                                                           cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                         <  2:2705258/148‑1 (MQ=255)
                                                             tCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGatat                                                                                       <  1:2940659/148‑1 (MQ=255)
                                                             tCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGatat                                                                                       >  2:1707656/1‑148 (MQ=255)
                                                              cAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGTTAACAAACGCCCGAGAGATATc                                                                                      >  2:1935206/1‑148 (MQ=255)
                                                              cAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATc                                                                                      >  2:2931047/1‑148 (MQ=255)
                                                              cAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATc                                                                                      >  2:1166491/1‑148 (MQ=255)
                                                                 gAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGaa                                                                                   <  2:3633540/148‑1 (MQ=255)
                                                                  aaCGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGaaa                                                                                  >  1:266516/1‑148 (MQ=255)
                                                                    cGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATg                                                                                <  2:2849606/148‑1 (MQ=255)
                                                                    cGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATg                                                                                <  1:2678258/148‑1 (MQ=255)
                                                                      gCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGt                                                                              >  2:270335/1‑148 (MQ=255)
                                                                       cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                             <  1:3105000/148‑1 (MQ=255)
                                                                       cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                             <  1:2996759/148‑1 (MQ=255)
                                                                         aTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttat                                                                           >  1:991190/1‑148 (MQ=255)
                                                                         aTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttat                                                                           >  1:2380276/1‑148 (MQ=255)
                                                                          tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCCGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAACTGGttctt                                                                          >  2:2433204/1‑148 (MQ=37)
                                                                          tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                          >  1:1728312/1‑148 (MQ=255)
                                                                          tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                          <  1:19979/148‑1 (MQ=255)
                                                                          tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                          >  1:450721/1‑148 (MQ=255)
                                                                               ccGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAAt                                                                     <  2:2247501/148‑1 (MQ=255)
                                                                                   ttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCa                                                                 >  2:1410406/1‑148 (MQ=255)
                                                                                    tttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAg                                                                <  1:1032881/148‑1 (MQ=255)
                                                                                     ttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGc                                                               <  1:909710/148‑1 (MQ=255)
                                                                                       aaTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGccc                                                             <  1:1680294/148‑1 (MQ=255)
                                                                                       aaTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGccc                                                             <  1:2788036/148‑1 (MQ=255)
                                                                                       aaTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGccc                                                             <  1:3447973/148‑1 (MQ=255)
                                                                                        aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                            >  2:1974053/1‑148 (MQ=255)
                                                                                        aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                            >  2:1973984/1‑148 (MQ=255)
                                                                                        aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAAGGGTTCTTAGAATGGCAGcccc                                                            >  2:2154395/1‑148 (MQ=255)
                                                                                         tctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCa                                                           <  2:137071/148‑1 (MQ=255)
                                                                                          ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCTCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                          >  1:3813524/1‑148 (MQ=255)
                                                                                          ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                          >  1:1003518/1‑148 (MQ=255)
                                                                                          ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                          >  2:1858157/1‑148 (MQ=255)
                                                                                          ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                          >  2:2910970/1‑148 (MQ=255)
                                                                                          ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGGTATTAGAATGGCAGCCCCAg                                                          >  2:2956108/1‑148 (MQ=255)
                                                                                           tcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAgg                                                         <  2:3175250/148‑1 (MQ=255)
                                                                                            ccAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGt                                                        >  2:3314683/1‑148 (MQ=255)
                                                                                               gCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAg                                                     >  1:678437/1‑148 (MQ=255)
                                                                                                  gATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGTTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCAAATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGccc                                                  <  2:1282840/148‑1 (MQ=255)
                                                                                                   aTTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcccc                                                 >  1:3626242/1‑148 (MQ=255)
                                                                                                          gTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgt                                          >  1:3294343/1‑148 (MQ=255)
                                                                                                           tCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgtc                                         >  1:3718303/1‑148 (MQ=255)
                                                                                                                 gACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAg                                   <  2:191308/148‑1 (MQ=255)
                                                                                                                      tCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTgcg                              >  2:1244696/1‑148 (MQ=255)
                                                                                                                      tCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTgcg                              >  1:2491347/1‑148 (MQ=255)
                                                                                                                        gCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGcc                            <  2:3415828/148‑1 (MQ=255)
                                                                                                                           ctcctGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCcaca                         <  2:1002809/148‑1 (MQ=255)
                                                                                                                           ctcctGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCcaca                         >  1:2881755/1‑148 (MQ=255)
                                                                                                                                gCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGAc                    >  1:1392439/1‑148 (MQ=255)
                                                                                                                                    aaTCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCa                >  2:2166851/1‑148 (MQ=255)
                                                                                                                                     aTCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCAc               >  2:107521/1‑148 (MQ=255)
                                                                                                                                                aTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGACAGCCCCCAGACGTCGCAAAGCTGCGCCACACGGACTGAACCGGAAAGCGTc    >  2:851627/1‑148 (MQ=255)
                                                                                                                                                aTCGATAGAGGTAGGATTACCCGGTGCCACCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTc    <  2:1726000/148‑1 (MQ=255)
                                                                                                                                                aTCGATAGAGGTAGGATTACCCGGTGCCACCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTc    <  2:1883047/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCA  >  CP009273/2912867‑2913160

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: