Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,138,006 T→G 5.3% T344P (ACG→CCG)  gldA ← L‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,138,0060TG 5.3% 267.3 / 4.6 89T344P (ACG→CCG) gldAL‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (46/38);  new base G (2/3);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGTTT  >  NC_000913/4137857‑4138155
                                                                                                                                                     |                                                                                                                                                     
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGt                                                                                                                                                       >  2:1259956/1‑150 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAACTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCTACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTc                                                                                                                                                      >  2:1140995/1‑150 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTcg                                                                                                                                                     >  1:248180/1‑150 (MQ=255)
      ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgcc                                                                                                                                                 >  1:997796/1‑150 (MQ=255)
           tAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAG                                                                                                                                            <  2:747507‑M1/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGG                                                                                                                                           <  1:278199‑M1/150‑1 (MQ=255)
             aGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGC                                                                                                                                          <  2:341853‑M1/150‑1 (MQ=255)
              gCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCA                                                                                                                                         <  2:929894‑M1/150‑1 (MQ=255)
                gcagcGCAGCAGGCAATTTTGCGCTCACCCTCCCGCACCACCCGCGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTTTCGTCGCGCCGCCAGGCATG                                                                                                                                       <  2:283446‑M1/148‑1 (MQ=255)
                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAGatcagg                                                                                                                                      <  1:1132587‑M1/150‑7 (MQ=255)
                     gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTtgtg                                                                                                                                  >  1:572899/1‑150 (MQ=255)
                     gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTtgtg                                                                                                                                  <  1:177775/150‑1 (MQ=255)
                         cAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCTGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATg                                                                                                                              >  2:677321/1‑150 (MQ=255)
                         cAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATg                                                                                                                              >  1:1398135/1‑150 (MQ=255)
                          aGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATgg                                                                                                                             >  2:1143848/1‑150 (MQ=255)
                          aGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATgg                                                                                                                             >  2:557601/1‑150 (MQ=255)
                            gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGtt                                                                                                                           <  1:605626/150‑1 (MQ=255)
                              aaTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTc                                                                                                                         <  2:578459/150‑1 (MQ=255)
                                 tttGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCAcc                                                                                                                      <  1:427263/150‑1 (MQ=255)
                                 tttGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCACGCCGCCAGGCATGTTGTGAATGGTTTCAcc                                                                                                                      <  2:1330374/150‑1 (MQ=255)
                                  ttGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCt                                                                                                                     <  1:11252/150‑1 (MQ=255)
                                      gTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTg                                                                                                                 <  2:863953/150‑1 (MQ=255)
                                        tCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCa                                                                                                               <  1:933415‑M2/150‑2 (MQ=255)
                                         cAAACTCCCTGACAAGCCTTGAGTTTGGAGTAGGTTATTCCCACTCTTGAAGGAAACTATGACCGTACTGGTATGCTACCATCAGAGCGGCGTAAAACTTATCTGTCTTCTAGCCTCCAGGCATTTTGTGAATGGTTTCACCTTCTGcac                                                                                                              <  2:933415/150‑1 (MQ=255)
                                           aaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCacat                                                                                                            >  1:1045599‑M2/1‑146 (MQ=255)
                                           aaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGATTCACCTTCTGCacat                                                                                                            <  2:979471‑M2/150‑5 (MQ=255)
                                             cTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCacatgc                                                                                                          <  1:951497‑M2/150‑7 (MQ=255)
                                               cccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCg                                                                                                        <  2:1165347/150‑1 (MQ=255)
                                               cccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCg                                                                                                        >  2:338470/1‑150 (MQ=255)
                                                ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGc                                                                                                       <  1:1215800/150‑1 (MQ=255)
                                                  ggACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTgcagaaggtgaaaccattcacaacatgcctggc                                                                                                     <  1:1510478‑M1/150‑34 (MQ=255)
                                                  ggACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCtt                                                                                                     >  2:246974/1‑150 (MQ=255)
                                                   gACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTgcagaaggtgaaaccattcacaacatgcctggcg                                                                                                    <  1:1141023‑M1/150‑35 (MQ=255)
                                                     cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTgcagaaggtgaaaccattcacaacatgcctggcggc                                                                                                  <  1:491774‑M1/150‑37 (MQ=255)
                                                        gCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCca                                                                                               >  1:109529/1‑150 (MQ=255)
                                                                  tGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAt                                                                                     <  2:997796/150‑1 (MQ=255)
                                                                   ggAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAtt                                                                                    >  1:1060811/1‑150 (MQ=255)
                                                                        aGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGc                                                                               >  1:1133746/1‑150 (MQ=255)
                                                                               tCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACa                                                                        >  1:1108945/1‑150 (MQ=255)
                                                                                cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAt                                                                       >  2:1324624/1‑150 (MQ=255)
                                                                                cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAt                                                                       >  1:1469481/1‑150 (MQ=255)
                                                                                cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAt                                                                       >  1:1127132/1‑150 (MQ=255)
                                                                                cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAt                                                                       >  1:642970/1‑150 (MQ=255)
                                                                                    ctctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAATCGCATTTTCGCCGGGACAtcttc                                                                   >  2:768636/1‑150 (MQ=255)
                                                                                     tctTGCAGGAAACGCTGACCGTCCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAtcttct                                                                  <  2:486380/150‑1 (MQ=255)
                                                                                       ttGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCttt                                                                <  2:178292/150‑1 (MQ=255)
                                                                                        tGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                                                               <  1:733672/150‑1 (MQ=255)
                                                                                        tGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                                                               >  2:279525/1‑150 (MQ=255)
                                                                                         gCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTaa                                                              >  2:582390/1‑150 (MQ=255)
                                                                                            ggAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAatat                                                           >  2:1016656/1‑150 (MQ=255)
                                                                                            ggAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAatat                                                           >  2:276711/1‑150 (MQ=255)
                                                                                            ggAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAatat                                                           <  2:1163406/150‑1 (MQ=255)
                                                                                            ggAAACGCTAACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAatat                                                           >  1:956711/1‑150 (MQ=255)
                                                                                              aaaCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATcc                                                         >  2:778111/1‑150 (MQ=255)
                                                                                                aCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAg                                                       <  1:582390/150‑1 (MQ=255)
                                                                                                 cGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGt                                                      <  1:1233187/150‑1 (MQ=255)
                                                                                                 cGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGt                                                      >  1:605384/1‑150 (MQ=255)
                                                                                                    tGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTGGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGa                                                   >  1:316340/1‑150 (MQ=255)
                                                                                                     gACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAg                                                  <  1:1016656/150‑1 (MQ=255)
                                                                                                     gACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAg                                                  >  1:364050/1‑150 (MQ=255)
                                                                                                       ccGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGc                                                                                       >  2:229998/1‑111 (MQ=255)
                                                                                                       ccGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGc                                                                                       <  1:229998/111‑1 (MQ=255)
                                                                                                       ccGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCg                                                <  2:1127132/150‑1 (MQ=255)
                                                                                                         gTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag                                              <  2:1398135/150‑1 (MQ=255)
                                                                                                         gTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag                                              <  2:1323579/150‑1 (MQ=255)
                                                                                                          tACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgaga                                             >  2:1104513/1‑150 (MQ=255)
                                                                                                           aCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGagag                                            <  1:422639/150‑1 (MQ=255)
                                                                                                            cTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGt                                           >  2:367000/1‑150 (MQ=255)
                                                                                                            cTGGGCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGt                                           <  2:20402/150‑1 (MQ=255)
                                                                                                              ggTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTa                                         >  2:831629/1‑150 (MQ=255)
                                                                                                              ggTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTa                                         >  2:1085273/1‑150 (MQ=255)
                                                                                                              ggTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAGGTTCTCAATGGTTTCACCTTGTGCACATGCCGCATCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTATAATATCCAGTTGAGCGAGAGTTa                                         >  1:535034/1‑150 (MQ=255)
                                                                                                                   gCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGc                                    >  1:139320/1‑150 (MQ=255)
                                                                                                                     tACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCaa                                  >  2:505028/1‑150 (MQ=255)
                                                                                                                      aCCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCaaa                                 >  2:514414/1‑150 (MQ=255)
                                                                                                                      aCCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCaaa                                 >  1:637983/1‑150 (MQ=255)
                                                                                                                        cagcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAAcc                               >  1:132302/1‑150 (MQ=255)
                                                                                                                        cagcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAAcc                               >  1:1353111/1‑150 (MQ=255)
                                                                                                                        cagcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAAcc                               <  2:830831/150‑1 (MQ=255)
                                                                                                                         agcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCt                              >  2:94125/1‑150 (MQ=255)
                                                                                                                                          catgccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATgg                     <  2:340269/137‑1 (MQ=255)
                                                                                                                                          catgccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATgg                     >  1:340269/6‑142 (MQ=255)
                                                                                                                                          cTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGc             >  1:1318494/1‑150 (MQ=255)
                                                                                                                                           tGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCa            >  1:1531753/1‑150 (MQ=255)
                                                                                                                                           tGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCa            >  1:41237/1‑150 (MQ=255)
                                                                                                                                            gATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCGCAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAg           >  2:1009029/1‑150 (MQ=255)
                                                                                                                                             gccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGc          >  2:531527/3‑150 (MQ=255)
                                                                                                                                             aTCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGc          <  1:831629/150‑1 (MQ=255)
                                                                                                                                             aTCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGc          <  1:1143848/150‑1 (MQ=255)
                                                                                                                                                tGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTAc       <  1:1104513/150‑1 (MQ=255)
                                                                                                                                                     tCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGttt  >  1:491781/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGTTT  >  NC_000913/4137857‑4138155

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