Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,013 C→A 16.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,0130CA16.1% 492.1 / 54.6 170D365Y (GAT→TAT) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (77/65);  new base A (10/18);  total (87/83)
Fisher's exact test for biased strand distribution p-value = 9.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 172 total aligned reads displayed.

ATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCA  >  CP009273/2912867‑2913160
                                                                                                                                                  |                                                                                                                                                   
ataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCg                                                                                                                                                    <  1:2826683/148‑1 (MQ=255)
 taGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTTATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGa                                                                                                                                                   >  1:327018/1‑148 (MQ=255)
 taGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGa                                                                                                                                                   >  1:327095/1‑148 (MQ=255)
  aGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGAt                                                                                                                                                  <  2:1784165/148‑1 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTCTCGATAg                                                                                                                                                >  2:1147133/1‑148 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAg                                                                                                                                                >  1:2229621/1‑148 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAg                                                                                                                                                >  1:565918/1‑148 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAg                                                                                                                                                <  1:242647/148‑1 (MQ=255)
    gCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAg                                                                                                                                                <  1:2734537/148‑1 (MQ=255)
       tCATGTGCCTTTTGTACCTGGTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAgg                                                                                                                                             <  2:2441587/148‑1 (MQ=255)
             gCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGAt                                                                                                                                       <  1:255192/148‑1 (MQ=255)
             gCCTTATGTACCTGTTCCCAGTCCAGTTACACGCCCAGACCGGGCGCATCAGGAACGCAAATTTTTCCGTTTTTAATATCCAGCGCATTATGGGTCAGGAGACAATCGCCCTCCAGCCAAATCCAGATGGTATCGATACAGGTAGGAt                                                                                                                                       >  2:2591042/1‑148 (MQ=255)
                 tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTAcc                                                                                                                                   >  2:105373/1‑148 (MQ=255)
                 tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTGCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTAcc                                                                                                                                   <  2:1140422/148‑1 (MQ=255)
                     tACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGt                                                                                                                               >  1:994912/1‑148 (MQ=255)
                     tACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGt                                                                                                                               >  1:2918357/1‑148 (MQ=255)
                      aCCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTg                                                                                                                              >  2:1304096/1‑148 (MQ=255)
                       ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgc                                                                                                                             >  2:467189/1‑148 (MQ=255)
                       ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgc                                                                                                                             >  2:469387/1‑148 (MQ=255)
                          gTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGGATAGAGAGAGGTAGGATTACCCGGTgccgc                                                                                                                          >  2:1880405/1‑148 (MQ=255)
                              ccagtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGTCGACAATCGCCCTCCTGCCAAATCAAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCcgcg                                                                                                                      >  2:948857/1‑148 (MQ=255)
                                agtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGcc                                                                                                                    >  1:293238/1‑148 (MQ=255)
                                   ccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCcac                                                                                                                 <  1:2521897/148‑1 (MQ=255)
                                      gtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATg                                                                                                              <  2:1603225/148‑1 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             >  1:862021/1‑148 (MQ=255)
                                       ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                             >  1:1313538/1‑148 (MQ=255)
                                           aCGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaa                                                                                                         <  1:535245/148‑1 (MQ=255)
                                            cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                        >  1:737538/1‑148 (MQ=255)
                                            cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                        >  1:1967485/1‑148 (MQ=255)
                                            cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                        <  2:21489/148‑1 (MQ=255)
                                            cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                        <  2:1698355/148‑1 (MQ=255)
                                               ccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACAt                                                                                                     <  1:1976470/148‑1 (MQ=255)
                                               ccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACAt                                                                                                     <  1:2271246/148‑1 (MQ=255)
                                                 aGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCg                                                                                                   <  1:2940586/148‑1 (MQ=255)
                                                 aGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCg                                                                                                   <  1:2217465/148‑1 (MQ=255)
                                                  gACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGc                                                                                                  >  2:272927/1‑148 (MQ=255)
                                                   aCCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGcc                                                                                                 >  1:2199921/1‑148 (MQ=255)
                                                       ggCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagag                                                                                             <  1:1695061/148‑1 (MQ=255)
                                                         cgcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagagag                                                                                           <  2:2843982/148‑1 (MQ=255)
                                                          gcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAgagaga                                                                                          >  1:1793289/1‑148 (MQ=255)
                                                           cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                         <  2:2767461/148‑1 (MQ=255)
                                                           cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                         <  1:2000611/148‑1 (MQ=255)
                                                                ggAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGa                                                                                    >  2:1186387/1‑148 (MQ=255)
                                                                ggAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGa                                                                                    >  2:641929/1‑148 (MQ=255)
                                                                 gAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGaa                                                                                   >  1:272473/1‑148 (MQ=255)
                                                                    cgtcAATTTGTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGTCGACAATCGCGCTCCTGCAAAATCCTGTTGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATg                                                                                <  1:2591042/145‑1 (MQ=38)
                                                                      gCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGt                                                                              >  2:3023252/1‑148 (MQ=255)
                                                                       cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                             <  2:545810/148‑1 (MQ=255)
                                                                       cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTCAACATCGCCAGAGAGAGATCGAACTGGtt                                                                             >  1:2258238/1‑148 (MQ=255)
                                                                       cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                             <  1:428415/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  2:2674198/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  1:1097999/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  1:445028/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  1:2002082/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  2:3261479/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  2:3260776/148‑1 (MQ=255)
                                                                        aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                            <  2:270628/148‑1 (MQ=255)
                                                                          tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                          <  2:2097108/148‑1 (MQ=255)
                                                                           ttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatta                                                                         >  1:1864111/1‑148 (MQ=255)
                                                                            tttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAg                                                                        <  1:1231976/148‑1 (MQ=255)
                                                                             ttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGa                                                                       <  1:2098076/148‑1 (MQ=255)
                                                                               ccGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAAt                                                                     >  1:2713749/1‑148 (MQ=255)
                                                                                cGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATg                                                                    <  1:1702335/148‑1 (MQ=255)
                                                                                 gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                   <  2:2623584/148‑1 (MQ=255)
                                                                                 gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                   >  2:2186832/1‑148 (MQ=255)
                                                                                 gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                   <  1:2440768/148‑1 (MQ=255)
                                                                                  tttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGc                                                                  >  1:640897/1‑148 (MQ=255)
                                                                                  tttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGc                                                                  <  1:511809/148‑1 (MQ=255)
                                                                                   ttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCa                                                                 >  1:258172/1‑148 (MQ=255)
                                                                                        aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                            >  2:1924711/1‑148 (MQ=255)
                                                                                        aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                            >  1:3089605/1‑148 (MQ=255)
                                                                                          ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                          >  1:1736708/1‑148 (MQ=255)
                                                                                             cAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTc                                                       >  2:923167/1‑148 (MQ=255)
                                                                                               gCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAg                                                     >  2:2408024/1‑148 (MQ=255)
                                                                                               gCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAg                                                     >  1:2001956/1‑148 (MQ=255)
                                                                                                 ggATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcc                                                   <  1:1764640/148‑1 (MQ=255)
                                                                                                 ggATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcc                                                   <  1:1764694/148‑1 (MQ=255)
                                                                                                   aTTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcccc                                                 <  2:293238/148‑1 (MQ=255)
                                                                                                   aTTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcccc                                                 <  2:292275/148‑1 (MQ=255)
                                                                                                   aTTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcccc                                                 <  2:2016845/148‑1 (MQ=255)
                                                                                                    ttttGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTCTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGccccc                                                >  1:1554854/1‑148 (MQ=255)
                                                                                                         ggTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcg                                           <  2:2016340/148‑1 (MQ=255)
                                                                                                          gTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgt                                          >  1:1656929/1‑148 (MQ=255)
                                                                                                           tCAGGCGACAATCGCCCTCCTGCCAAATCCCGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgtc                                         >  2:2907190/1‑148 (MQ=255)
                                                                                                             aGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGc                                       >  2:3068087/1‑148 (MQ=255)
                                                                                                                  aCAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGc                                  <  1:929504/148‑1 (MQ=255)
                                                                                                                  aCAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGc                                  <  1:2584984/148‑1 (MQ=255)
                                                                                                                   cAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCt                                 >  2:773946/1‑148 (MQ=255)
                                                                                                                       cGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTgcgc                             >  2:239336/1‑148 (MQ=255)
                                                                                                                           ctcctGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCcaca                         >  2:241969/1‑148 (MQ=255)
                                                                                                                               tGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                     <  1:3023252/148‑1 (MQ=255)
                                                                                                                               tGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                     <  2:2424990/148‑1 (MQ=255)
                                                                                                                                          gTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGaa          <  2:1683902/148‑1 (MQ=255)
                                                                                                                                          gTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGaa          <  1:1186387/148‑1 (MQ=255)
                                                                                                                                          gTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGAGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGACCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGat          <  2:1683929/148‑2 (MQ=255)
                                                                                                                                              gTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCg      <  1:105373/148‑1 (MQ=255)
                                                                                                                                                aTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTc    >  2:983896/1‑148 (MQ=255)
                                                                                                                                                 tCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTcc   >  2:1859126/1‑148 (MQ=255)
                                                                                                                                                 tCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTcc   <  1:3009825/148‑1 (MQ=255)
                                                                                                                                                  cGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa  >  2:336522/1‑148 (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: