Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,127,934 T→G 6.5% T344P (ACG→CCG)  gldA ← glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,127,9340TG 6.5% 186.5 / 3.3 72T344P (ACG→CCG) gldAglycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base G (4/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCT  >  CP009273/4127788‑4128076
                                                                                                                                                  |                                                                                                                                              
aGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTcg                                                                                                                                              >  1:453964/1‑149 (MQ=255)
 gcccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTcgc                                                                                                                                             <  2:1157575/147‑1 (MQ=255)
   ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCgcgc                                                                                                                                           >  1:419921/1‑149 (MQ=255)
      gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccgc                                                                                                                                        <  2:416174/149‑1 (MQ=255)
            cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCa                                                                                                                                  >  1:269085/1‑149 (MQ=255)
            cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatc                                                                                                                                  <  1:424915/149‑4 (MQ=255)
               aGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                                                                                                               >  2:557407/1‑149 (MQ=255)
                gcgcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgtt                                                                                                                              >  2:788246/1‑149 (MQ=255)
                 cgcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGCGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAcgttg                                                                                                                             >  2:117971/1‑149 (MQ=255)
                     tCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCCGGCATGTTGTGaa                                                                                                                         >  2:325704/1‑149 (MQ=255)
                     tCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGaa                                                                                                                         <  1:946095/149‑1 (MQ=255)
                     tCAGGCAATTTGGCGTTCACACTCCCGAGCAAGCCGGGAGTTTGGAGTAGGTTACTCCCCCCCTTGCAGGAAAGCCTGACCGTACCGGGTGGCTACCCGCAGTGCGGCGTAAACCCGATCTGGCGTCGCGCCGCCAGGCATGTTGTGaa                                                                                                                         <  2:52211/149‑1 (MQ=255)
                           aaTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGttt                                                                                                                   <  1:126860/149‑1 (MQ=255)
                           aaTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGttt                                                                                                                   <  2:214381/149‑1 (MQ=255)
                              tttGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTGGGTTATTCTCACTCTTTCAGGAAACGCTGGCCGTAGTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCACGCCCGTTGTGGATGGtttagc                                                                                                                >  2:812763/1‑146 (MQ=255)
                                 gCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCtt                                                                                                             <  1:248510/149‑1 (MQ=255)
                                    ttCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGAGAGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGCCGCCAGGCAGGTTGTGAATGGTTTCACCTTCTg                                                                                                          >  1:36985/1‑149 (MQ=255)
                                    ttCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTg                                                                                                          >  1:554315/1‑149 (MQ=255)
                                      gaaccTCCCGTACAAGCCGGGAGTTTGGAGTAGGTTATTCCCCCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGCTCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGca                                                                                                        <  2:866363/145‑1 (MQ=255)
                                      caacgccccGGACACGCCGGGAGTTGGGAGTAGGTTATTCCCACTCGTGCAGGAAACGCTGACCGTACTGGTCGGCTCCCAGCCGAGCGGCGTAAACCTGGTCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCGCCATCTGca                                                                                                        <  2:13061/143‑1 (MQ=255)
                                      cAAACTCCCGGACAAGCCGGGAGTTTGGGGTAGGTTATTCCCACTCTTGCAGGAAACGCTGCCCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGca                                                                                                        <  2:375330/149‑1 (MQ=255)
                                      cAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGca                                                                                                        <  1:189127/149‑1 (MQ=255)
                                      cAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGCCATGTTGTGAATGGTTTCACCTTCTCca                                                                                                        <  1:250857/149‑1 (MQ=255)
                                           tCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc                                                                                                   <  2:348680/149‑1 (MQ=255)
                                           tCCCGGACAAGCCGGGAGTTTAGATTAAGTTATTCCCACCCTTGCAGGAAGCGCCGACCGTCCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc                                                                                                   <  2:199105/149‑1 (MQ=255)
                                           ccccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACAGGTCGGCTTTCTGCCGCGATTTATAATCCTCATTTGGCGACGAGCCTCCAGGCATGTTGTGACTTCCTTTTCCTTCTGCACATGc                                                                                                   <  2:570085/148‑1 (MQ=255)
                                             ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTAGTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCGGCACATGGCg                                                                                                 >  2:709891/1‑149 (MQ=255)
                                             ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCg                                                                                                 >  2:68358/1‑149 (MQ=255)
                                              cGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCGGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGc                                                                                                >  1:683454/1‑149 (MQ=255)
                                                gACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGCTCTGGCGTGGCGCCGCCAGGCATGTGGTGAAGGGATTCACCTTCTGCACATGCCGGtt                                                                                              >  2:660398/1‑149 (MQ=255)
                                                 aCAAGCCGGGAGTTTGGAGTAGGTTATTCCCAATCTGGTCGTAATCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTc                                                                                             >  2:1134508/1‑149 (MQ=255)
                                                   aaGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATATGGCGTCGCGCCGCCAGGCATGTTGTCACTGGCTTCTCCTTGTGCACATGCCGCTtgtg                                                                                           >  2:618868/1‑149 (MQ=255)
                                                    aGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGc                                                                                          <  1:978920/149‑1 (MQ=255)
                                                    aGCCGGGAGGTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGc                                                                                          <  2:275825/149‑1 (MQ=255)
                                                     gCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGATACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTGCTGcc                                                                                         >  2:573739/1‑149 (MQ=255)
                                                              ttGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGc                                                                                <  1:1054741/149‑1 (MQ=255)
                                                              ttGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGc                                                                                <  1:884011/149‑1 (MQ=255)
                                                                ggAGTAGGTTATTCCCACTCTTGCAGGAAACGCTTACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAt                                                                              <  2:60186/149‑1 (MQ=255)
                                                                 gAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAtt                                                                             >  1:994969/1‑149 (MQ=255)
                                                                    tAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTc                                                                          >  1:200615/1‑149 (MQ=255)
                                                                    tAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTc                                                                          >  1:455116/1‑149 (MQ=255)
                                                                     aGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGATTTGAAAGTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCg                                                                         >  1:648286/1‑149 (MQ=255)
                                                                     aGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAGGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTCCTGCCACAATTCGCATTTTCg                                                                         >  2:568813/1‑149 (MQ=255)
                                                                        ttATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCCTTTTCGCCg                                                                      >  2:543434/1‑149 (MQ=255)
                                                                        ttATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCg                                                                      <  2:947793/149‑1 (MQ=255)
                                                                        gtatccgcACTGTTGCGCGAAACCGTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAATCCTGATCTGGCGTCGCGCCGCCAGGCATGCTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCg                                                                      <  2:122695/142‑1 (MQ=255)
                                                                          aTTCCCCCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGTCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACATTTCGCATTTTCGCCggg                                                                    <  2:917268/149‑1 (MQ=255)
                                                                            tCCCACTCTTGCAGGAAACGCTGACCGTACTGGGCGGCTACCAGCAGAGCGGCGTAACCATGATCTAGCGTCGCGCCGCCAGGCAGGTTGTGAATGGTTTCACCTTCTGCACATGCCGCGTCTGCCACAATTCGCATTTTCGCCGGGAc                                                                  >  1:908056/1‑149 (MQ=255)
                                                                             cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACa                                                                 >  1:217631/1‑149 (MQ=255)
                                                                                aCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAgct                                                              >  2:1004023/1‑149 (MQ=255)
                                                                                aCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATATGGCGGAACGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACCTGCCGCTTCTGCCACAAGTCGCATTTTGCCCGGGTCAtct                                                              >  2:510464/1‑149 (MQ=255)
                                                                                    ttGCAGGAAACGCTGCCCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATcttctt                                                          <  1:155832/149‑1 (MQ=255)
                                                                                        aGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCGGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAat                                                      >  1:756010/1‑149 (MQ=255)
                                                                                           aaaCGCTGACCGTACTGGTCGGCTACCAACAGAGCGGCGTAACCCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTTCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATc                                                   <  1:107038/149‑1 (MQ=255)
                                                                                              cgcgGACCGTACGGGTGGGCTACCGGGAGGGCGGCGTAAACCTGATCGGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAg                                                <  2:217876/149‑1 (MQ=255)
                                                                                              cGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAg                                                <  2:1038305/149‑1 (MQ=255)
                                                                                                  gACCGTACTGGTCGGCTACCAGCAGGCCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGGCGCTTCTGACACAATTCGCATTTTCGCCGGGACATCTTCTTTAATACCCAGTTGa                                            <  2:634567/149‑1 (MQ=255)
                                                                                                    ccGTACTGGTCGGCTACCAGCAGCGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCCCCTTCTTCACATGCCGCTTCTGCCACAATCCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGc                                          <  1:962482/149‑1 (MQ=255)
                                                                                                       tACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag                                       >  1:525839/1‑149 (MQ=255)
                                                                                                               ggCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGGACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTg                               <  1:557407/149‑1 (MQ=255)
                                                                                                                gCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTgg                              <  1:188284/149‑1 (MQ=255)
                                                                                                                  tACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCa                            >  1:821381/1‑149 (MQ=255)
                                                                                                                  tACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCa                            >  1:234167/1‑149 (MQ=255)
                                                                                                                    ccAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTCATATCCAGTTGAGCGAGGGTTATTGGCaaa                          >  2:302590/1‑149 (MQ=255)
                                                                                                                         agagCGGCGTAAACCTGGTCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTAc                     <  2:195697/149‑1 (MQ=255)
                                                                                                                               gcgTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATTGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGGGAGGGTTATTGGCAAACCTACCGCATg               >  2:596535/1‑149 (MQ=255)
                                                                                                                                      ccTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATCTCCAGTTAAGCGGGTGTTATTGGCAAACCTACCGCCTGGCTAAgg        >  1:538361/1‑149 (MQ=255)
                                                                                                                                       cTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAggg       >  1:789312/1‑149 (MQ=255)
                                                                                                                                       cTGATCTGGCGGCGCGCGGCCAGGCATGGTGGGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAggg       <  2:32499/149‑1 (MQ=255)
                                                                                                                                          gccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAg    >  2:849717/3‑149 (MQ=255)
                                                                                                                                           tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCCCCGTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGACAGc   <  1:29666/149‑1 (MQ=255)
                                                                                                                                            cTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCt  >  1:1044966/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCT  >  CP009273/4127788‑4128076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑